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a - Inverse slope of load waterline, i.e. tangent of the angle 

between the centerline plane and the load waterline - tan V 
a 0 - Average distance between ice plates in the ice plate-brine 
pocket model 

Bo - Constant depicting ice location in conversion of air temperature 

to ice surface temperature 

A - Non-dimensional coefficient 

A = / £b * fcl 8 C b ) 

{ Cyj 4k 2 2 C W ' 2 / 

A„ - Apparent contact area between hull steel and ice surface - sq. ft. 

a 

A*. - Real contact area between hull steel and ice surface - sq. ft. 

At - Area of contact between hull end ice - sq. ft. 

^w - Load waterplane area - sq. ft. 

b - Assumed area of loading when considering stresses in a flat plate 
under a concentrated load - sq. inches 
b^ - Temperature compensating constant 
b 2 “ Width of a wedge at any point x - inches 
b 0 - Average spacing of brine pockets 
b© - Width of a wedge at X ; 1; i.e., describes angle 
b r - Brine content of sea ice 
B - Maximum beam at load waterline - ft. 

Bj_ - Factor related by poissons ratio 
Bg - Faotor related by poissons ratio 
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BKt - Height of transverse metecenter above center of buoyancy - ft. 

BTF - Blade thickness fraction 

c - Radius of circular area over which load P distributed - sq. inches 
CG - Center of Gravity 

CB - Center of Buoyancy 

Cb - Block Coefficient s Cp i C, 

C 0 - Coefficient of Dfag. 

Cj^ - Transverse inertia coefficient 

CL - Coefficient of Lift 

CiL ” Longitudinal inertia coefficient 

Cp - Prismatic coefficient 

C v - Vertical prismatic coefficient = c »/c« 

G» - Waterplane coefficient 

C x - Midship Section coefficient 
D - Propeller diameter - ft. 

3 

D^ - Flexural rigidity of a flat plate z ^ /12(1 - H^) 

D 0 - Depth Amidships - ft. 

£ - Drag - lift ratio s ^o/cl 

e - Coefficient of restitution 

g - YovmCt'j ModwIuS * e^ L 

EHP - Effective Horsepower - Hp 

f - Frequency of v/quat/ov 

f 0 - Dynamic coefficient of friction = ^ 0 
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Statio coefficient of Friction = 

Vertical bow reaction force - tons 
Vertical bow reaction force - lbs. 

Average area of a brine pocket in surface plane 
Friction force - lbs. 

Force tending to move vessel through ice - lbs. equals T^ 
in steady state when inertia energy dissipitated 
Acceleration of gravity - ft/sec 2 
Length of brine cylinders 
Average spacing of brine pockets 
Transverse metecentric height - ft. 

Longitudinal metecentric height - ft. 

Thickness of ice sheet or plate - cm or inches equals t" 
Height of vertical forefoot - ft. 

Snow thickness - cm or inches 
Vessel full load draft - ft. 

Slant distance vessel penetrates ice due to compression, 
i.e., length of load waterline to penetrate ice - ft. 
Speed Coefficient - ^ 

Foundation Modulus - lb/in 
Inertia coefficient 

Nondimensional coefficient defined by S s k ,^/Z 
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- Shape factor of loaded area 

- Distance from base line to vertical center of buoyancy - ft. 

- Distance from base line to vertical center of gravity - ft. 

- Distance from base line to transverse metecenter - ft. 

- Distance from base line to longitudinal metecenter - ft. 
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Characteristic length z) 12k 
Length between perpendiculars - ft. s LBP 
Length overall - ft. 

Longitudinal center of floatation 
Load waterline length - ft. 
metecenter 

o 

added mass - lb/sec c /ft. 
vessel mass - lb/sec 2 /ft. 

Bending moment at upper fibers of wedge - inch lb. ’ . 

Distance from amidships to longitudinal center of floatation 
Propeller shaft HPM 
Normal load - lbs. 

Distance from forward perpendiculars to vertical portion of 
forefoot - ft. 

Plastic mean flow pressure - psi 
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- propeller pitch at the 0.7 radius - ft. 

P - Propeller Horsepower = 0.96 SHP « PHP - Hp 
P - Applied concentrated load - lbs. 
q - Uniformly distributed load on wedge - lb/inches 
Q, - Propeller torque - ft. lbs. 

Q,1 - Compressive force against the ships side at the bow ezerted by 
ice - lbs. 

■ Component of q! normal to the shell - lbs. 

- Component of Q,1 tangent to the shell - lbs. 
r - Ship resistance at any speed v - lbs. 

r^ - - Total distance vessel travels in horizontal direction - ft. 

r c - radius of assumed concentrated loadings of wedge 
r b - one half width of included brine pockets 
R t - Resultant force acting normal to the stem - lbs. 
s - Distance from intersection of stem end load waterline to the 
longitudinal center of floatation - ft. 
si - Vessel accelerating distance - ft. 
s b - Salinity of brine 

s^ - Relative insulation value of snow as versus ice 
b 0 - Salinity of ice 
s s - Relative amounts of solid salts 
S . - Shape factor of loaded area 
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SHP - Shaft Horsepower - Hp 

t - Thrust deduction factor 

t 1 - M /P3t 

t 11 - Thickness of ice sheet - cm or inches = h 

T - Propeller thrust - tons 

- Propeller thrust - lbs. 

Ti - Ideal propeller thrust 

T 0 - Propeller thrust at b ->****o condition - tons 

T 0 1 - Propeller thrust at conditions - lbs. 

u - tangent of hydrodynamic pitch angle = tan /d/ 

U - Non dimensional coefficient 

TJ f 8 - tan -< cos j 
f g tan ^ \ cos j 

y - Vessel speed - ft/sec. 

v^- - Relative brine volume 

v a - Propeller Speed of Advance - ft/sec. - v (1 - w) 

v e - Extensional velocity of sound wave through infinite solid medium - 
ft/sec. 

v 0 - Longitudinal velocity of sound wave in infinitely long bar - ft/sec. 

v r - Transverse velocity of sound wave over surface of extended solid 

medium - ft/sec. 

V - Vessel speed - knots 

w - Wake fraction 

W - Vessel displacement - tons 
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distance along wedge from apex 



Non-dimensional coefficient 



X = 



cot <*. cos /3 - / p 
cos 3 f f D cot 



nondimensional dub diameter 



hub diameter 

propeller diameter 



Vessel impact speed variable 

Horizontal vessel travel due to ice compression - ft. 

Distance from edge of flat plate to the centroid of loaded area 

Non-dimensional coefficient 

Y - cos 

cos 3 ♦ f;> cot « 

Non dimensional coefficient 

cos 3 cos 2 ■ _ y cos 2 * 

COS i f COt.x 

angular vi® of forefoot - degrees 

Angle between center line plane and the normal to the shell plate 



at the bow - degrees 

tan“ x /jtno 

Non-dimensional parameter = 

% 

Hydrodynamic pitch angle 

Angle measured in the plane normal to the ships side at the bow, 
between the centerline plane and the ships side - degrees 
Specific weight of sea water lb/ft^. 

Non dimensional parameter s ®Vgo 

Average ratio of length of brine pocket to the width of brine pocket 
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- Hull efficiency s " V 1 - w) 

- Mechanical efficiency » 0.98 

- Blade efficiency 

- Ideal propeller efficiency 

- Open water propeller efficiency 

- Relative rotative efficiency 

- Propeller coefficient 

^ ;r — P D S 
K 6J / - H 

- Dynamic coefficient of friction - f D 

- Static coefficient of friction s f^ 

- Poisson’s ratio 

- Density - lb. sec 2 /ft. 

- Non dimensional parameter 
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- Propeller coefficient 
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- Tensile strength of sea ice - psi 

- L~ p i *>/c AlL f determined ice strength in the region of ice deterioration 
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- Compressive strength of ice - psi 

- Tensile strength of fresh water ice - psi 1 

- Tensile strength of ice without internal cavities - psi 

- Bending stress in upper fibers of wedge - psi 
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Tensile strength of material - psi for ice 
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Non dimensional distance of loading 

^ a 
l = — , 

u 

shear strength - psi 

X 

Non dimensional distance measured along wedge - < 

Angle between the centerline plane and the load waterline, i.e. 
a r tan ^ 

Function representing the reduction in strength of ice due to 
brine inclusions 

Angle between the centerline plane and the ships side at a 
transverse bow section 
Change in vessel draft - ft. 

Change in vessel trim - radians 
Air temperature - °C 

Temperature at the interface between ice and snow - °C 
Water temperature - °C 
Propeller Coefficient 

y ") _ — J l J 

^ ^ i , , 

Volume of water displaced - cu. ft. 

Advance coefficient • tan > » ~~- 
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Part 1 
Introduction 

The true icebreaker, as distinguished from the vessel strengthened 
for ice, is a unique type of sea going vessel. It must be capable of 
continuous operation in areas of the world seldom seen by the conventional 
vessel and under environmental conditions which are often extremely severe. 
The vessel must be strong and powerful so as to be able to withstand and 
overcome its element - ice, and it must be self sufficient so as to be 
capable of long sustained operation far from any possible source of supply 
or material assistance. 

Primarily, icebreakers are vessels designed to force passage throu^i 
ice filled waters so as to accomplish several allied functions. One of the 
most important of these functions is that of escorting or leading conven- 
tional thin skinned vessels through the ice field by breaking a channel 
in the ice wide enough for the escorted vessels to pass through. Another 
of these functions is that of keeping a given waterway or channel free of 
solid ice to such an extent that conventional vessels may traverse the 
waterway without danger. The icebreaker is often called upon to free 
vessels beset by ice in the various waterways, rivers or lakes and at 
times even in polar regions. Lastly, the icebreaker must be capable of 
long periods of independent operation in polar regions for the purpose of 



carrying out valuable hydrographic or oceanographic survey and also to 
provide logistic support to remote ice-locked stations. 

The history of ice breakers is comparatively short -when compared with 
the history of navigation. The first true icebreaker was probably 
•EISBRE.CHER I* which was built in Hamburg in 1871 and worked between Hamburg 
and Cuxhagen. It had a displacement of 500 metric tons and an engine power 
of 600 IHP. Since that time many changes and innovations have occurred in 
icebreaker design practice until today there exist such goliaths as the 
8000 ton, 21000 SHP "Glacier** in the United States and the 16000 ton, 

40000 SHP atomio powered "Lenin" in Russia. 

Many new theories as regards icebreaking were developed during this 
evolutionary period. The spoon shaped bow of the "MURTAJA" in 1890 was a 
revolutionary change. Its effectiveness in solid ice was gratifying, 
however, in snow covered ice or ice slush it proved a failure since the 
spoon shape tended to push the snow and slush forward with the bow quickly 
absorbing the vessels forward energy. At about the same time, on the 
great lakes, the effectiveness of forcing a passage through ice by backing 
the vessel through it was observed. The adaption of bow propellers was a 
natural consequence, beginning in 1888 with the ferry boat "ST. IGNACE" 
with a stern shaft power of 2000 IHP and a forward shaft power of 1000 
IHP. The primary function of this bow propeller was to wash water and 
ice from the forward end of the vessel thereby reducing the friction between 
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the ice and the vessel. Unfortunately, the wash of the forward propeller 
was unsymetrical due to direction of propeller rotation and the single 
forward propeller always produced a side force which made steering and 
maneuvering to one of the sides more difficult. On many subsequent in- 
stallations this condition was alleviated by the use of twin forward screws. 
The use of these bow propellers proved very effective in bays, rivers and 
lakes. In later yeers, however, when they were employed on polar ice- 
breakers serious shortcomings became evident. In these polar regions 
where the ice is very thick and hard and breaking by ramming was a 
necessity, considerable propeller damage resulted. Consequently, bow 
propellers are seldom used on present arctic breakers. 

In 1926 oil fired boilers began to replace the previously used coal 
fired installations, thus greatly increasing the icebreakers cruising ' 

range. At about the same time heeling tanks were introduced which could 

fi*- 

heel the vessel several degrees by means of pumping water between opposing 
wing tanks, thus helping to free a vessel beset by ice friction. 1932 
saw the introduction of direct current diesel electric propulsion with 
its inherent high SHP/weight ratio and greatly increased vessel maneuvera- 
bility through feasable pilot house engine control. 

Despite these various evolutionary changes in icebreaker design, it 
is also true that basic icebreaker hull form has changed little through 
the years. Vessel displacements and powering have been increased and 
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propeller shaft configuration and number have been varied but through 

it all the basic hull shape and configuration -has remained about the seme. 

In fact, to the present day, little, if any, factual test or experimental 

data exists as to the effect of variation of hull form on ice breaking 

effectiveness. Limited personal experiences as to an existing vessel's 

effectiveness have been about the only guide available to the designer. 

Consequently, each new icebreaker is no more than a scaled up model of 

some existing vessel which it is f.elt has proven effective. No doubt the 

present hull form is effective, but is it as effective as it could be 

for its given displacement and power 7 It is amazing that even today, 

little, if any, theoretical evaluations or actual test data concerning 

icebreaking and ice breaking vessels is to be found in the literature. 

Has no one ever considered changing the present hull form in search of 

a more effective one^ perhaps a hull form which broke ice upward which 

* 

is considerably easier to do that downward, or a hull form that attempted 
to conserve the major portion of potential energy developed by ramming 
so that a quasi-constant vertical bow reaction force is maintained against 
the ice surface and only a small additional propeller thrust is necessary 
to break even heavy ice in a continuous forward motion. This latter method 
would at least keep the vessel's bow from falling back down when the ice 
broke after each ram and would be efficient which cannot be said of the 
present hull shape end method. Perhaps these things and others have been 
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thought of but since comparatively few icebreakers are built there is 
little incentive for publication. Since World War II, however, interest 
in the polar regions has materially increased. More nations of the 
world are finding that they have increasing interests and commitments of 
a military, commercial and scientific nature in these areas. Yet the 
majority of existing icebreaker type vessels are obsolete and unable to 
successfully cope with the present demands of operations in the Arctic 
and Antarctic. New, more powerful and more effective icebreaking vessels 
will, no doubt, be built. To do this most effectively it would seem that 
a more positive program of theoretical analysis, research and test is 
necessary both in the field of ice mechanics and >ct6& efitez. Naval Architec- 
ture. 

This paper attempts to present a discussion of the mechanics of ice 
breaking plus a means of approaching the preliminary design or design 
study of the conventional type of icebreaking vessel. It is only one of 
the many possible approaches and by no means purported to be the best, 
most logical or most correct possible. However, it is an approach and 
one which attempts to make effective use of what little is known about 
ice and ice breaking. Throughout the paper many assumptions are made 
and conclusions drawn. To the author they seemed reasonable and justifiable. 
Others might tend to disagree or object violently. However, to lend 
authority or credence to an objection or opposing view, justification must 



be presented to substantiate the objection. If this paper accomplishes 
even this end it shall have, to a degree, achieved its purpose, for only 
through published critical discussion and differences of opinion will the 
elements of icebreaker design be critically re-evaluated so as to stimulate 
a real desire towards determining the most effective means of breaking ice. 

There axe many types of icebreakers and each is designed so that it 
may best cope with the ice and environmental conditions of its particular 
service. It would be impossible for this paper to treat the particular 
problems of each of the various icebreaker types. Consequently, to bring 
the problem within reasonable limits, we shall consider the case of the 
large polar icebreaker (arctic or antarctic) with no bow propellersand with 
two stern propellers. Invariably, many of the items discussed herein will 
apply with equal validity to any type breaker and one with any powering 
configuration, however, in those instances where definition must be made, 
reference will always be to the twin stern screw polar icebreaker. It 
should also be noted that throu^icut the paper, full load displacement 
condition is used synonomously with design conditions and that when reference 
is made to vessel displacement, draft, beam, etc., it is always to the 
full load or design condition unless specifically noted to the contrary. 

Since only the twin stern screw vessel is considered, all references to 
vessel length are to length between perpendiculars. One final comment 
is in order. Snow cover on an ice surface will greatly effect the physical 



and elastic properties of that ice. Since so little is actually known 
quantitatively about the effects of the snow cover all references to ice, 
its strength and properties, is in every case to ice with no snow cover. 
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Pert 2 

Icebreaker Preliminary Design Considerations 

In order to commence a logical preliminary design of an arctic ice- 
breaker, there are certain items of information concerning the proposed 
vessel, its service and general area of normal operation which must be 
known to the designer. These items of information need not be known with 
great accuracy for the ship will invariably be a compromise as regards its 

service, area of operation and characteristics. However, in order to allow 

% 

for as realistic and accurate preliminary vessel parameter estimates as is 
possible, it behooves the designer to acquaint himself as thoroughly as is 
possible with these requirements which will either restrict or control his 
selection of vessel parameters. 

Generally speaking, those items which will have the most immediate 
effect upon the design are: 

(a) Physical size limitations imposed on the vessel 

1. Draft restrictions 

2. Minimum beam considerations 

(b) Normal operating area 

1. Arctic 

2. Antarctic 

(c) Operational requirements 

1. Convoying ability (high power and high length to beam ratio) 
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2. Independent operation (low length to beam ratio) 

3. Maximum thickness of ice which it is desired the 

vessel be able to overcome. 

4. Month of year during which it is desired vessel be 

r \HiT 

Jsst able to successfully penetrate ice. 

5. Mean air temperature typical of area of operation 

during period described in 4 above. 

The desired draft of the proposed vessel must be considered early 
in the design stage. Often there are draft limitations due to the restricted 
depths of water found in the vessel’s proposed area of operation. If no 
such draft restrictions exist then the draft should be selected large enough 
to accomodate large diameter propellers. This is desirable so as to develop 
maximum possible thrust yet still allow for sufficient propeller tip immersion 
so as to prevent the propeller tips from striking the bottom surface of a 
floating ice sheet and in addition, prevent air from being drawn down by the 
propeller. Unfortunately, however, most areas of operation in the arctic 
and antarctic impose some draft limitations. Consequently, a serious study 
of the vessels most probable area of operation should be made and the maximum 
safe operating draft determined. In most instances, for the arctic breaker, 
this limit seldom exceeds 29 to 30 feet. 

The minimum desirable extreme beam at the full load waterline is 
another item which must be considered early in the design. If the proposed 
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vessel is to be capable of efficient convoy service in the ice, it is often 
desired that the extreme vessel beam be equal to or greater than the beam 
of the largest vessel which the icebreaker may have to escort through pack 
ice. The channel that an icebreaker leaves is often about a foot or two 
feet wider than the vessel itself, dependent upon the thickness and hardness 
of the ice encountered. Hence, some decision as to the minimum acceptable 
extreme beam must be made. Such a vessel will normally have a relatively low 
length to beam ratio, in the range of about 4 to 4.8. These values are rather 
low and tend to produce a vessel with either poor or mediocre directional 
stability. In order to escort vessels through the ice effectively and safely 
it is necessary for the escorting icebreaker to have good directional stability 
so as to be able to maintain as steady and straight a course through the ice 
as is possible. Such a requirement would dictate a rather high length to beam 
ratio, probably around 5.5 or 6. Unfortunately, however, this would restrict 
the permissable extreme beam. Consequently a compromise choice of beam must 
be made. Some authorities on icebreaking feel, and the author tends to agree, 
that the minimum extreme beam should not be such an overiding consideration. 

It is felt that in convoying it is much more vital to have adequate directional 
stability so that the vessel will not tend to veer off sharply when the bow 
strikes a piece of ic« a glancing blow, hence producing a channel in the ice 
which is straighter and easier for the escorted ships to follow. Due to this 
high length to beam ratio, however, the extreme beam may be limited to a value 
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less than the beam of the escorted vessels. However, since most vessel 
convoying in ice is done during the summer months, the ice will normally 

V 

be open, with a considerable number of leads and not under pressure. Under 
these circumstances, -«'as the icebreaker forces a lead or traverses a floe, the 
ice will be able to move aside giving way to the icebreaker and leaving a 
resultant channel astern the breaker which is considerably wider than the vessels 
extreme beam. If, perchance, the ice is not open, but rather concentrated 
with almost ten tenths coverage and possibly still under pressure, then 
practical experience has shown that it is best to use two icebreakers to 
escort,. "^he first icebreaker with good directional stability forcing the 
channel and the second trimming off the edges of the channel so as to insure 
adequate width of channel and straightening any irregularities caused by passage 
of the first breaker. There are considerable pro and cons on this topic of 
length to beam ratio and minimum extreme beam which is acceptable in a good 
convoying icebreaker. It is felt, however, that the tide of opinion is 
changing and that more and more people involved in icebreaking are becoming 
convinced that directional stability is of paramount importance in an arctio 
icebreaker, hence the proponents of higher length to beam ratios are increasing. 

If the escorting of vessels through the ice is not to be a primary task 
of the proposed icebreaker, then lower length to beam ratios are acceptable 
and in fact, desirable. This will produce a short stubby vessel capable of 
quick sharp turns which can wind itself through the ice, taking advantage of 



- 3 t - 



every lead in Its forward progress. Such a vessel will have poor directional 
stability. Maneuverability, however, is much more important when forcing 
independent passage through heavy ice so that complete advantage may be taken 
of leads, polynays and the evasion of areas of rafted or heavily snow covered 
ice. * 

There is no doubt that the selection of length to beam ratio and minimum 
extreme beam for an icebreaker vessel is a compromise. Consequently, it behooves 
the designer to understand, as completely as is possible, the needs of the 
proposed vessel before making the selection of values. 

In addition to the consideration of vessel convoying in ice, the beam 

must be sufficient to insure adequate intact and damaged stability. As a 

result of their normally large beams, this type vessel has a large transverse 

metecentric height. Consequently, intact stability is normally very good -so 

good in fact as to often prove a discomfort to the crew. This stability, 

however, must not be assumed to be adequate merely because of a high initial 

metecentric heignt. Icebreaker operations are often such as to make such a 

%ol£L¥ 

high metecentric height desirable and necessary and not se^ a y the by product 

of a beamy vessel. Most characteristic is the serious reduction in stability 

which occurs when the vessel rides up onto the ice, especially v.ith iced top 
h 

One other item which should be seriously considered prior to the commence- 
ment of an icebreaker design study is the probable area of normal vessel opera- 
tion and the periods during which the vessel will actually be operating in ice. 
Environmental conditions differ in the arctic as compared to the antarctic and 
vessels which can operate with success during normal periods in the arctic 
often experience serious difficulty in the antarctic. The month of year 
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daring which the icebreaker is first expected to be able to penetrate the 
ice will materially effect its requirements; the thickness of ice normally 
found in the area, past records of ice growth rates, history of ico pressure 
during various months of the year, rafting, snow cover experienced, etc., 
all should be considered so as to more realistically determine the probable 
needs of the vessel. Considerable study of ice conditions both in the arctic 
and antarctic have been made for a number of years by such organizations as 
the Hydrographic Office, U. S. Navy Civil Engineering Corps., U. S, Army Snow, 

Ice and Permafrost Research Establishment, Navy Electronics Laboratory, etc.. 

Each of these organizations has e wealth of ice data and information concerning 
past ice experience at various times of the year and under various operating 
conditions which is readily available to the designer. Consequently, if behooves 
the designer of an icebreaking vessel to avail himself of as much of this data- 
and past history as is possible so that he can more realistically realize the 
actual needs of the proposed vessel. 

In this paper we attempt to make direct use of several of these en- 
vironmental factors which can be reasonable predicted from past records. Among 
the items we consider are mean air temperature during the month of first 
expected ice penetration and probable snow cover during this time. Based on 
this data we can estimate the probable ice surface temperature which might be 
expected. Other items we consider are noraial ice thickness encountered in the 
past in the general area of operation and at the desired times of the year and 
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the type of pack ice experienced prior to natural weather breakup, i.e., 
perennial sea ice (old ice), normal one season sea ice or new formations 
of young sea ice. Based on these estimates of ice thickness and type of 
ice expected and on our previous estimate of ice surface temperature we may, 

with some degree of accuracy, determine the vertical bow reaction force which 

it shall be necessary for our vessel to develop in order to overcome the 
expected ice and force passage so as to accomplish its mission. 

These are but a few of the items of information which past ice history 

can provide to aid the designer. Much more exists, equally pertinent, and 
all of which can be used to aid the designer to more realistically visualize 
his problem and provide a vessel which will be more effective in coping with 
its expected environment. 
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Section D 



Ice Mechanics 
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Part 1 

Ice - General 

Ice may exist in any one of several allotopic forms, depending on the 
ambient temperature and the pressure exerted on the ice. These forms are 

f 

illustrated by the phase diagram, figure / . It can be seen that as tempera- 
ture decreases and the pressure increases from atmospheric up to approximately 
30,000 psi, water solidifies into the common ice form, Form I. As the tempera- 
ture is furthered lowered and the pressure increased to the 30,000 to 31,000 
psi range, Form I ice may be changed to Form II or Form III ice. At 50,000 
psi, Form III ice passes instantaneously through an unstable Form IV (shaded 
area) into Form V at temperature ranges from approximately -13° F to 2° F. 

At temperature ranges from about -13° to -30° F., increases beyond 50,000 psi 
cause direct conversion of Form II into Form V ice. 

Only the common form, Form I, is found under natural environmental 

* 

conditions and has, consequently, practical engineering significance. Figure I 

shows the decrease in the melting point of Form I ice in relation to pressure. 

It is this decrease in melting point temperature that is responsible for 
« 

regulation of ice, i.e., the fusion of two pieces of ice at temperatures near 
the melting point brought together with slight pressure. 

Supercooling 

Even though the normal freezing point of water is 0° C. (32° F.), it 
may be appreciably supercooled before freezing begins, particularly if the 
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water is of high purity. In open water occuring in nature, the degree of 
supercooling varies from a few hundredths of a degree to about - 1.2° C. 

(29.8° F.). After supercooling has occurred, freezing may be initiated 
by seeding or by such measures as impact or splashing. The rate at which 
crystallization occurs depends on the degree of supercooling; in slightly 
supercooled water, such as is normally encountered in nature, the rate of 
crystallization may be sufficiently slow so that supercooled water will 
coexist with ice in quiet ponds or lakes for a measurable period of time. 

Freezing 

The formation of ice, as a rule, depends on the removal of heat from 
the upper surface by convection or radiation to the air. Generally, the 
air temperature must be appreciably below 0° C. for freezing to occur; 
otherwise the rate at which heat is supplied from the earth and from the 
lower layers of water will be greater than the rate of its removal. As 
heat is removed at the upper surfaces during and prior to the formation 
of ice, the colder .water sinks, and warmer, less dense water rises; this 
convection assists in cooling the entire body of water. However, because 
water has its maximum density at about 4° C. (39.2° F.) the convection 
essentially ceases after this temperature is reached by the entire body of 
water. Thereafter, the densest, warmest water remains at the bottom. As a 
consequence, tne rate at which heat will be supplied from the bottom will be 
markedly reduced and the freezing at the surface accelerated. Note here that 
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if the initial ice surface formed is covered with snow, the rate of increase 
in ice thickness will be markedly reduced due to the excellent insulating 
properties of the snow. 

Melting 

% 

Both warm air and solar radiation are important factors in the melting 
of ice occuring naturally. In the Polar regions where air temperatures are 
low, the solar radiation is by far the more important. The radiation is 
absorbed within the ice, particularly by impurities which accumulate at grain 
boundries, hence in fresh water, ice thaws along the crystal boundaries and 
becomes extremely weak while still of appreciable thickness. In sea ice the 
radiation is absorbed, particularly in regions where there is an accumulation 
of salt, causing thawing to start at these sites. Because of the internal 
absorption of radiation by ice, appreciable melting will take place when the 
surface temperature of the ice is still lower than the freezing point. This 
thawing along the grain boundaries and at impurity sites is due to the effect 
of the impurities in lowering the melting point of the ice since all substances 
that dissolve in water depress the freezing point and for dilute solutions 
the depression is proportional to the amount of impurity. 

Effect of Impurities on Ice 

In any body of natural water there will be some dissolved minerals 
and other impurities. Sea water at one extreme is an obvious example but 
even fresh water has impurities. When ice crystals form, they tend to reject 
these impurities which then form a layer of concentrated impurities about the 



crystal at the grain boundaries, and this layer has a lower melting point 

A 

than the crystals themselves. This may greatly effect the physical properties 

of the ice even though the amount of impUrity is small. Melting will ^always 

# 

begin at these boundaries between ice crystals and at a temperature below 
0° F. (32° F.). Fresh water ice becomes rotten as thawing, starting at the 
inner crystal boundaries, separates the ice into needles or columns. This 
ice is said to be candled. As a result of this crystal structure in bulk, 
the physical properties of ice may vary widely, especially when the temperature 
is near the melting point. 

In formation of sea ice, the effect of salt content is very pronounced. 

Freezing of the water does not begin until a temperature of 28.6° F. is reached 

* 

in undiluted sea water, and the structure produced is porous, containing 

P£g 

pockets of brine from which solid salt crystals begin to when 

the ice cools to about 17° F. The structure, and consequently, the strength 
properties' of sea ice improve with time as the concentrated brine drains and 
as the salinity of the ice becomes gradually smaller. Newly formed salt water 
ice of relatively high salinity is flexible and elastic as compared to fresh 
water ice which is harder and brittle. 

The general relation of sea ice strength as a function of temperature 
and salinity is shown in Figure Z . Here we see that as salinity increases 
the temperature must be considerably lowered before the ice has any appreciable 



strength. 
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Crystallography 

It is generally agreed that when ice forms at the surface of a body 
of water it crystallizes in the hexagonal system and that the ice sheets 
which form are usually composed of prismatic ice crystals oriented with their 
long axis (eptic or "c" axis) vertical and at right angles to the surface. 

Hence a normal ice sheet consists of long, vertical ice crystals, each 
usually several centimeters square in cross section. 

The mein difference between fresh water and salt water ice is the 
presence of brine inclusions in the later, forming systematically arranged 
hales. As diagrammed below, in figures -y and o , a sea ice crystal consists 
of parallel pure ice plates in the direction of growth. The brine inclusions 
are between the plates, the amount strongly dependent upon temperature and, 
at a given temperature, is proportional to salinity. An enlarged sketch of 
the brine inclusions shows then drawn as circular cylinders between the laminar 
plates of pure ice. In nature they do not necessarily have this exact shape, 
in fact an el^ptical shape might be more near correct, but this will suffice 
for the purpose of illustration. Note tnat the brine inclusions are shown 
as continuous cylinders although vertical interruptions are usually observed. 

Sea Ice 

Ice met at sea is for the most part of two kinds, icebergs originating 
from glaciers, and sea-ice formed on the sea itself by the freezing of sea 
water. Sea ice proper accounts for probably 99 % of the ice met with at sea; 
bergs are important inasmuch as they constitute a danger to navigation, but 
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they axe common only in a limited number of localities. Suffice it to say 
that from the ship point of view they are to be avoided and no attempt is 
ever made to overcome them. A certain amount of ice may also originate in 
rivers or estuaries as fresh water ice, but it is already in a state of decay 
by the time it reaches the open sea and its importance is no more than local. 

Fresh water freezes at 32° F. but sea water freezes at various lower 
temperatures, depending on the salinity; the greater the salinity, the lower 
the freezing point. Average sea water, with a salinity of about 35 parts 
per thousand does not freeze in the open sea until the temperature has fallen 
to 28.6° F. 

Ice forms first in shallow water, near a coast or over shoals and banks 
and in areas of low salinity. It spreads from these areas as centers. Such 
ice broken up and carried seaward by winds or currents starts further ice 
formation in deeper water. ,Ice in deeper water which has not melted during 
the previous season also acts in the same way. Wave action ordinarily hinders 
the formation of ice to some extent, however, the presence of old ice AP^ps 
down any sea or swell and, at the same time by cooling the water tends to assist 
the beginning of the freezing process. 

The first sign of freezing is an oily or opaque appearance of the water, 

due to the formation of ice spicules and then plates known as frazil crystals. 

% 

These consist of fresh water free of any salt, and increase in number until 
the sea is covered by sludge or slush of a thick soupy consistency. If low 
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temperatures continue, the ice crystals coagulate into small discs on the 
surface. By fusion at their edges these in turn assume the characteristic 

« 

pancake form. As the process continues, these pancakes adhere to each other 
to form continuous sheets of young ice. 

While the original ice spicules and plates are free from salt, a 
network of ice crystals forms as further plates are added, in the pockets of 
which, sea ice is trapped. The crystals grow at the expense of the enclosed 
sea water; the remaining water thus becomes saltier and therefore, heavier. 
Rapidly formed sea ice therefore, holds more brine between its crystals than 
ice of comparatively slow growth. Ice is a poor conductor of heat and the 
rate of its formation drops appreciably after the first 4 to 6 inches have 
formed; a snow cover if present, still further reduces the conductivity. 

The enclosed brine is itself seldom frozen except at very low temperatures; 
it tends to sink through the surrounding crystal network due to its density. 

Sea ice may grow to a thickness of 3 to 4 inches in the first 24 hours 

* ‘ * 

and from 2 to 3 more in the second 24 hours. The original surface of the ice 
is often below water level, hnd the original surface of the ice is rarely if 
ever, visible after the first few days. The ice frequently reaches a thickness 
of 5 feet and sometimes 7 to 10 feet in one season. 

Sea ice, other than fast ice in sheltered bays or along the coast, is 
continually in motion as a result of the effects of wind, tidal streams and 
current. In its motion the ice must open ana shut like a concertina; there 
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are always a certain number of leads and lanes present, otherwise the ice 
could not move. In summer these leads normally remain open, but in winter 
they are soon frozen over with young ice. Apart from the movement of the 
ice, the occurence of swells also tends to break up the ice. Due to these 
various motions, the ice is continuously breaking up, even throughout the 
winter. 

In other instances the ice may not be broken up but rather, subjected 
to pressure as moving floes are driven together or moving floes pressed against 
fast ice. Rafting or hummucking occurs according to the degree of pressure. 
Definite ridges may thus be found, the longer the pressure lasts the greater 
the ridges. Pressure ridges may be as high as 50 feet when grounded against 
a coast, but in deep water away from land the greatest height is from 20 to 
30 feet and more usually from 10 to 15 feet. 

The release of pressure gives rise to lines of weakness in ice fields, 
in the form of cracks, lanes or leads. These are often parallel with pressure 
ridges but an ice field does not necessarily crack in its thinnest part and 
cracks are frequently found passing through ridges and hummucks of considerable 
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The Elastic Properties of Ice 



Ice is generally accepted to be a plastic material and consequently 
should be described by Maxwell's equation: 



in which 
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7" s shearing stress in psi 

tVt 

j- 3 = elastic limit in psi 

y = viscosity in lb/sec/in 

S a rate of deformation 

Jt 

In much of the work that has been reported, it has been assumed that ^ ° 

that is, that ice is truely a viscous material. Although the elastic limit 
is very small, an analysis of the flow of glaciers indicates that it is not 
zero. There is evidence that the viscosity of ice depends considerably on 
the orientation of the optic axis relative to the applied stress. Thus, 
results of tests show ice to be elastic if the stress is applied perpendicular 
to the #ptic axis but otherwise to behave as a plastic. In the case where 
plastic behavior should occur, tests also indicate that if the rate of loeding 
is great enough, i.e. greater than 0.5 Kg/cm^ - sec, the ice will behave as 
an elastic substance rather than plastically. Hence, since the normal mode 
of loadings of ice is parallel to the optic axis, i.e. normal to the ice sheet 
all tests are conducted at rates of loading well in excess of 0.5 Kg/cm^ - sec 



-v>- 



A large number of determinations of the elastic properties of ice have 
been reported, but there is a marked lack of agreement, especially among the 
early results and they must be considered unreliable. Most of the early 
determinations were made using static loading and assuming ice to be truly 
elastic rather than plastic. More recently, measurements have also been made 
using sonic (dynamic) means. In these instances, consistent results have 
been attained and, although possibility of error does still exist, the values 
obtained do represent the most reliable to date. 

Some question does exist as to the assumption that ice is or behaves 
as an elastic isotropic material when subjected to stresses. As mentioned 
above, the elastic behavior may be /Viozcb by the rate of loading, and in 
addition, it has been shown that if crystallization of the ice is sufficiently 
rapid, the ice will consist of many crystals variously oriented and could, for 
practical purposes be considered isotropic. The truth of this last statement 
was borne out by the fact that although specimens for a variety of tests were 
cut in many random directions from ice block samples, no great variation in 
elastic properties were noted. Such minor variations as were observed could 
be explained by internal flaws which were present in some of the specimens. 

The principles involved in the sonic determination of Youngs modulus, 

E, and poissons ratio, V, are based on the principles of the velocity of 
sound waves through a solid medium. It can be shown that the velocity of* a 
sound wave in an infinitely long r&A , called the longitudinal velocity, 
is given by 

r e , 

Ax, = y P 
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also, the velocity of a sound wave through an infinite solid medium, called 
extensional velocity, is given by 
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It can also be shown that ^ the transverse wave over the surface of an 

rum 

extended solid is related to the shear waves which occur in their plates. 
Vibrations of this type were first discussed by Rayleigh in his Scientific 
Papers , Vol. 2, Cambridge University Press, 1900 and are commonly called 
"Rayleigh Waves". A mathematical treatment is given by Love in his 
A Treatise on the Mathematical Theory of Elasticity , 4th Edition, Cambridge 
University Press, 1927. This velocity may be written as 

i/r = £ y . 3 ^ ( / / J.'-J 

Therefore, if we can determine these three velocities for ice experimentally, 
we can determine Youngs Modulus and poisson ratio for the material since 
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once T"" determined from the above £ may be obtained from 



- yy- 



and checked by determining K from 
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Tests have indicated that the values of Youngs Modulus varies somewhat with 
temperature, especially in the range - 30° C. to 0° C. The results of these 
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experiments are shown in Figure 3 . Values of Poissons ratio remain relatively 

constant over this temperature range and may be taken as 

A"'- r.JAs' + o.oo P 

from 0 t©-15° Centigrade. 

These values in Figure 3 are representative of the variation of Youngs 
Modulus for fresh water ice. For the case of salt water ice there is evidence 
that the values of Youngs Modulus should be considerably lower, more a function 
of temperature and definitely dependent upon brine content (temperature and 
salinity). Unfortunately, however, little if any, experimental data exists 
in the present day literature to corroborate this evidence and show just what 
the range of values of Youngs Modulus for sea ice should be. 
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In addition, since the paper deals with the icebreaker problem, some 
question must be raised as to the accuracy of the sonic methods for the 
determination of Youngs Modulus. The sonic stresses are very small while the 
working stresses due to the operation of an icebreaker are very large and it 
is possible that the sonic method gives somewhat inaccurate information of 
stress-strain relationships for use in solving the icebreaker problem. 

We must, however, accept these probable inherent inaccuracies since 
the data presented above represents the most accurate and reliable published . 
to date. Possibly, future experimentation will provide a more accurate and 
realistic representation of the values of Youngs Modulus for sea ice. For 
the icebreaker problem, however, it is felt that the values quoted above are 
adequate. Since the ice sheet is considered as a semi infinite flat plate 
on an elastic foundation, and since only the fourth root of Youngs Modulus 
enters into the determination of stresses within the plate, it is felt that 
what errors exist in our above listed values of Youngs Modulus could not be 



such as to materially effect the determination of stresses within the ice sheet 
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Part 5 

I 

Ice Friction 

The experimental relations known as the classical laws of dry friction 
have been known- for several hundred years. These laws as stated by Coulomb 
are: 

1. Frictional force for a constant load is independent of apparent 
area of contact. 

2. Frictional force is directly proportional to the load, that is, 
to the total force which acts normal to the sliding surface. 

E/iufric 

3. Static friction force differs from hea*«44c friction force; 

however, the friction is independent of the velocity 

of the slidea. 

4. Frictional force depends upon the nature of the materials in 
contact. 

The second of these laws of dry friction may be expressed as 

= t/ AJ ' 

/ 

where 

Pr s frictional force 
AJ ' - normal load 

- ■ A7 KlfT/t 

s coefficient of friction - either static or hs r et-l c 
These static or coefficients, according to law four, are dependent 

upon the nature of the materials in contact, but, according to law one, 
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independent of the apparent area of contact. 

The characteristics of hydrodynamic friction, on the other hand, are 
quite different. In many cases it is exactly opposite to the characteristics 
of dry friction. For a Newtonian fluid the force of resistance is related 
to the apparent area of contact, the velocity gradient, and to the viscosity 
according to the relation 

F = ^ J L 

A x 

where 

JF a Tangential resistance force 
on = Absolute viscosity 
A ^ r Apparent contact area 
J- r /j i z velocity gradient 

Measurement of the frictional force, its variation with load, 
apparent area and velocity enable one to distinguish between the presence 
of dry and fluid friction. These criteria have been applied to friction 
on ice. Experimental measurements do indicate that the frictional force 
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is nearly a linear of the load, and related to the materials of 

the two surfaces. The variation of frictional force with apparent area and 
velocity is very small. These several relations indicate that friction of 
materials on ice is primarily of the dry friction variety. 



The resistance experienced in the process of sliding one surface over 
another, which is generally termed friction, is not a simple phenomenum but 
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probably arises through the action of several different mechanisms. 

One mechanism of dry friction is undoubtedly related to the presence 
of asperities on the surface. Several theories have been proposed to explain 
this relationship. According to this surface roughness theory, the asperities 
are considered to be rigid with an average slope angle equal to ft . The 
applied tangential force supposedly has a component parallel to the slope 
which is sufficient to lift the normal load over the rigid asperities by the 
relation 

1 1 -9-0 / — AJ ^ 

hence coefficient of friction 

c/ - — r / y/ 

AJ ' r 
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In 1920, however, Hardy pointed out that the lifting mechanism cannot 
be the principal mechanism of frictional resistance since experiments with 
flat clean metal surfaces indicated friction practically independent of surface 
roughness. He proposed adhesion as the dominant force of friction. Holm in 
1938 suggested that this adhesive force acts only over the real area of contact 
which is a fraction of the apparent or geometric area of surface contact. 

As can be shown, this real area of contact between a sliding body and a 
stationary surface composed of a single material is a function of the normal 
load applied and may be expressed as 




where s plastic mean flow pressure of softest material. 
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The resistance force due to friction therefore equals the product of the 
shear stress in the contact region times the real area of contact, i.e. 

- 7j 



where 



Therefore 



7^ * shear stress of weaker material 
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f ^ r 

This indicates a linear relationship between friction force and normal load 
which satisfies the second classic law of friction. Since we said coefficient 
of friction 
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where both shear strength and mean flow pressure are bulk properties 

of the weaker material, and explain classic law four regarding the dependence 
of friction force upon the nature of the materials in contact. 

Two factors which effect the strength of a material and hence 
and are temperature and time. It was reasoned, therefore, that the 

temperature variation of the coefficient of friction should be a critical 
test between surface roughness and adhesion theories. Since the geometric 
slope of the asperities is independent of temperature there should be no 
change in with temperature according to the surface roughness theory 
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but since both "Tj and ^ are temperature sensitive, a temperature variation 

of M would be expected by the adhesion theory. Experimental measurements 

* 

indicate that the temperature variation of was very similar to the ratio 

of the variations of Tj' and , hence , tending to substantiate the 

adhesion theory. 

If we apply this general theory of friction to the case of friction 
on ice we have 
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The mean flow pressure f ^ is often taken as 1.1 T where Y - yield stress 
or proportional limit of the material. In the case of ice, the elastic limit 
has been reported as lying between 7 and 14 psi. These values are low and in 
conflict with the assumption of ice as an elastic medium with elastic behavior 
up to the yield point and with a tensile strength of about 200 psi. Some 
justification lies in the belief that these reported values of 7-14 psi are 
but the limits of proportionality and that the elastic limit must be higher. 
However, as an approximation, if we assume 

^ = J O - P - i 
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Then r 
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This calculated value is much higher than experimentally obtained values. 

1 

The* lower actual values of friction on ice could indicate error in the measure- 
ment of ^ and or, as is more reasonable, indicate the presence of an 






intermediate lubricating fluid. The mechanism by which this lubricating 
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fluid - water - was formed or converted from ice was not, however, easily 
explained. 

( 

In the 19th century Thompson presented papers regarding the possibility 

of ice melting due to the lowering of the melting temperature by the applica- 

( 96 ) 

tion of pressure. Reynolds, in 1901, supported this theory by noting that 
several investigators had experimentally found a deorease in the slipperiness 
of ice at low temperatures. Reynolds interpreted this fact as indicating 
insufficient pressure available to lower the melting point to the low ambient 
temperature. Quantitatively, however, it must be noted that approximately 
100 atmospheres of pressure are required to lower the melting temperature 
by 1* C. The tremendous pressure necessary to supply this melting temperature 
depression places the pressure melting mechanism in doubt as the source of 
lubricant supply. As regards Reynolds observation of temperature dependency, 

in the case of sea ice, it can be attributed to the precipitation of solid 

(#>) 

salts at the lower temperatures as explained by in 1958. 

(•n>) 

Bowden and Hughes agreed that water must be the lubricant but felt 
that its formation was the result of frictional heating of the ice-slider 
interface. Depending upon the relative thermal conductivity of the slider 
and ice, part of the frictional heat is lost to the slider and part is used 
in melting ice. Thus, theoretically, the coefficient of friction between 
brvisi and ice should be greater than between ebonite and ice. Experiments 
substantiated this fact. In 1947, however, this theory was cast in doubt. It 
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was reasoned that if this heat melting theory were true, direct heating of 
the sliding runners would facilitate melting and further reduce friction. 
Experimental results did not show the expected frictional reduction. 

(9 b) 

Bell in 1948 also questioned this friction melting theory and proposed 
instead a plastic flow theory in which pressure and frictional heat caused 

increased plasticity of the solid ice surface, rather than a complete phase 

(/z(J 

change to liquid. However, work of Nakaga and Matsumato in 1953 shows 
experimental evidence of the existence of a liquid water film on ice surfaces 
below temperature of 0° C. This liquid film is supposedly present at all times 
and is not supercooled water, but rather, water in equilibrium with its vapor 
phase on one side and with ice crystals on the other. What effect this film 
has on frictional resistance or whether it is the one theorized to have been 
due to pressure or frictional heat is not known, however, despite Bell’s 
contention of plastic flow, and the inability to conclusively demonstrate 
experimentally the pressure or friction heating theories, it does seem 
pretty conclusive that what lubrication does occur during sliding on ice is 
due to an interposed liquid film. 

Eespite all the theories, and published experimental results, very 
little actual reliable data exists on the static and dynamic coefficients 
of friction of steel on ice. The variation of these coefficients with 
temperature, though recognized, has not been determined with any degree of 
reliability. The effect of normal loading is also elusive. Previously we 
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stated that the real area of contact could be represented as 

which indicates that it is independent of apparent area of contact as 
measured by standard geometric means and directly proportional to the 
first power of the normal applied load. This is in agreement with the 
first two classic laws of friction but is only applicable for a low ratio 
of loading to mean flow pressure. At higher ratios of loading to mean 
flow pressure, the real area of contact is an appreciable portion of the 
apparent area. Since the real area cannot exceed the apparent area, the 
ratio A>- must approach a limit of unity at extreme loads and may 

be expressed analytically as ; 

f) r 

ff. - • - * 

where B is a constant 



hence -*.i 

-* <3 /V r » O 

Al' » A r ~ 

In addition, for lew values of normal load the expression 

- 4 \! ' 

C X, / - 



Therefore 



A r 
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and real area is a linear function of normal load as required by the second 
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classic law of friction. 

3 = 



In this instance of low normal load 

^r_ / 

ck_ X so ' 



but 

AJ 

£) r ~ jo 
* /+1 

Therefore 

6 ’ 4^ 

Now this real area - load relation is recognized but has not, unfortunately, 
been confirmed experimentally over* the entire loading range. Reliable deter- 
mination of at low normal loads and of the constant d at higher 

loadings are lacking. Hence, though we recognize the effect of normal loadings 
as it effects real contact area and hence friction coefficients, our expressions 
still require reliable experimental verification. 

In addition, it is known that a snow cover will materially effect 
the coefficients of dynamic and static friction. However, since it has been 
shown that the effect is dependent upon the type of snow cover, water content, 
temperature, normal load and velocity, what data does exist on the effect of 
snow on friction, is unreliable and difficult to use. 

Consequently, for the ship problem we must attempt to select those 
values of coefficients of friction which appear to be most reliable and con- 
sistent with published experimental results. It must be kept in mind, however, 
that the values are only representative and not necessarily applicable to all 
engineering situations. For the case of the ship problem, the most represented 

( 96 ) 

tive values are those published by Arnold and Alabrieff (1938) which are given 



for polar salt water ice as 



coefficient 

coefficient 



of static friction 0.30 to 0.50 

of dynamic friction 0.10 to 0.20 
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Part 4 

The Tensile Strength, of Sea Ice 

The most comprehensive study of the strength of see ice to date is 
that made by Dr. ft. AaoH of the U. S. Army Snow, Ice and Permafrost Research 

Establishment in his paper Composition of Sea Ice and its Tensile Strength , 

(So) 

(1958) . In this paper, which is partially paraphrased here, Dr. A ^ >->£- attempts 
to analyze the strength of sea ice in terms of brine content rather than in 
terms of salinity and temperature as has often been attempted in the past. 

The laminar structure of sea ice crystals, shown in figures 4 and -f , 
where plates of fresh water ice are separated by layers with brine inclusions 
has been known since the turn of the century. The earliest attempt to study 
the relation of strength to internal cavities in sea ice is that of Tsunkov 
(1939) whose relation 

<5 / = *3. t ' - * £ ' ^ ' 

where a strength of ice without pores 

m porosity - internal cavities without ice - 
Tsunkov used air content 

is based on the assumption of cavities in the form of uninterrupted circular 
cylinders. 

Further progress was made by Weeks and Anderson in their study of the 
geometry of brine inclusions, leading to the formulation of sea ice strength 
as 

^ ' = y~ c / - a ; 
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where ~z>^ = strength of fresh water ice 

X e stress concentration factor for brine inclusions 
and proposed two models for 8. i.e. 

x = j 

rr 

^ D • Z//6 ~lT /- ra . &&// 

cd, o jed 




where - width of brine pockets * 

b 0 m spacing of brine pockets 

The first of these models is based on the assumption of circular cylinders 
as was Tsunkov’^. The second is elliptical and assumes changes in the brine 
pocket shape within a plane only. 

THU 

The formulation for sea ice strength was a novel approach to sea ice 

physics in that it is similar to the approach used to compute strength of 

perforated plates and it does try to explain the known fact that the flexural 

strength of sea ice measured while in contact with the water is far below 

* 

that as measured in air in the laboratory. However, this difference in 

tensile strength between fresh and salt water ice explained in terms of stress 

concentration factors is not supported by experimental evidence. In addition, 

neither this nor any previous formulation attempted to explain the more com- 
1 \ 

plicated case of reinforcement of the ice due to percipitated salts at low 

temperatures, i.e. below about - 8*C 

The present development by Dr. Ass^e, therefore, is based on the con- 

cept of 

'= y, ( / - 



where 






; 










1 



i 

I 
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'8 - strength of sea ice 

Z' s friction representing the reduction in strength due 
to brine inclusions 

•v s c 2s , ^ -- 3 basic strength of 

sea ice computed directly from 
sea ice data 

Z strength of fresh water ice 
C s coefficient 

~$p z strength of fresh water ice at given temperature </> 

4 

"2 ^ ) a strength of fresh water ice at a temperature of -10° C. 

Note the factor — is needed in order to reduce the strength of 

fresh water ice to ->10° C. so that the effect of temperature upon the strength 
of fresh water ice bridges (i.e. platlets) in a sea ice crystal need not be 
considered, and what remains is solely the effect of brine inclusions. The 
choice of -10° C. here is an arbitrary one. It is interesting to note here 
that for most instances the value of C < / . This can be explained by the 
fact that all tests for the strength of fresh water ice are conducted in air, 
on core samples, by the ring test, i.e. loading a hollow cylindrical core 
perpendicular to its axis until failure. Now we know that ice in the presence 
of its liquid phase is weaker than ice tested in air, therefore C should be 
less than one since the sea ice has brine inclusions, hence the interstitial 
liquid weakens the minute ice bridges. 

In order to determine Z' = j. { r'j a general three dimensional 
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geometrical model is developed to account for the effect of brine volume on 
the ice strength. In order to do this, of course, the brine volume must be 
determined as a function of temperature and salinity considering the percipi- 
tation of salts and numerically deriving the phase relations in sea ice. It 
is interesting to note the radical effect on sea ice produced by percipitated 
salts such as J • /o J and c c ■ z a! l o 

It is because of these salts that we can no longer maintain, as many publica- 
tions still do, that sea ice is three times weaker than fresh water ice. 

It is true that sea ice can be weaker, but at low temperatures, i.e. below 
about -8° C., where these salts begin to precipitate the sea ice can e%ual 
or even exceed the strength of the fresh water ice. The present theory, 
therefore, attempts to explain the wide variations in sea ice strength in 
terms of brine volume and precipitated salts. 

Based on many actual tests it was found that at a given temperature 
the reduction in strength of sea ice is proportional to y salinity. But 
since brine volume at a given temperature is proportional to salinity then 
the reduction in strength is proportional to f v' . Therefore, in the 
selection of a geometrical model for ice crystals with brine inclusions it 
seemed desirable to have one which represented the reduction in ice strength 
due to brine inclusions as a function of i 

If we introduce the notation as shown on figures o' f 6 
-\r ' w relative volume of brine 

a" 3 average area of a brine pocket in the & c 7 plane 
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<2-0 5 average distance between plates 
b „ = average spacing of brine pockets 

« length of brine cylinders in direction 
s average spacing of brine cylinders 
Then we may say that 
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( 1 ) 



Introducing the dimensionless parameters 

x -- i* 



equation (1) becomes 



tr 



A 



(2) 



The most convenient way to measure ^ is to find the maximum length 2 r 6 
and the maximum width ^ C of individual brine pockets in the ^ plane 
and averaging 












under the assumption of an elliptical shape which is a good approximation. 

If a tensile stress is applied in the direction of the Q axis, then it acts 
upon a reduced cross section of 
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If the basic strength of ice itself is designated as 25^ then the 
effective strength due to the reduction in cross section becomes 



^ 0 - 7 . 






(4) 



Now all depends on how z and ^ vary with the brine volume -v ' 

If we assume that the average length and spacing of the brine cylinders 
and therefore i stay constant then the changes in brine volume are reflected 
only in the cross section in the d -<r r plane. If we assume also that this cross 
sectional shape of the brine pocket stays geometrically similar then all linear 
dimensions of the brine pocket in the cross section change proportionately to 
the square root of the brine volume, therefore 



2 ^ ■=■ 1 ~ i / v ' y 

Based on many tests with actual sea ice this assumption appears most reasonable 
and does produce a model which satisfies our previous criteria that the ice 
strength decreases as a function of Tyf’ 

If we introduce another non dimensional parameter 



A 



then 



i!L 

v? 




and from (3) 

^ ^ Ci 

A' 0l o 




also from (1) we see that 
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so that 
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Therefore the non-dimen sional criteria 



(5) 



t o y 

i — / 

has to be constant if the strength of ice is to reduce proportional to y-v-' 
how if we consider the brine pockets as elliptical cylinders and let 
be the average ratio of the length of the elliptical brine 
pockets (in the d direction) to the width of the pockets (in the <r direc- 
tion), then 

- rrr • - 77" as before. 
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Since 
we have 
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Therefore, equation (5) becomes 
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From microstudies of ice crystals, reasonable values of /3* 
and 6 are determined so that the petrographic constant ^ >_ may be determined 

'T = ^ » 1.5958 

The ro fore 

^ — ) (?) 

r y t (/- /• 

In order to determine it was necessary to evaluate the results 

of a considerable number of ring tests performed on sea ice core samples in 
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the field so as to obtain the trend in variation of tensile strength as a 
function of brine volume. This, in essence, proved to be a straight line ***<*> 
for the region above -8.2° C y (prior to salt precipitation) veriattna, 
and since we've shown that 

% 

2 ' - X - tC ^ /V' 

we can equate the Blope of this straight line variation to 2 C / t and thus 
it was determined that 

- ^<l■6 6 t O.yo 

j t cr>. j / *4/ c „ <- 

hence 

~6 - //as (/-/.O /7j ( 8 ) 

for that portion of the sea ice condition above -8.2° C. 

The above test results also indicated that as the ice temperature 
continued to rise a point was reached where the ice strength fell off very 
rapidly (deterioration - complete loss of strength). For this region, the 
test results indicated the e.^ fi ir^c^L relationship 

t r -1 (9) 

*, =• H> l ' - Vo ( -U-' - o. 

As previously mentioned, as the brine volume in sea ice drops - tempera- 
ture decreases - various salts begin to precipitate out which act to strengthen 
the ice (i.e. below -8.2° C.). Figure represents a model of reinforced 

brine pockets. For simplicity, circular pockets are drawn. The initial radius 
>7 shows the brine pocket when the salt started to precipitate, r" is the 
radius after salt is precipitated, but the pockets are not drawn to scale. 
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If the salt particles together with ice are deposited on the walls of the 
brine inclusions, there is effective reinforcement at the place of highest 
stress concentration (in the B-B plane at the walls). The small bubble in 
the center is a symbol for stress concentration or local defects on the B-B 
plane. 



If Youngs modulus of this salt ice mixture does not differ much from 
its value for ice, all that matters is the reinforcement, which will act as 
soon as a small amount of salt is precipitated. Since a rapid transition 
from a high stress concentration to normal stress occurs near the wall, the 

f 

width of the salt-ice reinforcement does not matter, once a certain minimum 

I 

value is exceeeded. This model holds well for the precipitation of ■/& » x o $ 

•u l.< i » L c as evidenced by considerable test results. Using a test evaluation 
procedure similar to that used to determine ^ we may evaluate ^ as 

7 c /& ? - =>■ <7/ *~3/c '■ 



Therefore 
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for the region between -8.2° C. to -22.9 C. where oO y • / j h u 

salt precipitates out. The actual increase in strength, therefore, due to the 



presence of the salt is 

s&. 
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or an increase of one third. 

Now let us consider the strength behavior of ice when a LC l precipitates, 
which it does in much larger quantities than . Test results indicate 

that the presence of NaCl (in the form of NaCl' increases the ice strength 
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appreciably though, it was at first thought, in another manner than does 

. Later indications seem to indicate however, that the model 
used for precipitation of J y •'' <,// t oalso applies to the precipitation of 

-zh^o and that the value of could be evaluated as for the 

case of ^<> v -> 0 y • /o/v L (personal correspondence). This produces a value of 



<i J. 6-i' 



V/c 



therefore 
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for the region below -22. 9° C. where <i ■ z H u o salt precipitates out. 

The actual increase in strength, therefore, due to the presence of this salt is 

^■bj' 

/«■ • 4 o ' 

or an increase of 2/3, 

Having determined the relationship of tensile strength of sea ice as 
a function of included brine volume, it was necessary to determine the varia- 
tion of brine volume with temperature so as to be able to show variations in 
tensile strength of sea ice with temperature. 

The only attempt to compute the brine content of sea ice was that 

( S'”) 



reported by Malingren (1933) and £uaos (1945). 
in common, that is, that brine content K is 









where 



These attempts had one thing 



a salinity of the ice 
e salinity of brine 

This equation states that all the salts measured in s are dissolved in the 

brine inclusions, which is not so, and the correct relationship should be 

- 
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where S 3 a relative amount of solid salts. 

Tiie assumption that ->j - is not reasonable nor permi3sable. 

In order to determine the relative brine volume in sea ice, therefore, 
it was necessary to determine the exact concentration of different ions and 
salts in sea ice, after the freezing of sea wator, as a function of temperature 
at the salinity of cornel sea water so as to determine the phase diagram for 
sea ice. For given sea ice the values shown in the phase diagram for salts 
and brine have to be altered in proportion to the salinities, however, from 
this information it is possible to determine the amounts of precipitat d salts 
as a function of temperature. Where the amount of precipitated salts were 
computed, the amounts of ions in solution could be determined by subtration 
from the initial amounts known to exist in sea water. Once the absolute 
amount of ions were known, their total was obtained and the relation 







could be applied to obtain the amount of brine. Dividing the amount of brine 
of 34.325 (the salinity of normal sea water in parts per thousand) gives the 
amount of brine in standard sea ice of one part per thousand salinity as a 
function of temperature. Since theoretical densities of sea ice vary little 
depending on temperature and salinity, an average value of 0.926 was assumed 
so that the relative brine volume could be computed as strictly proportional 
to salinity. In this way the variation of relative brine volume could be 
computed as a function of temperature for a salinity of 1 part per thousand 
and multiplied by the salinity to get variation with temperature at any other 
salinity . 
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Knowing this variation of relative brine volume as a function of 
temperature for a given salinity, it is possible to apply equations (8), (9), 
(10) and (11) in order to obtain the variation in tensile strength as a 
function of temperature. These results are shown in figure 7 for representative 
salinities met in sea ice. The strength is given in terms of a computed basic 
strength of 14.20 Kg/cm^ (202 p>si) equivalent to the tensile strength of an 
imaginary sample if all the brine pockets are filled with ice. It should be 
recalled that values shown in figure 7 represent the tensile strength a3 a 
function of ice temperature ana not ambient air temperature. 

For the portion of figure 7 below -8.2° C. and -22.9°C. the sudden 
increases in strength is due to the presence of precipitated Nag 
below -8.2* C. and precipitated A/ " 1 Cl ' ° below -22.9° C. The presence 

of these salts, even in small amounts, increases the strength of the ice by one 
third in the case of the ^ >. and by two thirds in the case of the 

^ c u • i //,. o . There is very little change in strength with temperatures 

in this region. Note that although the strength increase is shown to occur 
instantaneously at -8.2° C. and -22.9® C., there is in reality a finite transi- 
tion temperature zone. However, since the mechanics of this transition is 
complicated and the transition temperature zone small, we may, for practical 
engineering application, ignore it and consider the transition instantaneous. 

In the case of perennial sea ice, it is suspected that a great part 

fitiovC a 

of the <5y . /j wq o does not redissolve -sbeai -5.2 C. 

This explains the fact that the strength of perennial ice is about one third 
above the strength of comparable one season ice even if the effect of salinity 
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is considered. Hence the higher values for perennial sea ice in the tempera- 
ture range of -1.3° C. to -8.2° C. Above -1.3° C. all the J 5 / ■ / & & >. o 

goes into solution and the ice strength breaks down completely. 

As previously mentioned, the temperature involved in the sea ice strength 
variation depicted in figure , is the ice temperature and not air temperature. 

Since we wish to use this data to predict icebreaker performance in the preliminary 
design stage it would be advantageous to be able to convert seasonal air tempera- 
ture to ice temperatures so as to facilitate use of figure 7 . In addition, 
since for our icebreaker problem we will consider the ice sheet as a semi- 
infinite flat plate on an elastic foundation which is caused to fail due to 
bending stresses in the outer fiber, and since we will subsequently show that 
it is the strength of the top most fiber of the ice sheet which dictates its 

<r 

breaking we can see that we are really interested in the ice surface temperature 
for application to figure 7 * Doctor in a subsequent work (Unpublished - 

personal correspondence) has developed a relationship between ice surface temp- 
erature and air temperature which is ideal for our problem here. This relation- 
ship is given for conditions of negative air temperature as 

-ac V (12) 

where 

-O; r temperature at the interface between ice and snow 

s air temperature (°C.) averaged over a number of days 
depending on For iuf 3^ - 60" it is 

4 days* 



-?s - 



s water temperature (°C.) which is dependent upon salinity but 
can usually be taken as -0.17°C. 

- constant dependent upon location of ice. For measured 
at a shore station = 0.84. For ^ measured above the ice 

sheet as for an icebreaker i.' ^ ©.<9o 
n - i° e thickness - cm 
i m snow depth - cm 

= relative insulation of snow versus ice. •?<. is 5 to 6 for normal 
conditions in the arctic, 4 to 5 under extreme cold for wind- 
swept snow, 8 to 12 for old snow and 12-19 for a normal snow 
cover during the winter in moderate latitudes. 
t> » a temperature value which is about f 1.0°C. in the middle of 
the winter, increasing during cooling periods and decreasing 
during warming periods. During late winter b’ decreases to 
about -0.4°C. 



For the purpose of the paper, however, we shall not consider the effect of a 

F Slc-ri oki 

snow cover since so little is known of its effect on ice f-use^ien, its com- 
pressive strength, etc. Consequently, we shall consider ice as without a 
snow cover, hence h s = 0 and our expression above becomes 

v. ^ o °C (13) 

Now assuming the values suggested by Dr. A a ^4 for and b and using 
^ we have 
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How it is true that there is some question as to the choice of b. However, 
since the modern arctic icebreaker attempts to get through the ice as soon 
as is possible at the close of winter it seems logical to use that value 
for b which is characteristic of late winter. It should be reiterated that 
the above expression does not hold for positive air temperatures and breaks 
down, in general, in late spring. It is desired, however, that our icebreaker 
be capable of penetrating the ice as early in the season as possible, preferably 
no later than early spring? Consequently, temperature conditions in the 
Arctic and Antarctic during this time will be such as to make expression (14) 
applicable. 

As a guide to the mean air temperature conditions existing in the 
Arctic end Antarctic during various months of the year, table l.aad 2 show 
representative mean temperatures experienced at various stations throughout 



the year. 
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Weather Data for the Antarctic 
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Argentine 

Islands 


Cap© Adare 


Cape 

Denison 


Debenham 
I elands 








LAT oo-loS 
LONG 64-16W 


71-18S 

170-09E 


67-OOS 

142-40E 


68-08S 

67-06W 


rroxanemt 

78-383 

163-37W 


WAN 


MAX 


39 


1 37 


34 


39 


27 




MIN 


31 


30 


26 


30 




|FLfl 


MAX 


- 


- 


29 


56 


1 




MIN 


- 


- 


21 


27 




MAR 


MAX 


36 


22 


18 


34 






MIN 


28 


14 


10 


25 




APR 


MAX 


30 


16 


6 


> 

30 


12 




MIN 


23 


5 


-2 


1 20 




MAT 


MAX 


28 


3 


6 


14 


-04 
- 4 




MIN 


21 


-11 


-2 


1 


C/\ 


JUNE 


MAX 


27 


- 5 


3 


10 


“O? 

13 




MIN 


21 


-19 


-3 


-10 


-7^ 


JULY 


MAX 


18 


. - 1 


0 


» 


fw 

9 




MIN 


7 


-18 


-8 


-6 | 


-66 


AUG 


MAX 


5 


- 4 


6 


17 


-12 




MTN 


-7 


-25 


- 3 


It 

-5 | 


O A 


SEPT 


MAX 


18 


- 4 


4 


25 


-74 

15 




MIN 


5 


' -19 


-6 


13 




OCT 


MAX 


30 

f 


5 < 


8 


± 2 . 


— 04 
16 




MIN 


22 


-11 


-3 


% 

19 


_ah 


NOV 


MAX 


32 


24 


22 


36 


TV 

23 




MIN 


23 


9 


21 


26 


.IQ 


DEC 


MAX 


36 


36 


31 


40 


•lo 

32 




MIN 


30 


27 


23 [_ 


29 


2 



TABLE 1 (Cent. 
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Cr * ° 


Laurie 

Island 


Little 

America 


Port 

Charcot 


Port 

Circoncision 


Water 

Boat 

Point 


66-026 


60-453 


78-343 


65-033 


j 

. 65-103 

i 


64-483 


89-54E 


44-43W 


163-56W 


64-OOW 


64-12W 


62-43W 


34 


42 


31 


35 


i 

38 


50 


24 


24 


18 


30 


33 


27 


30 


43 


11 


36 


38 


49 


22 


24 


-19 


27 


31 


25 


23 


43 


21 


37 


38 


42 


10 


18 


-45 


23 


30 


20 


10 


40 


10 


36 


27 


47 


-3 


6 


-49 


10 


19 


20 


11 

i 
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Part 1 

Icebreaker Propulsion Machinery 

The selection of main propulsion machinery for icebreaking vessels is 
completely governed by the severe requirements of ice breaking service and 
the long range operation it requires. The need is for a propulsion plant 
which is extremely reliable, rugged and capable of withstanding the frequent - 
and often extended periods during which the propeller and shafting are sub- 
jected to serious and sudden shock from both solid and broken ice. 

When breaking a channel through an ice field or forcing a lead, the 
propellers often operate at bollard condition while at the same time often 
sustaining hard serious contact with pieces of solid ice forced into the pro- 
peller stream by the ice breaking action of the bow. While operating against 
heavy pack ice, the engines are continuously reversed rapidly in order to 
allow the continuous maneuvering necessary to keep the vessel from becoming 
beset and to augment the rudder in forcing a straight channel through the ice. 

The plant must be capable of developing maximum thrust {hence maximum power) 
at zero speed of advance, requiring the development of practically full power 
at reduced shaft HPM. In addition, due to the nature of icebreaking, the propul- 
sive system must be amenable to direct, positive and simple control from the 

bridge. 

c 

Few kinds of propulsion machinery offer the degree of flexibility, 
ruggedness and reliability required for this icebreaking service. Recipro- 
cating steam engines have been successfully used in the past while Diesel 
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Electric systems are most often the choice for present day applications. 

Both types of plants offer the necessary characteristics of high continuous 
torque over a wide range of shaft HPM. Each system has its inherent advan- 
tages, and -although steam is seldom used on modern day arctic icebreakers, 
it must be stated that this type of machinery has sn excellent record of ice 
breaking performance in the many vessels in which it has been used. Among 
the advantages enjoyed by a reciprocating steam plant when in comparison with 
a Diesel Electric unit are: lower first cost, lower maintenance 'cost, a 

requirement for less qualified personnel and the fact that the plant is 
reliable and simple to operate. Among the disadvantages may be listed higher 
total plant weight, higher weight of fuel and water, greater plant space 
required, less flexibility in a multiscrew multiple engine installation and a 
larger marginfe of excess fuel capacity required for any given cruising range 
during which icebreaking work is anticipated due' to the greatly increased 
specific fuel consumption when breaking ice. It should be noted that for a 
comparable normal 15000/16500 SHP installation, the estimated all purpose fuel 
consumption for the two type plants is: 

Diesel Electric - 0.50 lb. at full power 

Steam - 0.86 lb. at full power for a 440 psi plant 

Admittedly, there are many points in favor of the_ steam plant installation. 
The present trend, however, is towards a greater general adaption of Diesel 
Electric propulsion. The most important reason for this trend in icebreaking 
vessels is the lower fuel consumption and correspondingly larger cruising 
radius. For arotio breakers, the multiplicity of power generation units is an 
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advantage since it is possible to cruise at lower speeds with the engines 
operating at maximum efficiency. Future developments may, of course, 

reverse this present trend. The adaption- of nuclear power to icebreaker 
power plants may well be the greatest stimulus yet foreseen. Fuel consumption 
would no longer be a controlling factor since tne use of nuclear power would 
allow continuous icebreaking operation for prolonged periods of time. Further- 
more, the fears of being beset by pack ice for the winter period will no 
longer plague operating personnel, since there would be no problem of fuel 
shortage to force drastic reductions in heating and other auxiliary functions, 
and the vessel could become fully operational as soon as ice conditions per- 
mitted. It must be noted, however, that even if nuclear power made a return 
to steam desirable, the propulsion system would in all probability be of a 
turboelectric d-c direct drive type. Hence, in reality, the present trend 
would not really be reversed, but rather, result in a combination. of the two 
systems. However, since only one nuclear icebreaker actually exists (the 
Russian "Lenin*) and since no others are presently contemplated, let us, for 

^ 4 

the purpose of this paper consider only the Diesel Electric installation. 

The Diesel Electric system is essentially a constant power system, though 
the torque and shaft KPM are still in direct proportion to the horsepower. 

The design speed for this type unit is that speed at which the propeller has 
been designed to absorb full horsepower at the highest effeciency. If the 
brake horsepower of a diesel unit is known, the shaft horsepower may be 
estimated witn an assumed efficiency of 94# for the generator and 92# for the 
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propulsion motor - hence, an overall efficiency of 86.4%. 

For those vessels with bow propellers it is possible by the sacrifice of 
a relatively small regulating effect to distribute the available total power 
in an arbitrary proportion between the forward and after screws. Normally, 
the power is equally distributed amongst all screws, however, when breaking 
ice, all power may be applied to the after screws. Several existing installa- 
tions are so flexible that as much as 75% of installed power may be applied 
to the forward shafts, which is often advantageous when operating in rafted 
lake or river ice. 

In the arctic icebreaker, bow propellers are seldom installed since 
service experience has shown that ttere presence is often a serious disadvan- 
tage. Consequently, all power is normally distributed between two after pro- 
pellers. If a triple screw arrangement is used, half total power is normally 
applied to the centerline screw and one quarter total power to each of the two 



outboard screws 
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Part 2 

% 

Propellers and Thrust Predictions 

Propellers for use on icebreaking vessels are normally designed so as 
to develop maximum thrust at a zero speed of advance. Some authorities ad- 
vocate a somewhat higher design spaed, in the vicinity of three to six knots, 
since it would be desirable to develop maximum thrust at the speed at which 
an icebreaker normally is able to proceed through one year pack ice. However, 
as shown by Bridewiser in a thesis entitled The Estimation of Tugboat Tow- 

(9 8 ) 

line Pull and Effects of Types of Propelling Systems , the effect of design 
speed on developed propeller thrust is very small for the diesel electric 
installation. For example, at a ship speed of three knots the percentage 
difference between thrust developed for a propeller design speed of 0 knots 
and 6 knots is only 2.2$ and at a ship speed of six knots this percentage 
difference is only 1$. This very small percentage difference is mainly due 
to the constant efficiency of the power plant and what little variation that 

does exist is due to the slight variance in propeller efficiency. Consequently, 

ot-> 

the propeller design speed is not too vital a criteria when considering the 
thrust developed at speeds of advance through ice pack. When breaking heavier 
ice by ramming, however, it is quite easy for the vessel to become stuck due 
to riding up too far onto the ice sheet. In such an instance it is highly 
desirable that the propeller develop the maximum possible thrust so that the vessel 
may extract itself. Heeling tanks are often a great help at such a time in order 
to break the vessel free from the ice. However, the thrust is still necessary 
in order to insure the commencement of astern motion. Since a difference in 
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propeller design speed of from zero to six knots c&n modify the bollard pull 
by as much as 11$ it would seem logical that a propeller design speed of zero 
knots is the most desirable for the arctic icebreaking vessel. 

The propeller diameter should naturally be chosen as large as possible, 
however, the propeller should never work below the baseline and at the same 
time should be so situated that the blade tip will clear the bottom of the ice 
at all conditions of trim. * 

t 

Professor Laurens Troost conducted a series of open water tests on modern 

I 

propeller forms in Holland, ^he results of which are published in the paper 

O'*) 

Open Water Test Series with Modern Propeller Forms. From these tests, a series 
of propellers was arranged and plotted on a new m. - ^ coordinate system which 
greatly facilitates propeller thrust calculations in the preliminary design stage. 
The constants involved are: 
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In the icebreaking problem, we can determine the force necessary to over- 
come ice of a given thickness (see section F-2) . Since the force which can 
be developed by an icebreaker at its bow is a function of the propeller thrust 



available and the kinetic energy developed in ramming, we must know the thrust 
which may be developed by a well designed propeller so that we may determine 
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the percentage of thrust developed which is applied to the development of 
ice breaking force. Consequently, if we consider the Troost data and a pro- 
peller design speed of zero knots we may reason as follows: 
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where P - propeller horsepower 

- x shaft horsepower 
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~ 7- 7a 

- mechanical efficiency ~ 98# 
m hull efficiency * 

y p s open water screw efficiency 
y r z relative rotative efficiency 
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If for the density of sea water at approximately 29 F. (when in equilibrium 
with sea ice) we assume 
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Now in the bollard condition it is found that the thrust deduction value drops 
to almost zero. If F 0 is the towrope pull and To the propeller thrust in this 
zero ship-speed condition, it is mostly found in model experiments as well as 
in full scale tests that for single screw vessels 
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Theoretically this should be zero, but the propeller sets up a flow around the 
stern such that the zero value is not quite attained. Hence the towrope pull 
obtained at zero speed is not equal to the thrust developed, but rather, must 

0,-M 

be modified to account for the above phenomenum. Since this thrust deduction 
is a function of water column flow along the hull we may assume without too 
much error that for a twin screw vessel the value is one half that of a single 
screw vessel, hence 

, , y 
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Therefore, because of this small remaining thrust deduction factor, the bollard 
pull becomes 

/ = (O. p j> -.77j / 
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Therefore, for a given P and D, T* is proportional to « 1/j end the highest 
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values for thrust T A will be found on the w- diagram for that part of the 



line 
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where is a maximum. 

Then since 
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we may determine the corresponding propeller HPM by substitution as 
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AJ = c y (34) 

In our icebreaker problem, however, we are interested not only in the 
thrust developed at zero speed, but also at various other speeds. As previously 
stated, for the diesel electric unit the power is constant (PHP s const). In 
addition, a large amount of testing has indicated pretty definitely that changes 
in local pitch from root to tip do not seem to have a marked effect on total 

> 

propeller efficiency, and that it is the pitch at about 0.7 radius that governs 
the blade efficiency, provided that the pitch changes at root and tip are not 

// VC / D/f A / Cf 

extreme, thus giving negative or extremely high angles of i-norodone-e. Conse- 
quently, for our analysis, it is sufficient to assume a constant pitch propeller, 
feeling assured that the results will apply equally well to a variable pitch 

blade. Therefore, the value of obtained from the .V- * diagram for the 

* 



point of maximum 



is constant also. Hence, since 
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if we assume values of we can calculate corresponding values of Q, since 
Q, s 53000 P 

3 rr aj 

Knowing Q, we can determine u since 
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and pick off values of & and 0 



With these values we may determine pro- 



peller thrust as 
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at a corresponding vessel speed determined by 
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Very little information exists as to the wake fraction of ice breaking 

c°) 



vessels. Dr. Jansson, in his paper on ice breaker design lists average values 
of 0.14 and 0.18 for two models which he tested and which had stern propellers 

only. Nordstron, Edstrand and Lindgren in their paper Model Tests with Ice- 

(*>) 

breakers reporting work done at the Swedish State Tank, list values of from 
0.123 to 0.105 covering a speed range from 9 to 14 knots for models which were 
designed with a bow propeller but which had the bow propeller removed and a 
dummy boss and covie substituted. Model tests at the David Taylor Model Basin 
for several classes of coast guard icebreakers list a value of 0.235 for the 
wind class breakers, 0.207 for the Ulaetoitaa, 0.205 for the Cactus class and 
0.200 for the Glacier. In the case of the wind class and the Mackinaw, each 
was built with three screws, one forward and two aft (wind class later modified 
to remove the bow propeller and replace it with a blank covering) and since 
the values given for wake fraction for these vfssels did not specify the status 



of the forward screw, and since the presence of a bow screw materially alters 
the wake fraction of the stern screws, these listed values axe not too reliable. 

In addition, it is known that the wake fraction varies slightly with speed 
over the normal operating speed range, but at the lower speeds, as the bollard 
condition is approached, this wake fraction varies more radically. Unfortunately, 
little if any information exists concerning the variation of wake fraction at 
these lower speeds. Consequently, for the purposes of preliminary design, we 
must be satisfied with the selection of a single representative value of wake 
fraction, which we must assume to apply throughout the entire speed range of 
operation. This, of course, admits to a certain inherent error, however, until 
such time as accurate test data does exist we must content ourselves with the 
fact that the error is constant and allows reasonable prediction of vessel 
performance. In this paper, since we are considering solely the case of the 
twin screw icebreaker without bow propeller, and since it is best to overestimate 
the value of this wake fraction, when it cannot be determined accurately, we shall 
therefore assume 

=0.18 

Therefore 

\J = p y ^/Jb J knots (37) 

In regards to the developed area ratio which we shall use to determine 
the specific Troost propeller series to use we need not be too concerned since 
the value is not too critical as far as the preliminary propeller design is 
concerned. It must be remembered, however, that blade area does effect overall 
propeller efficiency somewhat. In an ideal fluid, with no friction, it would 
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not matter what the blade area was since we would get the same lift forces 
from a thin blade as from a wide blade, hence the Barae thrust force and 
torque force. In a real fluid, however, we have a drag force as well as 
lift. The components of this drag force are such as to decrease the thrust 
obtained and to increase the torque required on the shaft to produce this 
thrust. Hence the larger the blade area, the more frictional resistance 
and hence drag. Therefore, a narrow blade is best for maximum efficiency. 
Careful consideration must, however, be given to the problem of cavitation 
since too narrow a blade will give rise to serious cavitation problems. The 
area of the blade, therefore, is normally chosen from a cavitation standpoint 
such that the blade width or area is as small as possible, yet sufficient to 
prevent the peak negative pressure from approaching the vapor pressure of 

water. It is interesting to note that in icebreaker propellers, the full 

«*• 

sections, necessitated from strength considerations, makes it impossible to 
avoid cavitation completely, especially at high loading and low speed. Con- 
sequently, for the purposes of the selection of preliminary characteristics 
the best course of action is to select a developed blade area ratio which is 
characteristic of present day practice. 

When considering the number of blades for use with a twin screw ice 
breaking vessel it may be said that there is very little effect on efficiency 
despite the number of blades. The most common arrangement is either a three 
or four bladed wheel. Some authorities hold that the general rule of thumb 
is that the three bladed wheel is more efficient than the four bladed one. 



However, for moderate or high blade area ratios it is possible to design a 



four bladed wheel which will have efficiencies equal to or greater than the 
corresponding three bladed propeller. Consequently, the choice of whether to 
use a three or four bladed propeller is almost entirely a question of vibra- 
tion and that the number of blade impulses per minute is the deciding factor. 

Of course, since the requirements of icebreaking differ greatly from those 
normally encountered on sea going vessels, there are many theories concerning 
the selection of number of blades. The Glacier and wind class were both 
designed with the*e bladed propellers. Both have shown blade failures. This 
may perhaps be taken as an indication that icebreakers should befitted with 
four bladed propellers. In this case, perhaps, large ice pieces could not get 
in between the blades and cause failure. Others, however, feel that the three 
bladed wheel is best. Each blade of a three bladed propeller is stronger than 
that of a propeller with more blades and the greater gap obtained with the 
three bladed propeller as compared with propellers with more blades, is more 
favorable from the point of view of ice passing through the blades. Conse- 
quently, the choice of number of blades is purely an arbitrary one and subject 
to the opinions of the designer. Despite the choice, however, the problem of 
vibrations must be considered and it must be insured that the blade frequency 
i6 different from criticals or harmonics thereof of the hull or shaft vibra- 
tion frequencies. 

For the purposes of this paper, therefore, we shall consider both the 
three and four bladed propeller. Experience dictates that the developed area - 
disc area ratio normally used in ice breaking type vessels falls into the 



- / O 0 



range of from 0.35 to 0.54 for the three bladed propeller and from 0.50 to 
0.60 for the four bladed propeller. For the purpose of this paper, therefore, 
we will Belect the Troost propeller series B3.50 and B4.55 as representative 
of present practice. 

By applying expressions (33) and (34) to these series propellers we 

4 

may determine the variation of maximum bollard pull attainable as well as 
the shaft HPM necessary to obtain this ballard pull both as a function of 
shaft horsepower available per shaft and propeller diameter. This has been 

done for a range of shaft horsepowers and propeller diameters and the results 

3 , /<=> /K / j 

are shown in figures 15, 16, IS and 19. By the use of this plotted data we 
may, for a known shaft horsepower per shaft and given propeller diameter, 
determine with reasonable accuracy the bollard pull which will be developed 
per shaft and the shaft HPM necessary for its attainment. 

By comparing the results obtained from these two series we may further 
substantiate our previous statement that the number of blades used is not 
very critical in the preliminary analysis. Since for the B3-50 propeller 

7 " i° 

and for the B4-5S propeller 
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we see that the percentage difference in predicted bollard pull between the 
three and four bladed propellers - all other factors remaining constant - 
is only about 2.7^ and the percentage difference in shaft HPM required is only 
about l.Q%. 

We may, in a like manner reaffirm another previous statement of ours. 

We mentioned that the propeller series used is not critical in the preliminary 
analysis so long as that used is typical of those normally found in a particular 
type service. For this purpose we may examine several propeller series, holding 
the number of blades constant, and observe the percentage variation in bollard 
pull obtainable and shaft HPM necessary for its attainment. For the three 
bladed propeller let us consider the Troost series B3-35, B3-50 and B3-65 
propellers and for the four bladed propeller consider the Troost series B4-40, 
B4-55 and B4-70 propellers. Considering these series we find that 
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B3-35 and B3-50 


a 0.94$ 


B3-50 and B3-65 


- 6.33$ 


B3-35 and 33*65 


s 7.23$ 


B4-40 and B4-56 


s 4.63$ 


B4-56 and B4-70 


= 1.11$ 


B4-40 and B4-70 


= 5.68$ 



and the percentage difference betxveen the shaft HPM required to obtain maximum 
bollard pull with the series 



B3-35 and B3-50 


s 5.22$ 


B3-50 and B3-65 


= 10.1$ 


B3-35 and B3-65 


a 16.2$ 


B4-40 and B4-55 


. 1.13$ 


B4-55 and B4-70 


« 2.78$ 


B4-40 and B4-70 


s 4.04$ 



Hence we can see that for the three bladed propeller a variation in developed 
blade area ratio as large as 85$ produces only about 7$ variation in obtainable 
bollard pull and about a 16$ variation in shaft KPM required. Likewise for 
the four bladed propeller a variation in developed blade area ratio of as 
large as 75$ produces only about 6$ variation in developed bollard pull and 
about 4$ variation in shaft KPM required. 
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Consequently, we see that the actual choice of propeller series is not too 



critical and that the results given in figures 3 , /o , / z , and IS for the 
B3-50 and B4-55 propellers are quite satisfactory for preliminary powering 
analysis. 

As previously mentioned, we are interested not only in the bollard pull 
which may be developed, but also in the thrust which we may expect at various 
speeds of advance. Consequently, we wish to investigate the variation of 
propeller thrust as a function of ship speed for a range of shaft horsepowers 
and propeller diameters. To do this, we proceed as previously outlined. 

The pitch diameter ratio is constant over the speed range since we have 
designed for b©llard condition, hence our values of '& and lie on a line 
of constant *°7 d • addition, since we are considering the. diesel electric 
installation, the propeller horsepower is also constant. Therefore, in order 
to obtain an orderly variation of propeller thrust with ship speed we may 

4 

proceed down the constant 7o line of the Troost series charts and select 
corresponding values of </ , > and ^ . From these we may determine our pro- 
peller thrust and vessel speed by our previously developed relationships 
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For the value of Q, we may combine 
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eliminating N to obtain 

p v h r 

0 , * zy.s - - 

If ?/e substitute the value of Q, into (35) and (37) we obtain our desired ex- 
pression for propeller thrust 

7-'= /j-6 l u/j ■’ _ ^ (38) 



and for ship speed 

^ >9 /• ^ --- < 39 > 

This investigation was carried out for the Troost Series propellers B3-50 and 
B4-55 for shaft horsepowers ranging from 4000 to 28000 per shaft and for pro- 
peller diameters ranging from 5 feet to 30 feet. The results of this investi- 
gation are shown in figures 23. to 3£. Interpolation may be used for inter- 
mediate values of shaft horsepower per shaft but for intermediate values of 
propeller diameter it is necessary to use the figure for that diameter closest 
to the propeller diameter desired. So doing will naturally introduce some error 
in thrust predictions, however, the error will be slight and usually insignifi- 
cant as regards preliminary powering analysis. 

If it is required to have more accuracy for intermediate propeller diameters 

/ <> <—(z> 

than can be obtained from figures at to K in predicting variations of propeller 
thrust as a function of ship speed we may proceed as follows. Returning to our 
original expressions 
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we can detemnine from the Troost data the variation of </ ^ as a function 
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of /v ,/j • Hence for a given propeller diameter and shaft horsepower per 
shaft we may determine the value of '/s for any desired ship speed from 
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We can then select the corresponding value of /^y L 6 and determine the thrust 
at the given speed as 

7“'» lbs/shaft 



For this purpose an investigation of the variation of ^*4 as a function of 

was made for the Troost series propellers B3-50 and B4-55. The results 
are shown in figures S3 and 3%. ' 

We can, if we wish, eliminate the use of figures S3 and 34 by means of an 
empirical relationship determined by curve fitting. If we carry this through 

w ® find that for the Troost series B3-50 propeller the straight line variation 

T = c.yffy - &./ 



gives a maximum possible error in of only about 2%. In like manner, 



for the Troost series B4-55 propeller the straight line variation 

*/y -&.//? 

gives a maximum possible error in only about 4#. Hence if we can 

accept these percentage errors in our preliminary powering analysis we have 
for the series B3»50 - • 
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And for the series B4-55 
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The thrust produced by a given propeller can be materially affected 
by the size of hub installed and the thickness of the propeller blades. 

This is especially true in the case of icebreakers which often employ detach- 
able blade propellers with their inherently large hubs and blades of very 
thick section so as to achieve reasonable -root stresses when operating in 
ice. Because of this fact, the values of propeller thrust predicted for a 
given blade from data presented in this section should be modified to reflect 
the influence of hub diameter and blade thickness. 

Let us first consider the effect of propeller hub diameter. Shultz, in 
his paper The Ideal Efficiency of Optimum Propellers having Finite Hubs and 

(/to) 

Finite Number of Blades reports on work done at the Taylor Model Basin on the 
effect of hub diameter on developed propeller thrust. His results are pre- 
sented as plots of the variation of ideal thrust Coefficient ( C± t as a 
function of advance coefficient ( a ) for various ratios of hub diameter to 
propeller diameter. In order to make use of this data, let us consider 
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Previously we stated that for the B3-50 series propeller = and 

for the B4-55 series propeller ^Vi> - *• Substituting for propeller pitch 

and using our selected value of = 0.18 we have 

>7 - ^ ^ 

B3-50 series 

X - £_■» • * /■>" ^ 

/U A 

7 . B4-55 series 

/U <D 

Noting that speed coefficient J” is defined. as 

r- JHL _ ^ 

' h ' /Ob 

we have 

s 0.318 J B3-50 and B4-55 series 

7 / - 1.855 J B3-50 series 

% 

7 = 1.94 J B4-55 series 

Over the normal range of J for icebreaker operation we may now investigate 
the variation of as a function of propeller hub diameter. This was done 

for the three and four bladed propeller and the results are presented as figure/* 

Z3 4 ' io 

35&and 35g. 



The data presented in this section on developed propeller thrust was 
obtained from the Troost propeller series wherein the test propellers used 
had values of hub diameter to propeller diameter of 

m 0.180 3 blade series 

& 0.«167 4 blade series 

The actual thrust developed by a given propeller, therefore, is 

' r ~ - ^ 

where K is a constant representing the reduction in propeller thrust due to 
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hubs whose is larger than those actually used by Troost. For the conven- 
tional solid propeller blades used on icebreakers 

~ 0.21 - 0.22 

For the built up or detachable blade propeller 

^ 0.32 - 0.35 

Due to these increased hub diemeters, the ideal thrust coefficient is reduced 
about 3# with the solid propeller and about 14# with the detachable blade pro- 
peller. Consequently, the thrust data presented in figures 3 , /,? ys ~ 

through 26 must be modified so that 

7^ crv „ t = 0.97 T solid propeller 
7 ^ t ruOL & 0.b6 T detachable blade propeller 
to take account of propeller hub diameter. 

Let us now consider the effect upon developed propeller thrust of varying 
blade thickness. Burrill in his paper Calculation of Marine Propeller Per- 

(9 9 ) 

formance Characteristics considers the effect of varying propeller thickness 
to cord ratio upon the lift and drag coefficients of a normal oc-r/saL pro- 

peller blade. By means of this data we can estimate the effect of blade thick- 
ness upon developed thrust. 

It can be shown that the propeller efficiency may be expressed as 

//> r % 

where 

/- z \ 

= blade efficiency a 

; / / * - y 

£ m Drag lift ratio * 1 °J Cl 

- Tangent of hydrodynamic pitch angle 

— , v». 

■ * /TTb 



at zero angle of attack 



( 

t 










1 













/ - 



If we consider first the Troost series B3-50 propeller whose sections all have 
a zero angle of attack, we have 
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We have previously shown that for optimum performance with this propeller 
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Based on data published by Burrill in his paper Calculation of Marine Propeller 
Performance Characteristics , we can determine the value of cr as a function of 
blade thickness-cord ratio for blades of the normal ogival shape. A reproduc- 
tion of this data Eire shown as figures 3/ Si - 

From data on the series B3-50 propeller we can determine that the thickness- 
cord ratio (Vc) at the seven tenths radius is 

*>C r c> . o y b ? 



Therefore, from figure 3 / we can determine that the lift coefficient for the 
blade section is 

C - 2> ■ z j > 



and from figure J ^ 
hence 



the drag coefficient as 
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and 
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For a typical three bladed propeller used on present day icebreakers 
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at the 0.7 radius. This large thickness is necessary from a strength viewpoint 
for operation in ice. In order for the blade section to develop the same lift 
as our B3-50 series propeller, it must be designed with an angle of attack of 

a - - 2 • (Figure J/ ) 

and a resulting pitch of 

- / 7 ~ 2 > ’to** ( - (kJ 

r o ■ Vo 6 ^ 

With this angle of attack and thickness-cord ratio, the drag coefficient of 

j i 

the section from figure 35S is 

C 0 8 0.017 

hence 



£ - 0.0724 



and 

y b c 0.762 

Consequently, by use of the thicker blade section, we have suffered a 11.5^ 
reduction in blade efficiency and hence propeller efficiency. Since thrust 
delivered is a direct function of propeller efficiency this percentage reduction 
will be reflected in the thrust delivered. 

In like manner, for the B4-55 series propeller we have 
r / c . 0.073 

Hence from figure 3 / for zero angle of attack 

^ = 0.355 

and from figure j z. 

m 0.0108 

hence 
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For this propeller, however, we have shown that for optimum performance 

f y 0 m 0.515 

therefore 

n 0.164 
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Based on our above value of e- we have for this propeller 
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There is little data available on the characteristics of four bladed pro 
pellers presently used on ice breakers. However, since from a strength view- 
point a propeller blade is similar to a cantilever beam, and since the blade 
sections of a four bladed propeller must have about the same unit strength as 
the three bladed type, it would seem reasonable to assume that 
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for the four bladed propeller sections at the 0.7 radius. In order for this 
blade section to develop the same lift as our B4-55 series propeller section, 
it must be designed with an angle of attack of 

cv * -1° (Figure j/ ) 



and a resulting pitch of 

Y> a = rr t> * 

- 7. Vi'/ <D 

With this angle of attack and thickness-cord ratio, the drag coefficient for 
the section from figure 3 1 is 



hence 



C 0 = 0.0128 
<r = 0.03605 



and 



/ 4 = 0.862 
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Hence, there results a 2.05# reduction in blade efficiency due to the use of 
the thicker blade section, which reduction is directly reflected in developed 
propeller thrust. 

Based on the above development, therefore, we may say that as a first 
approximation 



A - x. T" 

A c fu/Ox. re 

Where x>, is a reduction factor to reflect the effect of thicker blade sections 
The values of may be stated as 

a 0.885 3 bladed propeller 

s 0.98 4 bladed propeller 

In review, therefore, we note that we cannot use propeller thrust values 

directly as obtained from Troost data and as shown in figures 9 j ' z and 'S' 

2 -> 

through 32. Rather, these values must be modified so as to reflect the thicker 
blade sections and larger propeller hub used on actual icebreaker installations 
Consequently, in all our subsequent work, where reference is made to propeller 
thrust, actual thrust shall be used, where 
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3 bladed propeller 
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Section F 



Mechanics of Icebreaking 
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Part 1 

Theoretical Predictions of the Stresses In an Ice Field 

The nature of the action of an icebreaker in breaking ice is two fold. 

The ship forms a channel of open water in moving through the ice. The end 
of this channel is basically shaped in the form of the bow plan at the water- 
line. As the skip moves into this wedge, it exerts a force in the plane of 

* 

the ice sheet and, in addition, a force normal to the ice sheet as the ship 
rides up onto the ice. The nature of the stresses arising from these two sets 
of forces is different. 

The forces exerted in the plane of the ice sheet cause a state of plane 
stress to be established in the ice which ultimately produces failure in the 
form of long radial cracks emanating from the stem of the vessel. The second 
set of forces, normal to the ice field, produce bending in the ice sheet whioh 
eventually causes the ice to crack radially from the center of loading due to 
tension in the bottom most fibres of the ice. As this loading increases the 
'■ radial cracks lengthen and increase in number until finally a circumferential 
crack, roughly at right angles to the radial cracks, is formed caused by tension 
on the uppermost fibres of the ice sheet. The details of the icebreaking process, 
as well as the overall effects, depend upon which type of failure occurs first. 

For the case of relatively thin ice there is considerable photographic evidence 
to show that the radial cracks are sufficient to cause the ice to fail before 
any bending forces can be generated, hence failure of the ice is due to forces 
arising from the ramming motion of the vessel in the plane of the ice sheet. 

In the case of greater ice thickness, however, it seems pretty certain from 
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experimentation that failure is pri m a r ily due to the foroes normal to the ice 
field. Here what minor cracks are caused by the stresses in the plane of the 
ice sheet are lengthened and increased in number by the tension bending stresses 
in the bottom of the ice sheet, then as loading is increased, circumferential 
cracks are formed due to tension in the top of the ice sheet and the ice fails 
by breaking at these circumferential cracks. This later oase is characterized 
by the vessel riding up onto the ice sheet to generate normal forces and is 
most characteristic for the general arctic and antarctic conditions. It is 
this later case which shall be our primary interest. 

To illustrate the nature of the stresses which arise in an ice field 

f 

during the passage of an ice breaker, it is convenient to consider the ice 
field as a continuous, semi-infinite flat plate resting on an elastic founda- 
tion, which is subjected to a concentrated load (the vessel) applied at its 
free edge. In doing this we may neglect the case of failure due to plan® stresses 
and consider only that case which treats with the bending stresses in the ice 

field as the more critical condition and that upon which the design of an arctic 

or antarctic icebreaker should be based. Now as the vessel climbs up onto the 

ice, the ice shelf will deflect due to static and dynamic loading. Due to the 

deflection there wj.ll be set up a potential energy in the ice shelf. When we 
consider the expression for the deflection we must, by virtue of the dynamic 
loading, consider this a dynamic problem. Now it can be shown that the natural 
frequency of a flat rectangular plate may be expressed as 
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where 

/■yy i - * 1 

a. J i z plate dimensions 

* 1 , 3 

D' z r flexural rigidity of the plate 

u(/- -y' J 

d s weight per unit volume of plate 
in z thickness of plate 

hence if a. ■= h — * as would be true in an ice field then we may assume 

for practical application that 



/ 

Consequently, it would seem reasonable that for our problem, the ice shelf may 

be considered as a flat rectangular plate under a static concentrated edge load. 

eL * s * ,c 

The first study of a flat plate on an eiaetric foundation was by Hertz, 
who gave the solution for an infinite plate under the action of a concentrated 
load. Later Dinuik and Schleicher developed Hertz's solution for a round 
symmetrically loaded plate. Max Wyman of the University of Alberta also treated 

this problem of the round, symmetrically loaded flat plate on the elastic foun- 

(/>vj 

dation and applied his results to the problem of vehicle travel on ice. 

(m ■ ns) 

Westergaard, in the United States, gave the solution for an infinite plate 
loaded by a single row of concentrated forces. In addition, during his work 
on concrete pavement analysis he developed the solution for a semi-infinite 
flat plate on an elastic foundation loaded by a concentrated force at any 
point on the plate as well as for .the case where the concentrated load is at 
a free edge. Utilizing Westergaard' s method, N. M. Gersevanov at the Insti- 
tut*. of Foundations and Substructures studied the unlimited plate, the plate 
limited by one straight and two parallel straight edges and loaded periodically 
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by forces in a single row or in a rectangular grid. Shiparo, at the U. S. 
Army Snow, Ice and Permafrost Research, establishment, treated the problem of 
an infinite strip plate on an elastic foundation under the action of a point 
load using the Westergaard-Gersevanov method but employing Fourier integrals 
rather than series solutions. 

The above mentioned solutions to the semi- infinite plate problem, have, 
on occasion, been considered for the investigation of the icebreaker problem. 



Chenea, Ormondroyd and Naghdi attempted to do so in a joint work in 1950 during 
an investigation of structural stresses in icebreaking vessels. Unfortunately, 
however, they did no more than determine the solution of the flat plate problem 



of ice breaker effectiveness. 

For the purposes of this paper we will review and compare several solutions 
to the flat plate problem to see how effectively they predict ice sheet failure. 
As previously mentioned, we shall not consider stresses in the plane of the 
plate. 

In 1950, Professor Ormondroyd suggested the use of an expression first 
developed by S. Timoshenko. The expression holds for a semi-infinite slab 
supported on an elastic foundation and carrying concentrated loads spaced at 
intervals of "a" along the free edge. When M a" is made very large, the maximum 
stress on the under surface of the slab just under the load becomes 



• and neither they nor any others before them applied the results to the evaluation 
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where 
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c z radius of circular area over load P distributed 
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If c - then ? is a concentrated load and i> > 0.325 A . Therefore, if 

we use the value of poissons ratio previously developed and a value of 
A. - & 4 - &■ oj pos- we ma y substitute and rearrange our ex- 

pression to give 
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as the concentrated vertical force required to cause failure in the bottom 
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most fibers of the slab. 

A somewhat similar solution was developed by Westergaard in his work on 
concrete pavements for air fields. In this work he expressed the stress in 
the bottom fibers of the slab directly under the load as 
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5 ( £ m shape factors of loaded area defined by 

£ ' = i + .A ■£; a.- //. 6 OL.V h' - *■ *** 
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5 - - — . c^o z-&- for load distributed over a circular 
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area of radius r • a 



Q, t,' B t factors related by 
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If we consider the case of the concentrated load, i.e. r= J ® we have 
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If we again use our value of poissons ratio and foundation modulus A_ , 
substitute in (17) and rearrange we have 
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It is interesting to note that if this expression is compared to that quoted 



by Ormondroyd (16) we find that 
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This variation could be due to certain approximations in Westergaards develop- 



P,a 



meat. For instance, the expression stated for the determination of Bi and B 2 

are the result of a curve fitting process and may therefore contain certain 

% 

inherent inaccuracies. It is doubtful if this would account for as much as 
a 44# difference. 

An expression for the maximum stress on the bottom surface of a slab 
» on) 

calculated to that developed by Westergaard is given by Roark as 
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Substituting for poissons ratio and foundation modulus k and rearranging we 



have 
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Each of the above expressions may be used to predict the force P required 
to cause failure in the ice sheet due to tension in the bottom most fibers. 

( / O 3' ) 



Little use, however, has been made of them in the past. Kari, in a discussion 
on icebreaker design gives an empirical expression to use in predicting the 
resistance of ice to rupture. This expression is given as 



p - ? ( h + 2v) (21) 

Since this expression fails to consider the effect of temperature on sea ice 

strength, its reliability is quite doubtful. In a paper on the bow charao- 

((, 3 ) 

teristics of icebreakers, Simpson suggests the use of an expression originally 
developed for a circular flat plate freely resting on a circular rim. This 
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is given as 

P= (22) 

This expression too, is not felt to be too reliable since, although it has a 
temperature dependence in / it does not account for the elastic foundation. 

In addition, the circular plate configuration does not closely enough approxi- 
mate the physioal ice breaking picture. 

It must be recalled that in all of the above expressions for the solution 
of the flat plate problem certain simplifying assumptions have been made. 
Firstly, that the flat plate is elastic within the range of action with a 
single constant -value of Youngs modulus and poissons ratio. Secondly, that 
the thickness of the plate is constant and thirdly that the reaction at the 
base of the plate due to the elastic foundation is a vertical pressure equal 
per unit of area to a constant k times the deflection. In support of these 
assumptions we can say that assuming constant thickness and constant poissons 
ratio is not too unrealistic. As regards Youngs modulus it is true that it 
will vary with temperature and position in the ice sheet. The variation with 
temperature we can account for by selecting a value representative of that 
temperature at which we are working. The variation with position in the ice 
sheet we cannot account for, however, since only the fourth root of Youngs 
modulus enters into the above expressions, it is felt that the resulting 
error will not be unreasonable. The assumption that the ice sheet remains 
elastic within the range of action could cause difficulty since it is known 
that ice behaves plastically. However, as we have mentioned, if the rate of 
loadings is rapid (above about 7.11 lb /** 1 -- je<L ) the ice will behave elastically, 






and such may realistically be assumed to be true during ramming with a 
modern icebreaker. 

In order to better compare the previously given expressions for the force 
P required to rupture an ice sheet in bending, figure ->3 shows a plot of the 
force required as a function of ice thickness for the single case of perennial 
sea ice at the average mean ioe temperature of -12° C. This plot shows a com- 
parison of values obtained using each of our above expressions. 

So as to evaluate the adaptability of these flat plate expressions, let 
us select the expression given by Timoshenko as representative of the correct 
solution and investigate the forces required throughout a temperature range. 

For this purpose, let us seleot the strength characteristics of perennial sea 
ice and using the values of Youngs modulus as given in figure 3 determine the 
force required to overcome ice of any given thickness. These results are shown 
in figure 34 - . 

Now for a conventional type icebreaker it is often quoted that the static 
vertical force which may be developed at the bow is approximately 10# of vessel 
displacement and that under dynamic conditions this value may approach 20% as a 
maximum. If such were the case, the conventional wind class breaker should be 
capable of rupturing over 16 feet of perennial sea ice. Experience, however, 
dictates that such a feat is quite impossible. Consequently an error exists in 
our reasoning. 

This error was brought to light by the results of experimentation, on pack 
ice in the field. It was observed that when the floating ice sheets were put 
under load the ice would crack radially from the center of the loading due to 
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tension on the bottom surface of the ice. As the loading was increased, the 
radial cracks lengthened and increased in number until finally a circum- 
ferential crack was formed caused by tension on the upper surface' of the ice. 

** %'f 

How the theory of the flat plate on an elastic foundation will closely predict 
the radial cracks. However, even after the radial cracks have formed the ice 
sheet still has considerable local carrying ability, that is, the pieces of 
ice shelf will not break off. Thus in our ice breaker analogy, although we 
have predicted the cracking of the ice we have as yet not broken it. The ice 
does not break until shortly after the circumferential crack has formed. Now 
the flat plate theory does not accurately predict the circumferential crack 
because the plate is already cracked. Consequently, we must consider the pro- 
blem of a wedge on an elastic foundation in order to compute where and under 
what conditions the circumferential crack will occur. Note that although the 
formation of a circumferential crack does not determine the ultimate 
capacity of the ice, it does occur very shortly before failure. Thus it would 
seem that the solution to this wedge problem would more closely simulate our 
physical icebreaker problem and allow us to more accurately predict the vertical 
bow force necessary to actually break a given thickness of ice. 

The only useable solution to the problem of a wedge on an elastic founda- 
tion is that developed by Mr. D. E. Nevel of the U. S. Army Snow, Ice and Perma- 

(n?) 

frost Research Establishment. Mr. Nevel considered the case of a narrow, free, 
infinite wedge composed of an elastic homogenous isotropic material and resting 
on an elastic foundation. The term narrow indicates that the wedge width is 
sufficiently small so that tangential bending may be neglected. Model studies 
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have shown that this is true if the angle of the wedge is less than 60 degrees 
as is normally true of sea ice under loading. The term free indicates that 
there are no forces exerted on the sides of the wedge by adjacent wedges and 
although the wedge is not infinite, this case will approach that of a finite 
wedge attached to a plate. 

This problem involves the solution of the differential equation describing 
the wedge, i.e . 
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The general solution of this differential equation was determined by Nevel in 
1958. He then considered the case of the wedge with a concentrated load applied 

at the apex,- i.e. x = o and expressed the deflection of the wedge as 

p ’ rr V 

- — ; — - x — ZU (>) (24) 
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and the bending moment in these upper fibers as 

M -- p * .. l23) 
where ( >) describes a Forbenius power series and ^ its second 

derivitive. Now the bending stresses in the upper fibers of the wedge may be 

described as 



■** = r r ~777 (26) 

b B * ^ 

hence if we can determine the bending moment M we may evaluate the rest; 1 • ; 
fiber stresses. In an unpublished report, Nevel has, with the aid of el nic 
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computers, evaluated the series involved in the expression for M. These results 

M j 

he presents ee a plot of /Px versus the nondimensional distance along the 
wedge X for various values of the nondimensional distance of loading 'T' . 

In order to apply these results to the icebreaker problem we must first 
be able to determine the distance of loading of the wedge 7" - — . Now 

Yf J 

n _ £ h 3 , hence we may evaluate it as a function of the ice thickness 

* ' V /zA. 

h, foundation modulus k and Youngs modulus for the ice at the given temperature. 
In determining r, however, we must make some approximation of the area of 
loading. The physical picture of the icebreaker climbing an ice sheet does 

£ o *J Cf 

approach that of a c o nventiona l loading, however, the loaded area is finite. 

In order to approximate 'Q , therefore, we may use the method developed by 
Westergaard in his treatment of the flat plate on an elastic foundation wherein 
he assumes a circular area of loading described by a radius factor L defined 

as ^ =■ y /• & c w t - d> . (> y*, c < /.? ^ 

La ~ C C > /• /<? / A 



where = thickness of the plate 

q - radius of circular area over which load distributed 
For the icebreaker case wherein we approach a point loading we have 

^ — h “ o . j Z a k 

Hence we may approximate K solely as a function of ice thickness. 

Having determined the nondimensional distance of loading V we may enter 

\ /v7/ . ' 

Novel’s curves and determine the maximum value of /Py - t and its location 

along the wedge x . From this value ^ we may determine 

/yj = 1 3 X t ' 
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(27) 

The vertical force P required to break the wedge becomes 

(28) 

& £ ' 

where ^ 4x is the tensile strength of the upper fibers of the ice sheet. Now 

t 

the only unknown in this expression is the width of the wedge at a distance 
X - / hence the angle of the wedge. Field tests show that when the ice sheet 
cracks the wedges formed are less than 60°. Nevel, in his paper, pictures 
these segments as of approximately 22.5° each. Since the actual angle chosen 
has an insignificant effect upon the resultant value of P whioh must be developed 
by the icebreaker, let us assume wedges as depicted by Nevel of 22.5°. Hence 



^ = 0.398 and 
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Now this is the vertical force necessary to break the wedge. Since the ice- 
breaker climbs up on to the ice and exerts a vertical force on an adjacent 
group of wedges (which we assume to be 8 due to the angle of 22.5°) the actual 
vertical force which it must develo'p to overcome all of than and break the ice 
upon which it is resting becomes 



or 
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(30) 



In order to evaluate this expression a curve of the loci of the maximum values 
of ^ /py - c for values of t' in the normal icebreaker range is shown in 
figure 3 i' . The expression (30) was evaluated for perennial sea ice (salinity 
=r 2 i), normal one season ice (salinity* //7W) and young sea ice (salinity 

* for ice thicknesses ranging from 10 inches to 225 inches as a function 
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of ice surface temperature. These results are shown in figures 36 , 3? and 3& 
where vertical bow reaction force F necessary to overcome the ice sheet 
is plotted as a function of ice surface temperature for the above mentioned 



range of ice thicknesses. 
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Part 2 

Icebreaking Forces and Vessel Impact Velocity 

Our analysis is the result of work done by I. V. Vinogradov and published 

( 72 .) 

in his book Vessels for Navigation in Icebound Waters , 1946. The development 
is based on the assumption that no vertical forefoot exists, if one does exist, 
that it does not strike the ice surface. 

We assume an icebreaker, moving with a known velocity, to strike a uniform, 
continuous ice shelf such that the vessel, due to its inclined bow, rides up 
on to the ice until a downward pressure at the bow is developed of sufficient 
magnitude to rupture or overcome the ice. During the period of time that the 
vessel is climbing on to the ice shelf, the propellers continue to turn and 
exert a forward thrust which, at the *■"*»' , is equal to the bellard pull. 

In general, when the vessel is breaking ice of sufficient thickness to necessi- 
tate ramming, it may reasonably be assumed that at the instant the ice collapses 
the vessel's forward speed is, for all practical purposes, reduced to zero. It 
is true that since the ice just forward of the bow is ruptured and since prac- 
tically full ballard thrust is being developed, the ship will have a tendency 
to continue its forward motion. However, personal experience with heavy polar 
ice indicates that this forward motion, if it exists, is very limited and will 
be reduced to zero upon contact with the still unbroken ice sheet - a distance 
of only several feet. In fact, icebreaker practice in heavy ice normally 
dictates that the engines be reversed at the instant of ice collapse so as to 
insure that the vessel will not become beset due to the pressure often exerted 
by broken slabs of ice wedged between the unbroken ice sheet and the vessel's 



hull. Consequently, for the purposes of our analysis of the polar ice- 
breaker we may safely assume that forward, motion is zero upon collapse of 
the ice« 

■yhe quantity which we wish to determine in this analysis is F, the 
maximum value of vertical bow reaction force developed at the stem of the 
vessel. Since this maximum value of F is a function of the distance traveled 
by the vessel on to the ice, it will occur at the instant the ice collapses. 
Therefore, only a study of the dynamic forces leading up to the collapse of 
the ice shelf need be made and no consideration need be given to the downward 
drop of the bow which occurs when the ice gives way. 

The easiest approach to this problem is through a consideration of the 
energy forces involved and the principles of the conservation of energy. The 
energy expended by the vessel is its kinetic energy plus the energy of the 
propeller thrust acting through the distance of vessel travel. This total 
energy is dissipitated by impact of the bow on the ice shelf, change in vessel 
potential energy due to its being raised and changed m and by the fric- 

tion loss due to rubbing of the stem against the ice shelf. In equation form 
this may be expressed as 

r 

< + ^ = *S + f-y f «-> (42) 
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where 

♦ 

Ei r kinetic energy of vessel at instant of impact 
E2 * energy derived from ships propeller 
E3 m Energy dissipitated through impact 
E4 - change in vessel potential energy 
E5 s energy loss due to friction 
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Let W represent the displacement of the vessel in tons and v its velocity 
in ft/sec at the Instant of impact, with the ice shelf, then 

- jj xr '' ft tons (43) 

Next consider the energy delivered by the propellers as the vessel rides up 
on to the ice shelf. During this vessel motion there is a reduction in mean 
draft, designated ^4 , and the ship assumes en angle of t£> m t designated 
£\<P . The distance from the point of contact of the stem with the ice shelf 
to the longitudinal center of floatation (LCF) at the original waterline 
designate as S. The stem of the vessel is sloped at an angle °( from the 
horizontal. During the change in mean draft we may assume that the LCF moves 
in a direction parallel to the stem since the error resulting from such an 
assumption will be very small due to the small ^4 . 




Notice that the assumption is also made that the upper surface of the ice shelf 
coincides with operating waterline of the vessel. This, in reality, is not a 
true representation since there is a certain thickness of ice which is above the 
surface of the water, hence contact is first made at a point on the stem 
slightly above the original waterline. However, this difference is not enough 
to introduce any noticeable error. Consequently, from the instant of first 
contact with the ice until the ice collapses, the linear advance of the vessel 
is 

4 = 4 , f 4*. 
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and. since 



we have 
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therefore, combining fie have 
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If v<e let T represent the thrust in tons which is exerted by the propeller 



during this travel we have that 

- T - ( +- a (44) 

Since we are interested in the maximum vertical bow reaction force F, it 



would be better to express Eg as a function of F. Now if we assume that 
the change in draft, &A, is small (seldom exceeds 10$) the properties of 



the waterplane will not markedly change, hence 
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In addition, the change in trim resulting from impact with the ice is due 
to the re/~>sn/*-><f movement FS. 




If 4^ is small (normally it is less than 5°) we may say that 
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A<p — 

6^J Cy 

Substituting for and in (44) we have 
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Now waterplane area is equal to the product of length, breadth and waterplane 
coefficient, i.e. 



A, 
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and vessel displacement is the product of length, breadth, draft, block 
coefficient and specific weight of sea water, i.e. 

U) z ^ * 

a a /a> 

If we substitute these expressions for /? w and in (45) and introduce the 
non-dimensional coefficients JL x and 4 U • where 
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we have 
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Note that in establishing we make the assumption that longitudinal mete- 
centric height is essentially equal to the distance between center of buoyancy 
and the longitudinal metecenter. 

By the theory of Impact, when two bodies collide normal to each other, the 
loss of kinetic energy involved may be expressed in terms of their relative 
velocities and a constant for the materials called the coefficient of restitu- 
tion. Now in the case of our vessel, the stern does not collide normal to the 
ice sheet since the stern is inclined at an angle c\ with the horizontal. The 



( 



I 



I 



















' / & & 



velocity component normal' *‘to the ice shelf, therefore N Is a/- ^ . Hence 



is the vessel velocity at the instant of contact with the ice shelf, and if, 
as before, we assume that the vessel comes to a stop as the ice collapses, 
the loss of kinetic energy due to impact may be expressed as approximately 



where e is the coefficient of restitution between steel and ice. As far as 

is known, no actual data exists in the literature for the value of this 
coefficient of restitution. Consequently, in order to determine its approxi- 
mate value simple experiments were conducted, the results of which are recorded 
later in this section. For the present, however, we know that it varies between 
values of zero for perfect inelasticity to one for perfect elasticity. 

We can readily see that our elusive vertical reaction force F is a variable 

. - 

which increases from zero at the instant of contact to a maximum value F at 
the instant the ice collapses. Now if we consider the total change in vessel 
potential energy we have that the total rise of the point ’on the stem which 
first contacts the ice is equal to the change in vessel draft plus the rise 
due to the change in vessel re/^7 . Since & is normally less than 50 we have 



The gain in vessel potential energy, therefore, is equal to the vertical bow 
reaction force developed moving through the vertical distance'h, or 
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But earlier we said that 
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Therefore if we substitute ana integrate we have 
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Again using the above substitutions for F and F 3 plus our previously given 

definitions of 4 » £ > ^ and we have 
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Finally we consider energy dissipitated through friotion of the shell 
plating sliding across the ice. The value of the frictional force is equal 
to the product of the coefficient of sliding, or dynamic, friction for steel 
on ice and that component of bow pressure which acts normal to the ship' s hull 
plating. The resultant frictional force, , acts in a direction parallel 
to the stern of the vessel and may be assumed to be evenly divided on either 
side of the stern. 

The energy dissipitated by this frictional force is equal to the 
force moving through the distance the vessel travels in a direction parallel 
to the stern as a result of the change in vessel draft and re/,n , i.e. 

where 
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Therefore 
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Before integrating (52) we wish to know the expression for 5^ . If we disect 
or pass planes through the forepart of the vessel normal to the stem, the 

P 

sections for the submerged portion of the vessel will be in the form of wedges 
with flat or slightly convex sides, i.e. 





Since the vessel is being forced against the ice by its kinetic energy and 
propeller thrust, the ice will be deformed in the form of a wedgelike groove 
in which the bow slides when climbing on to the ice. Pressure is developed 
normal to the faces of the groove and friction forces along the faces of the 
groove and in the direction parallel to the stern. If we designate as 

the resultant force acting normal to the stem then the force acting normal 
to the plating on each side is 

C-<V3 , j 



resulting in a friction force 

where f 0 is the coefficient of sliding or dynamic friction of steel on ice. 
The total friction force, therefore, is given as 
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The value of the force^is related to the vertical bow reaction force and the 
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propeller thrust euch that 
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substituting into (52) we have 
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Again using our expressions for F and F 3 , i.e, 
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integrating and re-substituting for dJ and dc/> we have 
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Using our previously given definitions for A, t A. ^ ^ 



we have 



Isj 



f +r±> 
L <~<r>A 



' C ^~ <* 



] 



S- 7~ ~ 7~ o AJ J 



(54) 



i*o L lot./ ■ ^ y 

If we now substitute the expressions for Ex through E 5 back into (42) and 



rearrange terms we have 
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this becomes 
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Solving for the vertical bow reaction force F and noting that our solution 
must be real and that have 

^ ^ ' /7-^- t££? /^« c-^ -] (66) 

In order now to make use of expression (56) in out vessel evaluation, 
it is necessary to decide upon the value of propeller thrust (T) and vessel 
impact velocity ( *-) to use. T is the thrust acting during the period in 
which the vessel is climbing on to the ice, hence it is a function of ship 
speed and varies from some value T at the instant of impact to a maximum of 
bollard thrust T when yessel forward motion stops. The vessel impact speed 
«r is a function of the vessel mass and the propeller thrust available to 
accelerate the vessel during the ramming period. Let us first investigate 
the expression for ship impact speed v. 

When a vessel accelerates through the water, there is a certain amount 
of water which accompanies the vessel and which must also be accelerated. 

This apparent increase "in vessel faiass - or added mass - was fij;st investigated 

* 

by Froude who first realized the concept and attempted to estimate its effect 

on ship motions. Attwood and Pengally, in their book, Theoretical Naval 

(9V 

Architecture , assume a value of added mass equal to 20% of the vessel displace 
ment. Sinoe such an assumption does not consider the effect of vessel finenes 
or form in terms of its length, breadth and draft - which effects may be pro- 
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nounced - the results are not entirely reliable. Consequently, to more 
closely approach a realistic value for this added mass tie must consider 
the hydrodynamic problem of the motion of a solid through a fluid. 

When a vessel moves in calm water, the inertia forces experienced in 
accelerated motion are 



where the first term is the inertia force due to the ship itself and the 
seoond term summarizes the influence of the surrounding water. The assump- 
tion here is that the fluid is perfect, which is reasonable since as far as 
is known, viscosity has only a slight effect on the values obtained. 

The theoretical solution of the problem of computing added mass was 
first presented by Lamb in his Hydrodynamics . The application of this con- 
cept to Naval Architecture did not occur until relatively recently. For 
ships, exact mathematical solutions are not available, however, approximate 
solutions, sufficiently accurate for practical purposes, can be carried out. 

A comprehensive treatment of this subjeot is given by Weinbliyn and 



St. Denis in their paper On the Motions of Ships at Sea . In this paper, 



the authors present the virtual mass of a vessel in terms of £r©a-dimenslonal 
inertia coefficients wherein for linear motions the additional mass is 
referred to the mass of water displaced by a circumscribing ellipsoid. For 
linear motions along the longitudinal axis of the vessel, the inertia co- 
efficient^, may be expressed as 
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where 



added vessel mass - 
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a . s one half ships length at LWL ■ L/2 - ft. 
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b s one half ships beam at LWL * B/g _ 
c b vessel draft = H - ft. 

Therefore the added vessel mass may be expressed as 

w ' = j &' rr pi. d# 

Values of the inertia coefficient A may be obtained from figure U which 

shows the variation of 4' as a function of the ratios B/ffl and tyfi. Based 

on the above expression, therefore, the total vessel mass which must be 
accelerated may be given as 
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Also, since block coefficient is defined as 
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our expression for vessel mass may be written as 
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For extremely long fine vessels, the L/b ratio becomes large and, con- 
sequently, A approaches zero. In this instance the second term above 
approaches zero and vessel mass is given solely by 
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For the icebreaker type vessel, however, which is quite short and fat, we 
must investigate the value of A more closely. If „ e consider data from 
approximately forty past and present icebreaker type vessels designed for 
ocean service we find that the value of the ratio &'/ varies from about 

0.091 to about 0.153. Based on this range of values we may select a value 
Of 
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as typical of the polar icebreaker type vessel and representative of present 
practice. Using this representative value of ^/ c and assuming a density of 
sea water at 29® F. of 1.995 lb sec 2 / ft 4 we have for total vessel mass which 
must be accelerated 

4 

Z 7 r- r-/ ( 57 ) 

It is interesting to note that the assumption given by Attwood and Pengally, 



ie. 
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is not too bad an elementary assumption for the icebreaker type vessel. 

Knowing the vessel mass which must be accelerated, we may now proceed 

to determine vessel impact speed ^s- . Let propeller thrust acting be given 

\ 

by 7"''and ship resistance at Je*s*ae speed which is below design maximum speed 
be given by ^ . Then the force which is acting to accelerate the vessel in 
open water is 
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therefore 
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Acceleration is defined as 
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where s' is the distance in feet ever which the vessel accelerates. 

The total resistance of a vessel is a function of ship speed and the 



common assumption is that resistance varies as the speed squared. This 
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assumption, however, does not hold too well for the icebreaker type vessel. 

In order to approximate this relationship more closely we investigate the 
variation of total vessel resistance plus appendages for the Glacier and 
wind class vessels for which some data is available. This variation is 
shown in figure up to a ship speed of 12 knots. We are not too interested, 
at present, in the variation beyond this limit since we are considering only 
the velocity to which the vessel will accelerate during a ramming charge, 
and this velocity will seldom, if ever, exceed about 8 knots in open water. 
Also shown on figure <Y° is an empirical relationship between total ship 
resistance and ship speed which approximates the typical icebreaker variation 
pretty closely. We shall assume this empirical relationship in our develop- 
ment rather than the more common y . Therefore 

y ~ - jjf 4 t /,4/ 



and 



a'- — — — (58) 

Tbe propeller thrust which may be developed by a given propeller with 
a given shaft horsepower is a function of vessel speed. This variation for 
the Troost series B3-50 and B4-55 propellers was developed in section E-2 
and was given as approximately 
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If these expressions for r~' are substituted into our expression (58) we may 



integrate and solve for vessel impact speed . The result, however, is so 
unwieldly and difficult to apply as to make it useless in preliminaiy design. 
Consequently, in order to circumvent this difficulty, an assumption is in order 
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If we visualize an icebreaker which is backing down in order to oommence 
another oharge at the ice sheet, we find that after having backed down about 
a ship length to u ship length and a half the engines are reversed full power 
ahead so as to accelerate for the charge. The actual time required for the 
shafts to reverse direction of motion is really unimportant to us here except 
that by the time astern ship motion is stopped and forward motion begins, the 
propellers are turning over at full power and we may, with reasonable accuracy 
assume that the thrust being developed at the commencement of th,is forward 
motion is equal to the ballard pull. Once forward motion begins, it is true 
that the propeller thru3t developed will decrease as a funotion of vessel 

acT/vCr 

speed. Consequently, the thrust during the period of acceleration 

may be visualized as a constant and equal to the arithmetic average of 
bollard pull and thrust at impact velocity. When breaking ice by ramming, 
hov»ver, personal experience dictates that the impact velocity will seldom 
exceed 8 knots, and in most instances will be much less. Therefore, the 
actual thrust decrease over this speed range as 3een from data in section 
E-2 is relatively small. For the present, therefore, let us assume that the 
thrust acting during acceleration is constant and equal to the bollard pull 

t 

/ 

~T 0 . If we do this rather than use tne average thrust acting, the percen- 

tage error in thrust acting becomes 
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In section tf - z. we developed the relationships for 7^ and t for both 
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the series B3-50 end B4-55 propellers, i.e. 
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(33B) 

( 41) 



Therefore for the B3-50 propeller by substitution 
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and in like manner for the B4-55 propeller 
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In order to evaluate our percentage error we oould investigate the range of 

p '/i / «/. 

variation of the ratio /o for past and present icebreakers. This ratio, 
however, has little physical meaning and its trend would be hard to predict. 
Consequently, let ub introduce the ratio T /&- which is a propeller loading 
factor in psi and has physical significance. Since we wish to determine the 
expression for ^ '^/^i in terms of our propeller loading factor we return to 
our expressions for propeller thrust from section E-2. For any given shaft 
horsepower, the ratio ? / D must be constant with ship speed, hence for 
convenience let us evaluate it at the b^dlard condition. From section E-2 
we have, 
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Sinoe the difference between the three and four bladed propeller is so small 

(about 1.5J&) we can use an average value expression to represent them both, 
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Therefore, percentage error in T becomes 
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Since the maximum error will occur at the impact velocity, and since this 
velocity will seldom exceed 8 knots (13.5 ft/seo), the maximum percentage 
error will be approximately 
% error 
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From an Investigation of past and present icebreaking vessels we find 
that the ratio of ^ / cT per shaft varies from about 490 to 765. This is a 
considerable variation, however, we may evaluate present and possible future 
trends in powering so as to arrive at a reasonable value of 77 /o v which 

will be representative of future designs. The value of 490 is characteristic 
of river icebreakers and the smaller, lower powered arctic icebreakers up to 
and including the Mackinaw and wind class vessels. These breakers, as regards 
polar service, are somewhat outdated and underpowered. The higher value of 
765 is representative of the Glacier and several of the more recent large, 



foreign built, polar icebreakers and typifies a trend through the gears toward 
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higher shaft horsepowers per given diameter propeller. 

Now as regards the so called first line polar icebreakers, it may be 

* 

inferred that future designs will go to even higher shaft horsepower so as 
to increase developed thrust and icebreaking capability. Such an inference, 
however, does not seem to be really justified if we review the problem a 
little more closely. In any icebreaker it is desirable to have as large a 
propeller diameter as is possible so as to obtain maximum thrust. Due to 
limitations in water depths in the areas of operation in polar regions plus 
the desire to provide maximum propeller tip immersion so as to minimize pro- 
peller damage by contact with floating ice, there is a practical limit for 
the draft of a polar icebreaker. This necessarily limits the maximum useable 
propeller diameter. This limit becomes somewhat evident if we examine present 
trends. The Glacier, a twin screw vessel with 21000 installed shaft horse- 
power has a full load draft of about 28* 7" and a 17' 6" propeller while in 
contrast, the wind class vessels, all twin screw (modified) vessels with 
10000 installed shaft horsepower and built about eleven years before the 
Glacier have a full load draft of about 28* 0" and a 17' 0" propeller. 

Further, the atomic powered Lenin, a triple screw vessel with 39200 installed 
shaft horsepower has a full load draft of 30' 3" and a design study by 

0. Oakley of BuShips and S. W. Lank of the U.S. Coast 'Guard, presented in a 

C*0 

paper Application of Nuclear Power to Icebreakers , propose a twin screw 
vessel of 30000 installed shaft horsepower with a 29' O'* full load draft and 
an 18' 0" propeller. Of course, installed shaft horsepower may still be 



continuously increased in future designs despite restrictions and propeller 
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protection considerations. If this is done, developed propeller thrust will 
increase, however, such increases in shaft horsepower would be extremely 
inefficient since the gain in thrust with increased power is severely limited 
by the maximum propeller disc area which can be provided. This inefficiency 
would be magnified by the desire to provide as much propeller tip immersion 
as possible so as to protect the propellers from damage from 'floating ice. 

From the above discussion, therefore, it would seem that a value of 
To ^ s 765 per shaft is more representative of present day practice with a 
possible increase to a somewhat maximum efficient limit of about 1000 for 
future large designs. If we seleot a value of 

/sf- - S S' o p er shaft 

as representative of present and possible future trends in ioebreaker powering, 
we find that about the highest percentage error introduced in 7" due to the 
use of 7T rather than average thrust over the accelerating distance becomes 

# error s f 9.12# 

B3-50 

# error s ♦ 10.2# 

B4-55 

Consequently, for our development of icebreaker impact speed, if we consider 
the thrust acting during the period of acceleration as constant and equal to 

7 "'= o.?/77' (59) 

We shall have simplified our problem with a maximum introduced error of about 
0.5#. It is interesting to note that if powering trends do continue as dis- 
cussed and the ratio of /0 //> v do increase to as high as 1000, the error intro- 
duced will still be less than 1.5#. 

Based on our above development, therefore, expression (58) becomes 
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(60) 



•then we have 
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where 
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The expression for F(K^) cannot be integrated by elementary methods but may 
be reasonably approximated by a method of numerical integration such as 
Simpsons rule. This has been done for values in the ship range and the results 
presented as figure y/ . Based on the values in this figure we may determine 
our vessel impact speed for acceleration over a given distance in open water. 
The value of F(X*) may be determined from (61) as 
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wb if ballard pull is expressed in tons, as 
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Knowing FU 1 ) we can determine i 1 from figure Y/ and the vessel impact speed 



as 
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The vessel seldom, however, operates in clear open water when engaged 
in breaking ice. The water soon becomes filled with pieces of ice varying 
in size and thickness which are the remains of previous ramming attempts. 

This ice filled channel which the vessel generates as it progresses, offers 
considerably more resistance to the passage of the vessel than that resulting 
in open water. 

Consequently, when we consider the vessel accelerating in an ice filled 
channel we must, in addition to normal water resistance, include a considera- 
tion of the added resistance offered by the pieces of ice in the water. This 
ice resistance is the result of several forces acting against the hull as it 
accelerates. The largest of these forces is the inertia force which occurs 
when the pieces of ice are accelerated from zero to the speed of the vessel. 
This force can be expressed in the form 

•p _ r -1— 

where h is the thickness of the ice pieces and K a constant dependent upon 
hull form. In addition to this inertia force there also exists a displace- 
ment force of the ice pieces and a friction force of the ice pieces as they 
pass against the steel hull. Thus the total ice force is a combination of 
forces partly dependent upon the thickness of the ice and partly upon the 
speed of the vessel. 

Since it appears from experimentation that the resistance of ice has 
a positive value even as zero speed is approached, the total ice resistance 
may theoretically be expressed as 

= ( C /- C <r'')8A 
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The value of Ci will, in most instances be very small due to the ability 
of the ice pieces to move as vessel motion begins. Its value will inorease 
as the ice pieces become larger or as the channel becomes clogged with a 
heavier concentration of ice pieces. In either case, however, unless we 
approach the ultimate wherein the ice is solid and not broken into pieces, 
the value of Ci will be small enough so that we can neglect it. In this 
case we may express total ice resistance as 

= C. 

If we generalize this expression for the normal ioebreaker hull form, we have 
as an approximation f 

Zee = /V2-3 

• which experience andt experiments indicate is a fairly reasonable approximation. 

The thickness of the pieces of ice in the broken channel vary in thickness 
from very small to quite large. We must, in order to be able to make reason- 
able preliminary predictions, assume some sort of statistical average thick- 
ness. Since we are considering the case of the sea going polar icebreaker, 
we may assume an ioe thickness of 24 inches as a reasonable average for the 
thickness of the pieces resulting from our ramming operation. If we do this, 
we obtain 

l£ lce r 3 */ 2. o V ^ Lbs ( 64 ) 

Since the resistance is in addition to the total water resistance of the 
hull, we must add it to the values of figure in order to obtain a new 
representative curve of total resistance versus ship speed. This is shown 
in figure V z . In order to more easily handle our ship acceleration problem, 
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an empirical carve of total ship resistance versus ship speed is also given. 
This empirical relationship may be expressed as 

v h<j/ lbs 

or if ship speed is expressed in ft/sec 

^ror, L W lba (65) 

If we now return to our expression for vessel acceleration (58) and 
substitute (65) for vessel resistance, we have 
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we have 
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where 

fix') - 

We again evaluate the expression for FjX 1 ) by means of numerical integration 

so as to obtain the variation of F(X^) as a function of X 1 in the ship range. 

yj 

The results are shown in figure As in the case of acceleration in open 

* water we may determine F(X^) by .« 
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or if b©llard pull is expressed in tons 

FU 1 ) . 20.35 - -—^ rJTT (67) 
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Knowing F(X A ) we may determine X A from figure 4$ and vessel impact speed as 
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( 68 ) 



It should be noted here that the values of bollard pull 7T as obtained 
from figure 9 and are open water values, i.e. free of vessel influence. 
Consequently, the values should be modified by a factor whioh represents the 
influence of the hull so as to obtain the thrust actually developed by the 
propeller when in the "behind the ship" position. The factor we speak of 

here is that commonly referred to as the relative rotative efficiency. 

♦ 

The relative rotative efficiency which reflects the dissiramilitude of 
flow conditions to the propeller when behind the hull as compared to when 
in open water, is a difficult item to determine for the icebreaker type vessel. 

Doctor Schoenben, in a paper on propulsion and propellers lists a value of 

(/etr) 

0.985 for twin screw vessels. This, however, is for merchant ship forms 
and may not be too reliable for the case of the icebreaker. If we consider 
several icebreakers for which test tank data is available, we may obtain 
representative values for this quantity. 

By definition 

.= 

'Tv- Y? 7* 7-i 

The mechanical efficiency, , which considers the losses occuring in the 
shaft bearings and stern tube stuffing boxes is normally about ninety eight 
per cent. Therefore 

The hull efficiency, ^ , is defined as 

'% = — 

Therefore 
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For the Glacier, at design speed this becomes 
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For the Wind class vessels, at design speed it becomes 

- 0.982 

If we try to generalize more completely we find that the ratio of effective 
horsepower to shaft horsepower, or propulsive efficiency, remains relatively 
constant throughout the vessels speed range and for the icebreaker type vessel 

averages about 0.58. Therefore 

o ■ £ e Q_ - - 0 - 1 
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Propeller open water efficiency ranges between about 0.50 to 0.65 at design 
speed, hence we may assume a value of z 0.60 as representative for our 
purposes. In section E-2 we discussed wake fraction for this type vessel and 
determined that a value of W * 0.18 was realistic and reasonable for the twin 
stern screw icebreaker with no bow propeller. As regards the thrust deduction 
coefficient, -t , we must again examine past and present vessels for which data 
is available so as to determine a reasonable value to use. Dr. Jansson, in 
his paper on icebreaker design lists a value of 0.27 as characteristic of a 

* yV 

twin screw vessel with no bow propeller. In a paper by Nordstrom, Edstrand 

and Lindgren, reporting work, done at the Swedish State Tank, values O&0.218 

% 

to 0.245 over a speed range of 9-14 knots are given for models which were 
designed with bow propellers but which had the bow propeller removed and a 
dummy boss and come substituted. Model tests at the David Taylor Model basin 
with several classes of Coast Guard snd Navy icebreakers resulted in values 
of 0.27 for the wind class, 0.27 for Mackinaw, 0.24 for Cactus and 0.27 for 
Glacier. In the case of Mackinaw and the Wind class vessels, each was built 
with a bow propeller, although in the wind class vessel this propeller was 
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later removed and replaced with a covering blank. Values listed for these 
vessels, therefore, are not too reliable since the data did not indicate 
the status of the bow propeller which could materially effect the resulting 
value of i « For our purposes, however, we may select a representative 
value of t = o-z s’ 

Substituting our values of ^ and ^ we have for representative 

relative rotative efficiency 

=/•<£> 7 ? 

Since this is only a representative value and determined from so few data 
it would seem advisable to ignore the possible effects of relative rotative 
efficiency and assume 



Tr - /■* 

Hence, for all practical purposes, the thrust developed by the propellers 

* 

when installed behind the ship is essentially equal to open water values. 

In our expression for vessel impact speed, as well as all others wherein 
propeller thrust is a factor, we may use the open water values as tabulated 
in Section E-2. 

Returning now to our expression (56) for the vertical bow reaction force 



f = XT' + Ja x r" - (t-e'J -j] 

we see that the thrust referred to is that which acts after vessel impact with 
the ice. This thrust increases from that which acts at the moment of impaot 
to a maximum of bollard pull at the instant of ice collapse when vessel speed 
becomes zero. This, in essence, is the same condition as existed during the 
period of vessel acceleration, only reversed. Consequently, the reasoning 
used to determine the average thrust acting during the period of acceleration 
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applies here also and. we may re write (56) as 



o.?/ rr. ] (69) 

where at as used here applies to the impaot velocity as determined by (63) 
or (68). 

Previously we alluded to the term (1-e^) in the above expression where 
e is the coefficient of restitution between ice and steel. It was mentioned 
that a value for e could not be found in a search of existing literature. 

Mr. F. L. Ferris, of BuShips, in a paper on icebreaker design uses a value 
of e : 0.95. This value, however, seems excessively high, so i-fc=w»e ‘'decided 
to perform a simple experiment so as to determine a little more accurately 
the range of values for this coefficient. 

Consider two bodies, one with mass rt, and velocity ^/ ( and the other with 
mass ^ and velocity ^ which collide so as to produce direct central contact. 
Let kr; , be the velocity of nf, end Cthe velocity of * /I -h after the impact. It 
may be shown that 

sir - rfl = e ( - *t. ) 

where e is the coefficient of restitution between the materials involved. 
Rearranging 

-tH - All 

e 

If a body be allewed to fall vertically upon the horizontal surface of a 
large mass, we have a special case of direct central impaot. If the mass upon 
which the body drops is rigidly supported, its velocity at all times may be 
assumed equal to zero. Therefore 
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where^ and sr are of opposite sign to signify opposite direction. If the 
body be allowed to fall through a distance H, then its impaot velocity is 

•W - - 1 / h 

and if it rebounds a distance A , its rebound velocity is 

y\r ~ i/z^ X 







Consequently, 




For this experiment, therefore, large pans of fresh and salt water were 
frozen. A steel ball weighing 5.44 grams was dropped from a known height 
through a guiding glass tube sufficiently large in diameter to pass the ball 
with little possibility of friction loss and the height of rebound was 
measured. In all instances the temperature of the icd ranged from -3° C. 
to -10° C. Due to the fact that a controlled freezing atmosphere could not 
be maintained during actual testing, it was necessary to perform a test 
immediately upon removal of the ice sample from the freezer so as to insure 

rnsiriKtOj 

thet subsequent wettin g of the ice at its surface would not result in 

i 

erroneous readings. This process was repeated a considerable number of times 
so as to insure consistent readings. The result of these tests may be 



* 












- ) ?'/- 



summarized as follows 



<? = 0.189 



for fresh water ice 



e = 0.120 



for salt water ice 



la the case of salt water ioe it is known that the results will be effected 
slightly by the salinity of the water and for both fresh and salt water e 
will no doubt vary slightly with ice temperature. However, it is doubtful 
whether these variations would be of such magnitude as to materially effeot 
the above results. Assuming the above stated values of e , therefore, as 
representative for the impact of steel on ice, we have that 



Since the values of e*" in both instances are very mail, we may for all 

« 

practical purposes assume 

e'' - 0 

Therefore, expression (69) becomes 



u 0.0359 for fresh water ice 



e'" s 0.0143 for salt water ice 




(VO) 



where 




Co-o^j -/ j- a C<rt <x 



(55) 



=■ y 



<^0/3 / o< 




(47) 
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Part 3 



Teasel Advance due to Impact 

When an icebreaking vessel strikes an ice field it can overcome the ice 
by means of a bending stress generated when the vessel climbs up on to the 
ice or it may do so by means of compressive stresses as a result of impaot 
forces and propeller thrust acting in the plane of the ice sheet. The forces 

l»v -THE 

with plane of the ice sheet cause a state of plane stress to be established 
in the ice which ultimately produces failure in the form of long radial cracks 
eminating from the bow of the ship. The bending forces, normal to the ice 
sheet, produce stresses in the ice which ultimately cause failure along lines 
roughly at right angles to the radial cracks. The details of the ice breaking 
process, naturally, depend upon which type of failure occur first. For 
relatively thin ice there is considerable photographic evidence to indicate 
that the radial cracks precede those due to bending. In the case of thick ice, 
with whioh we are concerned here, there is considerable reason to believe that 
the compressive stresses in the plane of the ice sheet are insignificant and 
that the ice is broken solely by means of stresses in the plate due to bending. 

In order to investigate the effect of these compressive forces in the 
plane of the ice sheet let us consider, in an elementary manner, how they 
develop . 

Let us again consider the bow section of the vessel to have the shape 
of a wedge and that the various waterplanes and bow seotions are formed by 
straight lines inclined toward the centerline plane. As we have already 
mentioned, such a simplifying assumption will not lead to appreciable errer 



♦ 






since most sea going breakers do have ?>aterplanes and sections of small 

* 

a t their forward end 






/cf >Hger 






Since the vessels waterlines are symmetrical about the centerline or ^ axis 
we need consider only one half of the vessel. If we take the X and ( axis 
as shown we may express the load waterline in the form 



where 4.may be determined by an investigation of existing vessels. The length 
of waterplane which has penetrated the ice sheet may be expressed as 




which, for our case gives 






a 
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Therefore if we consider the bow section which has penetrated the ice in 
three dimension, we hare 





where A = ^ °< . This restriction is necessary to account for the 

H 

condition when forward motion due to compressive forces is such that A > £ . 

In this instance, the above pictorial representation does not hold since A 

if 

becomes a constant and equal to i the ice thickness. For the present, 
however, we shall consider only the case where h 6 since it is the more 
severe condition and typical of polar conditions. 

Now it is also true that the relationship 

does not hold if forward motion due to compressive forces is such as to cause 
the vertical portion of the forefoot to strike and penetrate the ice. In this 
instance h will again be constant and equal to vessel draft. We need not 

t 

concern ourselves with this restriction, however, since in normal icebreaking 
by ramming the vessel will never penetrate the ice so deeply due to compression 
in the plane of the ice. If such penetration is accomplished, the ioe is thin 
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enough not to require ramming since continuous forward motion can be main- 
tained. For thick ice, where continuous forward motion cannot be maintained 
and ramming is necessary, the energy available to compress the ice will soon 
dissipitate itself in crushing the ice locally and forming a wedge-like 
groove in which the vessel will climb on to the ice. However, we shall not 
belabor this point now, but rather, accept it and in a later section its 
truth will be substantiated. 

Now when the vessel strikes the ice the forces acting cause the ice to 
fail by compression in the r immediate vicinity of the stern so as to form a 

9 

wedge-like groove. If we consider that portion of the vessels side which is 

i 

actually in contact with the ice in this groove we have a triangular section 
such that 




? 9 - 




Therefore the area of the triangle becomes 

Area — r _L \ I~, „ v ( ^ v , »- v »/ , v */ 

^ - Z y y la t i- a. C. / It c - 4 -< _c 

For our oase above, however, 
a= j 

b = |/ ^ 

c = v*j '-oS +U C 

Therefore for one side of the groove 



Area 



= A. j/otav 



'<a- 



but 



h s ^ tzu* o( 



Therefore 



Area per side = J_ t ^ 



A* 



Now if we consider the energy available for compressing this ice as 
well as its dissipitation, we have 



El f Eg s Eg MB 

where 

Ej_ s Component of kinetic energy acting normal to the 
stem and lost in impact with the ice 
Eg s Energy due to the propeller thrust acting through 
the distance of compression 
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E 3 = Energy dissipitated in crushing the ice 
In section F-2 we stated that the energy lost due to impact with the ice 
could be stated as 

~ - Li (48) 



where 



yy i 



vessel mass 



Z 2 Vd7 l*J 

~T 



In addition we showed that 

at* - 3.1 x 77 0,621 ft/ sec 
in open water, and approximately 

at = J ./2 x ' r o °' n v ft/sec 
ice filled water. Therefore 



lb sec 2 /ft 



Ejl = 3 34- A '"77 /- us ^ 

in open water and 

, , ^ / o VS 

J ^ ^ i 6 



(63) 



( 68 ) 



(71) 



(72) 



in an ice filled channel. 

The energy due to propeller thrust acting may be stated as 



^ - n 

where 

j - total horizontal distance vessel advances during ice 
compression - ft. 

Now we know that T' is a variable since, as the ice compresses, vessel speed 
will decrease, end hence, propeller thrust increase. In section F-2 we 
demonstrated how this variable thrust could be replaced by a constant average 
value of 0.91 T , where 77" is the beLlard pull in tons, with little intro- 
duced error. Hence 

Eg s ?o 4 o 77 j ft lb (73) 

The energy dissipitated in crushing the ice may be deduced from the physical 
problem. As the vessel strikes the ice the forces normal to the shell act on 



- 2 o/ ■ 



a very small area of ice, hence the compressive stress is very large and 
the ice will fail and deform in compression. At no time, however, does 
this stress exceed , the compressive strength of ice since as this 
stress is reached the ice will fail and relieve itself by continuously 
deforming and enlarging the area of contact with the hull. This will con- 
tinue until such time as the foroes acting normal to the hull act over an 
area of contact large enough to reduce the applied compressive stress below 
. At this instant the ice will resist further deformation and the 
remaining vessel energy will dissipitate itself in climbing up on to the 
ice. Therefore, we may say that 



E3 ■ Force x distance 




( 74 ) 



where 



J ■ compressive strength of ice - lb/ft 

p 

- area of contact between hull and ice - ft 



- distance normal to shell through which force 



moves - ft 



Now previously we showed that 

l ' _ 

fi t ~ j / ^ CK f o( 

In addition! v»e can determine that 

4 - h >^7 H 
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but 



hence 

or 



h : j ^ therefore 



ck 



J ^ C<*svt c\ 



C>{ (J J * n °< J ^ 

r] v r 

^3 ' J cL ^ ^ ^ j/ a . w xAd, 'V / dA 

<> * ii- J ^ 7 /ETv ^ v 



In section F-6 we show that 



z^t y - ‘/j 

^A^yy C\ 



therefore 



y> 



✓ u- v^ y> - — 



and 









_J C j J 






^ ^ ^ J/' ci V 4-^-C 'V 

/ *Xo~! ^7" r 3-V, ^ - ' 

Now wo said previously that x s ay defines the waterline, therefore 



^ - j - C. 



hence 



U t-M ^ — 



CX- 



> r a. ~t , 



Substituting we have 



_ i< « 



/■”, - —1-2 12 <x 



’ J 

fi*ir 






fi-n 



(75) 



Considering j*tst the case of the vessel accelerating in open water and 
equating (48), (73) and (75) we have 



33*/ A ' 7^ / *co V / 3 ° V° /7 J = 



or 



^ J - (,/Zo -/0>c } 2 A ' P 7^ 



O / Vl 



p/ -z: O 



(76) 



where 



p - 



ci < (X. 



Let (77) 

^ = w * ^ 

where and /v- are unknowns whose sum equals a root of the cubic, then 
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- z ^ y - 



J -b -£'±opJ( V / <r~) ~ /ool- 



lox‘ /* 



*}^** <x =r _? 



(78) 



Since the two unknowns and ^ have imposed on them only the condition (77), 
they are not determined uniquely. Hence we may impose a second condition, 
which, in order to simplify take as 

j ^ -xr — o (79) 

We, therefore, have the simultaneous equations 

j ^ \T ^, / J P — O 

. i , w cw/i / 

-/OOWi />/, - o 7 



Eliminating v we have 



j , o 4 /i 1 J 

yiy _ / O o -V /* ^ -i ^ / 47 ' ^ 



(80) 



Solving as a quadratic 



>- o '•/» 



V J =- ro/ ^ -* ' 77 0 ' ’ y>. v t / ro/ ^ ■* ' ^ /"u~, V 



. V W*. 



<5^. so v p.' 



If we let 



* ' 7Z p V y- 7* p * W 






IT r / ^r ’ 4* P \y — ^ L f o/ UJ X ' / ( 



" pfx'o' 7 P 3 - ^ 



we can show that the roots are 

ir^TT v, 4— ^ 



3 

c/ =r / 



\r r 



where 



J - V / 

^ / O 

j ^ 



-z ' 



- i / 

_ - - _ - z , fj 



Now if we impose condition (79), then the only combinations which will 
satisfy it are 

p- pf ° hf u w /<p y Pr- 

Substituting these pairs of values in (77) we have as roots of (76) 



Z c?V* 



Jv - .j //r 4- -o rtf' 

y— J - 

Ji - M' + oo t #■ 



however, since ours must be a real root, we may take as the solution of (76) 



^ VL T 



o , j y *■ 



- r". ^ ^ y o 



Y AS * ) / L 



-f- iO ( X 
/ 



- c 'n>/ ^.-t ' /7 



o • * ✓*- 



r^.j" ^ r /^ 



^ : 






Substituting the value of P we have 



J* 



^ jTo / v * ' *>- 7 ^ ^ * w 

^ j X<^> c* 



r . . , / w r a ^ v' 

jj?/ ^ /j 



Z l 



7 



<Vu CX 



Z o -/ >» /V J * 
t e^c ^ *<i~**,/ 



y. 



j oi ^ ^ 



r , ^ 1 / c 

/ ^O/ C\J^ /J“ ^>v-w <x 



/ 






' <. oy^ r;_ V " A 

L _i«. A <x 



Now for the value of c f^ we find that little reliable and consistent data 
exists. Dependent upon the investigator, values have ranged from as low as 
about 250 psi to as high as about 1700 psi. For the purpose of this paper, 
however, we shall use values which are the result of work by T. R. Butkovich 

of the U.S. Army Snow, Ice and Permafrost Research Establishment as reported 

C m) 

in his paper Strength Studies of Sea Ice . The values given are the result 
of sixty unconfined compression tests including older near shore ice as well 
as younger ice. From the results of this work we may assume that 



C c 



s 1400 psi s 2U2000 



, 
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>v£ 

Now it is frequently claimed that the ultimate strength of sea ice increases 
as the temperature decreases, although many workers in the field were unable 
to detect any temperature dependence. Butkovich does indicate a 6llght 
temperature dependence. However, reliable and consistent data as to this 

dependence is non existent and what variation is indicated by existing test 

% 

results is very small and we can, for the purposes of this paper, consider 
constant. 

As regards the value of a, if we accept the load waterplanes of the 
Glacier, Mackinaw and wind class icebreakers as typical of modern polar vessels, 
then from figure we see that the slope of the waterline at the forward 

end may be taken as 1.835, hence 

a s 0.545 



and 



fas 110000 

° C 



Therefore, 



i 

cJ 



, ^ /. L </ 1 

0 0 118 77 ^ 



^ /. t ■/ 

& ^o.ojmr^r' 77 ^ 

1- 



^ £cxj - ( 



>. O J X d'j" 77 




Oo w ^ ^ \ ' , 



_/•** L 



~ L ^ 3 Ll $ ° x ‘ 



^ V. . 

" ax y f O • ^ ^ '-'XV 0 



For the polar icebreaker the term 

r I ^ 

• < 5 . o/SS~fT? of ) 

is very small when compared to the term 

W* 

( O. 00 ZZ& ’ '77 /•■'/*' J 

normally over 1000 to 1, therefore, we may reasonably assume that 

/ \ 3 

( 0 • ^ / 8 Sf 77 Co-V 1 cxj O 
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without sacrificing our accuracy. Hence 
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( 81 ) 



when acceleration was in open water. 

If we now consider the case where the vessel accelerates through a 
channel whose water is relatively clogged with pieces of ice and whose impact 
velocity is given by 

- / /a j ' 77 ly }+/,*<. l 6 ®) 

we have upon equating E^, Eg and E 3 

VJ* 6> A ’ y 9 r 1 <? */ j> / 4 j = ^ c\ 

" j 



or 
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where 
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Solving for as before we get 
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Again for the polar icebreaker, the second term above is very small as com- 
pared to the first term, about 500 to 1, hence we may assume that 

f o. o / £> j -y c^>V ^ j <2^ o 

and 
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for acceleration in an ice filled channel 
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In order to observe the order of magnitude of the horizontal advance 
due to compression, let us evaluate the expression (83) as the most typical 
case in polar operations for both the Glacier and wind class vessels. For 
the purposes of discussion we shall assume that the vessel accelerates 
through a distance of one ship length which is typical icebreaker practice, 
therefore 

s' s LHP 

and we find that 

y = 8.06 ft. Glacier 

y s 6.92 ft. Wind 

Now it will be recalled that our development was based on the assumption that 
A ^ h , i . e . 

+ " ^ i 

zr n - ^ x 

Therefore 

h s 4.65 ft. Glacier 

h . 3.99 ft. Wind 

Now practical experience dictates that an ice thickness of 4.65 ft for the 
Glacier and 3.99 ft for the Wind class is hardly such as to require breaking 
by heavy ramming except perhaps in those cases where the ice field is under 
considerable pressure. In the majority of instances, therefore, this is 
relatively light ice. Consequently, for all practical purposes we can say 
that 

A 

and horizontal advance due to compression in the plane of the ice sheet is 
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relatively insignificant. 

Now the question arises as to whether or not these compressive forces 
should contribute to the icebreaking force development of section F-2. To 
an extent they do. The component of the vessel kinetic energy which acts 
normal to the stefn is completely expended in impact due to the extremely 
low coefficient of restitution of steel on ice. This energy acts to deform 

v 

the ice in compression and its complete loss is accounted for in the develop- 
ment above. In addition, the thrust developed by the propeller during the 
initial impact stage also acts to deform the ice in compression. This force 
acting through the distance ^ constitutes an energy expenditure which is not 
considered in the development of seotion F-2. However, this is proper since 
its action also is limited to compressive deformation only. 

To reiterate, ice may be broken by compressive forces, i.e. cracking 
radially from the point of impact of the bow, or by bending due to the develop- 
ment of a vertical bow force when the vessel climbs on to the ice. Dependent 
upon the thickness of the ice, one of these methods will become dominant. 

In relatively thin ice where vessel forward motion is continuous, rupture 
of the ice is primarily by oompressive loading in the plane of the ice sheet 
(radial cracking). In this instance the vessel does npt climb up on to the 
ice. As the thickness of ice increases the vessel may begin to encounter 
difficulty in breaking by compression loading alone in which case the vessel 
will climb up on to the ice sheet until the combined compressive loading and 
bending due to developed vertical bow force overcomes the ice sheet. Here, 
the vessel forward motion may still be continuous, i.e. the vessel is not 
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« 

stopped, however, forward speed is not constant since the bow of the vessel 
describes a motion which is similar in character to a saw tooth wave. 

The problem we are considering, however, is the case of very thick ice, 
where vessel forward motion cannot be maintained continuous, but rather, is 
periodic wherein the vessel advances a certain distance with each charge. 

Hence the process of backing and ramming. In this instance the compressive 
forces discussed above do act, but only during the very beginning of the 
ramming process. These forces compress the ice locally into the form of a 
wedge due to the shape of the vessels stem. During this period some radial 
cracks may develop, due to impact, which act to weaken the ice. However, 
since the rate of loading is normally not too high, large cracks seldom occur 
in thick ice. The local compressive deformation will produce very small radial 
cracks which, though they weaken the extreme edge of the ice sheet somewhat, 
do not materially effect the strength of the ice sheet for any distance from 
the immediate point of contact. Consequently, if the ice is to be overcome, • 
it must be done by the development of a vertical bending force which ruptures 
the ice sheet due to large developed tensile stresses in the top fibers of the 
ice sheet. Now, as long as the vessel can develop forces in the plane of the 
ice sufficient to overcome the ice in compression no climbing on to the ice will 
occur. However, once the distance ^ (developed above) is reached, sufficiently 
large compressive forces can no longer be developed and the ice will resist 
further deformation in the plane of the ice sheet. Consequently, the foroes 
which are still acting at this time (propeller thrust and unexpended vessel 
kinetio energy) will expend themselves in forcing the vessel on to the ioe 
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and developing a vertical bow force to put the ice sheet in bending. Hence, 
for very thick ice the actions of compressive deformation and vertical bow 
force development may be considered two separate and distinct phenomenon. 

First, a partial local deformation due to compressive loading, then the 
development of vertical bending forces which cause ice to break in large 
slabs. Now there is no doubt that the forces in the plane of the ice sheet 
will tend to act in conjunction with the vertical bending force so as to 
increase the stress in the ice due to bending. However, due to the thickness 
of our ice sheet the effects of the compressive loading, once compressive 
ice failure has ceased, are small and our problem is basically one of pure 
bending. Therefore, we are justified in our development of section F-2 in 
considering bending only. In that development, the energy available is the 
vessel kinetic energy plus the propeller thrust acting through the horizontal 
distance moved during the process of climbing on to the ice. The energy due 
to propeller thrust acting through the distance ^ is not considered there 
since this energy is completely dissipitated in causing local compressive 
failure. As for the energy expended, however, we must consider the kinetic 
energy lost in impact since it too is completely expended in local compressive 
defomation and does not contribute to the vessels climbing motion. 

Hence we may conclude that although forces do act in the plane of the 
ice sheet and do cause local deformation, the effect is small and the energy 
developed is soon dissipitated leaving the remainder of the energy to force 
the vessel on to the ice to cause failure, if possible, by bending. The energy 
balance as shown in section F-2, therefore, is proper as shown. 
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Part 4 



Vessel Advance and Evaluation of the Vertical Forefoot 

Let us consider the horizontal distance traveled by the vessel from the 
moment of first impact with the ice until vessel forward motion is stopped. 
Thi6 horizontal distance r' may be given by 

r' s d 4- y 

where 



d a horizontal movement of vessel while climbing on to the ice 
y - horizontal movement of vessel due tc ice compression after impact 
Now in eection F-2 we showed that the distance JL could be expressed as 

cX — A d. C\ t A Co -4 tx j b 



where 



AJl z change in vessel draft due to climbing on to ice - ft. 

^ change in vessel /*"»» due to climbing on to ice - radians 

0 * distance from point of impact on stertt to horizontal center 

of floatation - ft 

01 = angular tase of forefoot 
In addition we showed that 

4a = ft 

3 ^ 
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since by definition 

A. 



■SC> 



to r 

Z ^ Vo 



and if we again introduce the non dimensional quantities JL, and 4.x. where 

5 = a, y. 



Cr^ - 



_ i" <ry l 



we have 
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^ , ££i ^ x 



LaJ 



(91) 



where A is the dimensionless coefficient 

r <y . C 1 

J 



^ r i c - ^ yC c, 






ln section #-l we show that the vertical bow reaction force F can be expressed 

as 

^ c.. >, -r- , o- f)> a y '■-7-'' ^ oo 

tons 



°-»z8 x'-rO + °° r 









In addition if we consider the normal icebreaking condition where the vessel 
accelerates through an ice clogged channel then we showed wrr r.hcn wgd in section 



F-2 that impact velocity could be expressed as 

> si y 
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(60) 



Therefore 
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In section F-3 we showed that the distance the vessel would advance due to 



ice compression J could be given as 

^ __ /• O VS 

[ r , <?• °o o^yv 4 CO ' o 
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(83) 
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Therefore, the total horizontal distance travelled is 



r = 



' Ji*. f 2* 
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In section 0-1 we show that the values of Of end i for optimum ice breaker 
effectiveness are 

0/ X 26° 



/* - 



x 45 s 



hence 



r = 1-26 

j? = 0.557 

we also show that for the conventional polar icebreaking vessel 

d ~ V 

Substituting into (93) above we have 



r* 

L 4 J / 

Now let us investigate the forward ship advance so as to determine the 
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(94) 



probability of the vertical forefoot, if one exists, striking the ice sheet. 

If no vertical forefoot is designed into the stern of the vessel, as is true 
in many modern icebreakers, (especially foreign built) our problem here of 
horizontal advance is somewhat academic with the exception of its effect on 
vessel stability when sitting with the bow up on the ice. This is a problem 
we shall consider at a later point in the paper. 

Assuming, therefore, that a vertical forefoot does exist, then our vessel 
advance consideration becomes of practical significance since from an ice 
breaking viewpoint it is probably desirous that the vertical forefoot not 
strike the ice sheet. Some , authorities hold that such a striking blow against 
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the ice sheet by the vertical forefoot is good and desirous since it tends 
to crack the ice and aid the icebreaking process. This is, in a sense, true 
for light ice or free floating floes of thick ice but for solid continuous 
ice this opinion is not subscribed to by this author since such a conclusion 
does not seem reasonable nor does practical icebreaking experience substantiate 
the claimed advantage. The bow of the conventional arctic icebreaker is 
rather blunt and spoon shaped. In addition, by the time the vessel has run 
up on to the ice sufficiently to cause the vertical forefoot to hit the major 
portion of the original vessel energy has been expended. Consequently, it 
is doubtful if the remaining energy acting through the vertical forefoot 
portion is sufficient to cause any appreciable radial cracks or even be of 
much assistance in breaking ice which has resisted failure due to bending 
stresses in the outer fibers. But even supposing that the ice did orack 
radially as a result of the impact - the icebreakers job has not been done 
since the ice has not broken off. We know from experimentation and the theory 
of ice mechanics that it does not take much force in bending to cause this 
ice sheet to orack, but it takes a considerably greater force to actually 
break the ice off. This theory we reviewed in section F-l when we considered 
the force required to break a wedge on an elastic foundation. Consequently, 
our cracking the ice does not materially aid the icebreaking process. 

Either the ice has already been cracked into wedges by the vertical bow 
reaction force F, but this force has not been sufficient to break the wedges j 
or else the vertical bow reaction force has not been large enough to crack the 
ice in which case it is entirely insufficient to ever break it. In either 
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case, additional cracking of the ice due to the impact of the vertical fore- 
foot will be totally ineffective. If, sufficient power is available in the 
vessel to propell it far enough on to the ice so that the vertical forefoot 
strikes^ it would seem desirable to move this vertical portion further aft or 
remove it completely so as to allow generation of a larger vertical bow 
reaction force F which is effective in actually breaking the ice. It would 
therefore seem that the most desirable arrangement would be either the com- 
plete absence of a vertical forefoot or its location far enough aft to permit 
the vessel to develop the largest bow reaction force F possible with the 
power installed. 

There are two limiting considerations to the above discussion. One is 
that the vessel should be able to extract itself if she becomes beset. Con- 
sequently, the powering of the vessel should be such as to satisfy the con- 
ditions outlined in section F-5 where this problem of vessel retraction is 
considered. Secondly, the vessel must not be so overpowered as to climb up 
on to the ice to such an extent that vessel stability is jeopardized. 
Admittedly the intact stability of this type vessel is inherently very large. 
However, when sitting in the inclined position with the bow on the ice this 
stability can be markedly reduced and care must be taken that the vessel not 
be able to propell itself so far as to become unstable. This problem is 
considered in a later section. 

In regard to these limiting considerations, it might be proposed that a 
vertical forefoot be used to prevent the vessel from advancing too far on to 
the ice thus limiting the vertical bow reaction force which may be developed 
in accordance with the above considerations of vessel stability and retraction 
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ability. Such a proposal is justifiable in that the vessel impact energy 
is a direct function of the impact velocity, and though the powering of the 
vessel may be such as to keep the vessel from climbing too far out of the 
water when accelerated through a reasonable distance (say two ship lengths), 
it is possible that increasing the accelerating distance or accelerating in 
open water could cause the vessel to climb so far as to be unable to retract 
or impair vessel stability. In such instanoes the vertical forefoot would 
act as a preventer. In addition, it might be desired to considerably overpower 
the vessel so that it could more effectively force leads, tow in ice or have a 
higher design speed. In this instance, even if a reasonable accelerating 
distance is used, the vessel could again climb so fax as to be unable to 
retract or impair vessel stability. Here again the vertioal forefoot could 
act as a preventer. In both these instances the use of a vertical forefoot 
would be desirable. Its location, however, must not be so far forward as to 
materially reduce the vertical bow force which may be developed. Consequently, 
if the vertical forefoot is used, it must be located far enough aft as to 
allow the vessel to develop the maximum vertical bow force possible with the 
power installed yet be able to retract in accordance with section F-5 and not 
have its stability endangered in accordance with section 0-2. It must be 
reiterated, however, that vessel energy remaining at the instant of vertical 
forefoot contact is completely lost and useless to the icebreaking process. 

Let us now consider the present day conventional arctic icebreaker so 
as to determine the probability of the vertical forefoot striking the ice 
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sheet 




Now let 



h' = JL 



where Jh is some constant. Therefore 

J2. = // C ' * &*) «* 

Now if v the vertical forefoot will strike the ice sheet with subsequent 

loss of remaining ship inertia energy. If we consider our selected value of 
c* s 26° we may state that the vertical forefoot will strike if 

o s’ V f/~X) r (95) 



Substituting (94) for r we have 
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or if we rearrange terms we have that the vertical forefoot will strike if 



X > C 



c= / - 
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In order to examine these relationships for present and past icebreaker 
designs let us assume that the Teasel is twin screw and accelerates through 
a distance of either one or one and one half ship lengths, i.e. 

s'= /•J'k. 

From an analysis of existing and past vessels we can determine representative 
variations of draft, ^/Dia 2 and r » /W as a function of ship length and also 
a representative variation of propeller diameter as a function of vessel draft. 
These are shown in figures ^ , v? , fo and . Based on these characteristic 
values we may evaluate the above expression for C. The results are shown as 
figure s'i. 

An investigation of past and present vessels which have a vertical fore- 
foot installed reveals that for the conventional arctic icebreaker with a bow 
propeller 

-JL o ■ bo 

and for those without bow propellers 

Jb. ~ ■ V 6 

Based on these values of J}_ we see from figure *2. that for the conventional 

arctic icebreaker using full installed power 

&. ^ c 

therefore when breaking ice by ramming the vertical forefoot will always strike 
the ice sheet. This fact can also be seen from figure r J wherein are plotted 
r' (vessel horizontal advance) and J (advance possible before vertical fore- 
foot strikes) as a function of ships length for the representative present day 
ioebreaker. Here again we see that JL must be almost zero in order for the 
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vertical forefoot not to strike. There is really no sound reason for «• 
ranging so high as from 0.46 to 0.60, and proposed future designs do indicate 
a trend towards reducing this ratio to about - 0.35. However, figure 
indicates that even this value is too high and with present powering trends 
will cause loss of inertia energy by striking of the vertical forefoot. It 
would therefore seem that the complete elimination of the vertical forefoot 
(i.e. JL- o ) or the installation of a very small one would be advantageous. 
Reiterating that the only restriction to such a move is the vessels ability 
to retract and the retention of adequate stability when sitting up on the ice. 

Now let us consider the condition where the vertical forefoot does strike 
and calculate the percentage of energy lost and the maximum vertical bow 
reaction foroe which may be developed. If the vertical forefoot strikes the 
ice sheet then ship horizontal motion due to climbing on to the ice (d) i6 
limited to ( £ - ^ ) where taken for most common condition of acceleration 

in ice channel. 

We have shown that £ could be expressed as 



J = 



(91) 



therefore 



r** co-V a - 3 



J 
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and the maximum bow reaction force which we may develop becomes 
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Again we may consider A - Y as characteristic of conventional polar ice- 
breakers. Substituting this and our selected value of K we have 
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This represents the maximum bow reaction force whioh may be developed if the 
vertical forefoot strikes. This must be compared to that which could be 
developed as expressed by 



where 



- </> ■?/ x ?/ o. ft s j- f ‘^Ll2fL- 1. 

Vj V 
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for the ship accelerating in ice filled water. If we substitute the value 
of at into and using our value of o< » 26° and 3 o 45° we have 
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If, for the purposes of comparison we again make use of the representative 
polar icebreaker size and powering characteristics previously depicted in 
figures *'•»" , 7/ , and »'/ we may evaluate expression (98) and (99) for 
a percentage comparison. The results of this comparison are shown in figure 
S’ 4 . Here we can see the amount of vertical bow reaction force which we may 
develop when a vertical forefoot of a given Jt_ value is installed ( f, ) as 
compared to that which could possibly be developed without the vertical fore- 
foot ( f L ) for a vessel accelerating in our ice filled channel for a distance 
of one ship length. It can be seen that for a vessel comparable to the w±ad 



dates where JL * 0.46 about 4$#- of the vessels potential icebreaking force 
(vertical reaotion F) is completely lost when the vertical forefoot strikes 
the ice. 



In order to better visualize the percentage of vertical bow reaction 
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force which is unattainable due to the vertical forefoot striking the ice 
sheet, let us consider the ratio of '■A. , where again 

/f s Vertical bow reaction force which may be 4eveloped 
if a vertical forefoot exists and strikes the 
ice sheet (98) 

s Vertical bow reaction force which can be developed 

if no vertical forefoot is installed (99) 

f, me 

The results are shown in figure ry" where^/1 fcs based on the vessel accelerating 

for one ships length in an ice clogged channel. Again it is seen that for a 

QlHcK & VS 7° 

vessel comparable to the wjf td -Q l a aa , only about -54^ of its potential ice break- 
ing force is realized before the vertical forefoot strikes the ice sheet. 

From our past discussion it can be seen that the use of a large vertical 
forefoot is expensive in terms of lost potential icebreaking force. Two 
alternatives exist. Reduce vessel powering to a point where full vessel ice 

breaking potential may be realized despite the forefoot. This is not too 

/ 

desirable an alternative since it seriously limits the vessels ability to 
foroe a lead thereby hampering her successful operation in field ice. The 
other alternative is to reduce the size of the vertical forefoot or remove 
it completely. In this instance the vessel may climb out on to the ice so 
as to realize her full icebreaking potential commenserate with her powering. 
This seems the better of the two alternatives with the reduction in vertical 
forefoot size preferred over its complete removal. It must be reiterated, 
however, that care must be taken to insure that the vessel does not climb up 
so far as to endanger transverse stability or that once on the ice she will 
be unable to retract. See Sections F-5 and 0-2. 
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Ab ility to Retract once Beset 

Suppose the vessel has climbed up on to the ice sheet but the ice falls 
to break and the vessel becomes beset. According to the theory of mechanical 
friction, the force which is acting to impede the retraction of the vessel 
is the static frictional force which is equal to the product of the coefficient 
of static friction for steel on ice and the component of the bow force which 
acts normal to the ships hull plating. The resultant mechanical friction 
force Tf acts parallel to the stem and may be assumed evenly divided either 
side of the stem. 

If, for the present, we assume no stepped forefoot, and if we pass planes 
through the bow of the vessel normal to the stern, the sections will be in 
the form of wedges with flat or slightly convex sides. 





Since the ice has been deformed in the fora of a groove in which the bow 
slides when climbing, pressure is developed normal to the faces of the groove 
causing friction forces along the faces of the groove and in a direction 



parallel to the stem 



As the 3hip sits in this beset condition there is a vertical bow reaction 
force A acting. 




The component * * aots normal to the stern, hence the force which acts 

normal to the shell plating on each side is 

^ a 

Z 

resulting in a friction force of 

^ - Of 

^ 2. y J 

Therefore, the total friction force becomes 

y- _ J- s ^ ^oo 
O-crt , '3 

where P r coefficient of static friction 

The force f also acts, but parallel to the stern and tends to over- 

come the friction force holding the vessel. In addition, if the vessel is beset, 
the propeller will be put full power astern in order to aid retract. For lack 

i 

of better data on astern propeller operation we may assume that full bollard 



pull is being developed. 
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Now the component of propeller thrust » acts to aid the vessel extract 

itself since it acts parallel to the stem. In addition, the component 
tends to reduce the normal force hence reducing the force of friction which 
must be overcome. If we consider the propellers acting then the total normal 
force becomes 

C\ — 77 " X 

resulting in a total friction force of 






£ COO C\ *“ *( 



and the force tending to overcome this friction force becomes 
consequently if 



/C j T7 c -00 or f i 



( ^ C<rO 



V 



( 100 ) 



V ^/ 3 

the vessel will be able to extract itself once beset* We may re wTite this as 



rr ^ /T 



Vj - c ^°/3 1 



h 



PC / c<rO 
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3 



(101) 



If we let 



we have 



U - 



•fi — o> 



■r: ^ (102) 

Now L> can be either positive or negative dependent upon whether # <^ro ^ 

is greater or less than . If u is negative, then 

(Tail o< ,-Oro i i 

If we consider the vessel in the static condition at the end of a charge 
without the propellers acting the forces acting are the friction force 

A ^ C^rO 0( 
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tending to prevent retraction and the force 

acting parallel to the stern tending to cause the vessel to retract. Now if 

/=- <* ^ + J _ ^ * 

the vessel will not become beset since the friction force will not be great 
enough to hold the vessel. We may reduce this expression to 

/3 fj 

But this is exactly the condition that exists if ^ is negative. Conse- 
quently, if we can, in the design of our vessel, assure that 

fc* X c "‘ r ' J 7 £ 

we need not worry about the vessel becoming beset at the end of a charge. 
This, however, is not always possible or feasible since the value of £ 3 
is an elusive thing and difficult to predict since it is so dependent upon 
existing atmospheric and ice conditions. 

In section G-l we show that the values of <A and / most conducive to 
maximum icebreaking effectiveness are 

^ ■ 26 ® 

fi = 45 ° 

Consequently, so long as £ £ 0.345 we can design for maximum icebreaking 
effectiveness and still satisfy the condition 

j- j 

However, in section D-3 we showed that £, could vary up to approximately 0.5 
and perhaps greater if temperatures are low, the ice very dry and the vessel 
sitting in the beset condition for a period of time eo as to allow adhesion 
between the hull and ice surface. In such a condition U is positive, i.e. 



II 
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Hence, let us reconsider our expression for the vessels ability to retract, 

Le. 

7T o 

Now in section F-2 we showed that the bow reaction force ^ could be expressed 
as 

& •?' r tona 

where A and e* are the dimensionless quantities 

^ C-o—/" ex' /3 — J~ $ 

J / J~ D cx 

£ __ d!-oO J 7 -d-o-0 0< 

Ccto J V /*> ^ w 



)■& = dynamic coefficient of friction 

In addition we gave expressions for /v dependent upon whether the vessel 
accelerated through open water or an ice clogged channel prior to impact 
with the ice. Since the severity of this retraction problem increases direotly 
with an increase in bow reaction force f and since £ increases as Impact speed - 
increases let us consider the more severe case where the vessel has accel- 
erated through open water. Therefore 

v, J . / 4 '4 ft/ sec (63) 

Hence, when is positive, i.e. 

the condition we must satisfy in order for the vessel to be able to free 
itself becomes 
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( 103 ) 
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As previously mentioned, in section G-l we establish values of X and * 
based on an and J which, would produce the greatest ice breaking effective- 
ness, i.e. greatest vertical bow reaction force F, and still be consistent with 
accepted icebreaker design practice. It does not seem logical at this point 
to alter these values of \ and j so as to minimize ^ since the vessels 
primary task is to break ice and this can only be accomplished by the greatest 
bow force F we can develop. Our only concern here is that the vessel will be 
able to extract itself after ramming under the normally encountered ice conditions. 
Consequently, we shall continue to accept our values of 

* = 26° , 

; i - 45° 

unless we can subsequently show that such values of \ and J will prevent 
the vessel from extracting itself. If such is the case, then we will modify 
\ and j and accept a reduced bow reaction force F so as to increase the vessels 
probable ability to retract. Consequently, accepting the above values of end 
S we have 

/ s 1.26 

t? = 0.55? 

and our desired condition becomes 

— r — ^ • ? s' s 3-0 ° J -k (104) 

' ■» ■> 

/y ( , - r 

Now as we’ve seen, the value of U is a direct function of the coefficient 
of static friction . In section D-3 we stated that the value of £ could, 
for steel on salt water ice, vary from about 0.3 to as high as 0.5. Existing 
experimental data indicates that the value of 0.5 is quite severe and would 
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probably be indicative of a condition where the vessel remains in the beset 
condition for a considerable length of time so as to allow the steel hull to 
completely freeze to or adhere to the contacting ice surface. In normal ice 
breaking this is not the case since the vessel attempts to retract shortly 
after the completion of the ram. Hence the length of time that she actually 
sits in the beset condition, assuming ice conditions are such as to cause her 
to become beset, is relatively short and probably only equal to the time it 
takes to reverse her engines to full power astern. Therefore, to evaluate 
our ship characteristics on the basis of a value of f a 0.5 would be overly 
severe and not at all realistic or consistent with icebreaker operating 
practice. To be more realistic and rational in our evaluation it would 
probably be best if we accepted an average value of 



f , ■ °- 4 

and reserve consideration of the condition where J » 0.5 for those instances 
when the vessel moors itself to the ice by ramming it and remaining beset for 
a considerable period of time. 

From figure ? we see that based on our previously selected values of 

c* and A and £ - 0.4 we have 

/ r -> 

a - 0.061 



Therefore, a reasonable condition which we must satisfy so as to ensure that 
the vessel will be able to extract itself once beset is 
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(105) 



For the very extreme condition where the vessel has remained beset for a 



considerable length of time and the temperature is very low, i.e. j- i z 0.5 
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we have from figure s? 

J - 0.163 

and the more severe condition to that 
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(106) 



Now, admittedly, these two conditions above are rather severe since they 

assume that the vessel has rammed the ice after accelerating in open water. 

This is seldom the case in heavy icebreaking for a channel is soon found, 

and when the vessel accelerates, it does so through this ice clogged channel. 

In this instance the impact velocity v is given by 

ft/sec (68) 

% 

hence the condition we must meet so that the vessel will be able to extract is 

(107) 



, ^ . jf i / 

ir= /■/ 2 - \ 
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// C ' -/■ ri- A oj 

Applying the same reasoning as before as regards ^ and j and based on 
j. } a 0.4 we have 
\ =1.26 
^ a 0.557 
U z 0.061 

Hence the condition w<s must satisfy is 



-r~ 0 . >j'<- ^ A ' 

A> • Oooo £S*f‘ — 77 — 

ff 

or tne more severe condition when £ e 0.5 that 

U s 0.163 

hence 



(108) 
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Now previously we mentioned the probable feasibility of modifying A and 
so as to make O - o , i . e . 
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and tnerefore ensure that the vessel would be able to retract solely by virtue 
of the component of the vertical bow reaction F which acts parallel to the 
stem. We could do this based on a static coefficient of friction of 0.4 
which would then leave that component of propeller. bffiLlard pull acting 
parallel to the stem as reserve power to account for the possibility of the 
static coefficient of friction rising above 0.4. 

If such a modification as described above is made, it might be be6t to 
vary of while keeping J constant since any variation in .3 would have to act 
to increase j and cause the vessels forepart to become excessively spoon shaped 
which is undesirable. Consequently, for our case where we have selected a 
value of - 45° we see from figure ^ that we need increase ex' to approximately 
29.5° to meet the condition. 



and make 0 = o . Now this is not too rash a modification and would only 
reduce the vertical bow reaction F by about 2-3^ which in reality is insigni- 

COULD ”1 > 

ficaht. Therefore, it seems that a modification -efeottid be made and wertt£d,in 




the endjbe beneficial. Hence it seems the best selection waalfl be 



= 45 ° 



= 29.5° say 30° 
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Part 6 



Evaluation of Optimum Flare at Side 

Vessels of the icebreaking type are usually provided with a generous 
flare aft of the stem which is carried back to amidships and normally varies 
from about 10 to 20 degrees from the vertical. There are several reasons for 
this practice. The one most often quoted results from the epic experiences 
of the *Fram" wherein the flare will cause the vessel to lift up rather than 
be crushed when subjected to the extreme pressures due to the compressive 
forces of drifting floes. In addition, the flare and the waterline can be 
laid out so that extensive portions of the buttock lines in the forebody will 
be parallel to the sloping portion of the stem. This is felt to be advan~ 
tageous since it causes the entire periphery of the forward end of the vessel 
to bear down uniformly on the ice so as to make the entire forebody, from bow 
to maximum beam, act effectively to break ice. 

The exact angle of flare which is used, or should be used, is subject to 
controversy. It is agreed, however, that some flare should be used since a 
normal ship with vertical sides amidships, fitted with a cut away bow gives 
comparatively poor ice breaking performance. Such a flat sided vessel will 
break ice continuously until the vertical sides are reached, where all breakage 
will stop and forward progress halt due to the friction between the ice and 
vertical sides. 

In many instances, angles less than the 10° - 20° quoted above have been 
used. As early as 1898, Admiral Makarev of the Russian Navy, expressed his 
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feeling that such a reduced angle produced a more effective vessel. To 

/ 

w 

this end he wrote " I find that sharpness of the bow outlines plays 

a large part in icebreaking and that excessively sloping frames, although 
useful in the initial breaking of ice, cause difficulties. It seems to me 
that with sloping frames ice does not get into the vertical position at 
which it can slide along the side. It only inclines somewhat and its whole 
snow covered surface tightly adheres to the vessel. This produces enormous 
friction which stops advance. Steamships with more vertical sides, to my 
observation, can cope more easily with ice." And he continued, "The more 
sloping is the outline of the frames the more the forepart of the ship . 

approaches a spoonlike shape which, in my opinion, is worthless in breaking 

in) 

ice." These observations were fully confirmed by the improved breaking 
ability of the Russian breaker "Eraiak" after her forepart was altered so as 
to give more near vertical side frames. This same view concerning the slope 
of an icebreakers frames in the forebody has also been expressed by many others 
since Makarev. In the present larger icebreaking vessels the flare angle is 
often small since a large flare angle amidships will result in a certain degree 
of flatness of the sides which causes some shouldering of the buttocks where 
they turn into the bow. The smaller flare angle results in an easier buttock 
curvature and the generally rounder lines so produced seem to result in better 
ice breaking performance with less tendency for the vessel to stick. In 
addition, the form which results due to high angles of flare amidships makes 
it necessary to harden up the bilges so as to provide adequate deck area in 
the engine rooms. This form will prevent the early rise of ice driven under 
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the vessel, causing more of it to come up into the propeller. 

4 

In making the choice of flare to use on a new design, therefore, we must 
consider the effect of this flare on the ability of the vessel to withstand 
crushing as well as its effect on the icebreaking effectiveness of the vessel. 
As regards crushing, it is felt that with the modern all welded, heavily con- 

t 

% 

structed icebreaker hull there is little if any possibility of the vessel 
actually being crushed by an ice field. Consequently, the practical choice 
of flare normally will evolve about its effect on the vessels ice breaking 
effectiveness. 

In order to more properly evaluate this point let us evaluate the 
influence of the separate elements of the vessels forepart on ^ts mobility 
in ice. To do this we shall consider movement in broken ice since this is 
where the effect of the shape of the frames is most critical. 

For simplicity, assume that the bow of the vessel has the shape of a 
wedge and that the load waterplane and< bow sections are formed by straight 
lines inclined toward the centerline plane. Such a simplifying assumption 
will greatly facilitate the solution of our problem yet will not lead to 
appreciable error since most sea going ice breakers have waterplanes and 
sections of small curvature at their forward 
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Whan a vessel of such configuration moves under the action of a force 'f" 
through ice, it is opposed by forces representing the pressure of the ice 
which acts normal to the sides of the vessel. The value of this normal 
force is 

Q = ----- 

where y' - angle between center line plane and load waterplane 
The frictional forces caused by the pressure Q' will be 

h 4 ' 



'y - /o 






where y s coefficient of dynamic friction between the ice and steel 

Let us dissect the bow of the vessel with vertical planes coinciding 
with the direction of force Q, as shown in^figure (planes ab and ac) . 

Force C? may be broken into component forces y t ' and tangent and normal 
to the sides. We may now say that once the ice is broken, the force tending 
to push the ice cakes down under the water so as to give way to the vessel 
is the force opposing this motion of the ice is . Therefore, 

when 

/ 

h 4 ( 110 ) 

the ice cakes pressing on the hull plating will slide down along the plating 
under waiter and give way to the vessel. This may also be expressed as 

£> ^ > ( 111 ) 
When the angle X is large, ice cakes may easily go down under the stem and 
bottom and, as experience shows may even, due to the resulting spoon shaped 
bow, adher to the sides. If this occurs, successive cakes of ice will adher 
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to the first so that the vessels motion is soon impeded. On the other hand, 
if the angle / is small and the tangent $ only slightly larger than the 
coefficient of dynamic friction, then the ice will be forced down, turn on 
its side and give way. 

Now we have been considering the section normal to the ships side. If 
we consider the transverse section at this point 




where o«-> - angle between the centerline plane and the ships side at 
the transverse section 
ck r angular rise of forefoot 

In order to develop a relationship between the angle i and ^ we may note 
that 



b - 
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also 

Therefore 

hence 

In like manner 
also since 

Therefore 



b = c co-o 

C r a £5-^ 



(X $ — <x, oo<j 

^ c 7^2-w 

<3- = cP tSLn ^ 

C = <2. T^a^-f C-O 

C - J ^ 



hence 
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C^J — 
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tfX~ 



If we combine the expressions forz*~ 7 and we have 

7^, X ' z rl 

hence our expression (111) becomes 
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From figure 64 we see that for maximum <? whioh is beneficial to icebreaking 
force 

22 ^ <x <c. 31.5° 

Therefore, if we assume our value of coefficient of dynamic friction of 

h = °- 15 

The condition we wish to achieve is 



1/ p( / 

hence as a minimum condition 






s> . /s 



(117) 
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but 0.0225 is extremely small compared to 1 when we consider the 

range of 22° < * < 31.5° hence the minimum condition becomes 






a 0.15 
^ a 8.5° 



(119) 



We may therefore say that the ideal condition would be when^ is sli^itly 
greater than 8.5° say about 10°. 

Now in an actual design the sides of the transverse section are not 
straight as assumed above, but rather, slightly curved. Our results 




a r <j - 



above, therefore, must be interpreted to mean that the angle between the 
tangent to the side at the load waterline and the vertical should be W s 10°. 
We may therefore conclude that the flare of side which is most condusive to 
effective icebreaking effectiveness is about 10°, 
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Part 7 
■ ■ 

Ability to Break Ice without Ramming 

Consider an icebreaker which is attempting to proceed through a contin- 
uous field of sea ice. Assume that the vessel approaches the ice from open 
water and has considerable momentum at the instant it Just strikes the ice 
shelf. This kinetic energy in conjunction with propeller thrust acts to 
force the vessel on to the ice breaking it by bending. Soon, however, the 
vessels kinetic energy is completely expended and the vessels propeller thrust 
must act alone to continue forcing the vessel up on to the ice so as to develop 
sufficient vertical bow reaction to overcome the ice. If the thrust available 
is such that it cannot force the vessel up on to the ice sufficiently to 
develop a vertical bow reaction force large enough to overcome the encountered 
thickness of ice, then vessel forward motion will cease and the vessel must 
back down so as to accelerate and again ram the ice. 

It seems logical, therefore, that when undergoing a design study of a 
proposed icebreaker it would be desirable to insure that sufficient propeller 
thrust will be available so as to allow the vessel to proceed continuously 
through a reasonable thickness of ice without the necessity for having con- 
tinuously to back down and ram. 

Let us imagine, therefore, that the vessel has just completed a charge 
and has broken sufficient ice as a result of that charge so as to reduce its 
forward speed to zero. As the propellers continue to turn at full power the 
stern touches the unbroken ice shelf and the vessel attempts to continue its 
forward motion through the ice. The ensuing events may be described in equation 
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form as 



E]_ z E2 4 E3 



(120) 



where 



E-^ s Energy derived from the ships propeller 
Eg = Change in potential energy of the vessel 
Eg z Energy loss due to friction 

In section F-2 we showed that the energy derived from the propeller could 
be expressed as 

jT / - r^' 1 ^ ft tons , 
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If we put this in terms of the vertical bow reaction force (F) this became 
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(46) 



where 

iz: f . (47) 

In like manner we showed that the change in vessel potential energy could 
be expressed as 

ft tons (51) 
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and the energy loss due to friction as 



A' - Fr U 

J ^ ^ <3 
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ft tons 



d - /5 

If we substitute the value of Ex» Eg. and E3 baok into (120) we have 



(54) 
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If we rearrange this expression and solve for the vertical bow reaction 
force F we have 

a"- ^ 7-. — * 
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but we have already defined the non dimensional quantity X aa 



Therefore 
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X , J -jp 

V o c-ttV ix. -/- :.<x3 J 



i, 7 s\ d E3 > 

Since the series of events which we are here describing take place after 

the complete loss of kinetic energy, that is at ship speed zero, the thrust 
developed by the propellers in causing the above described motion is prac- 
tically equal to ballard thrust and we may, since the vessels forward velocity 
is practically nil as it climbs on to the ice, assume that in our expression 
(123) __ 

' - K 

Therefore 

_ _ tons (124) 

“ ~ 2 . / 0 \ 

Sinoe the larger we can make F the greater the thickness of ice we can over- 
come, it behooves us in our preliminary design to make ~77 and a as large 
as possible. 7 ^ is a direct function of propeller diameter and installed 
horsepower hence there is' not too much we can do other than insure that we 
use the largest diameter propeller allowed by the vessels draft or beam 
limitations and attempt to have as high an installed horsepower as can 
effectively and efficiently be absorbed by the selected propeller. We dis- 
cussed this problem of power limitation in section F-2 wherein we described 
a propeller loading factor "y o v and discussed how from an economic point 
of view the highest we could reasonably expect this factor to go in the 



foreseeable future was about 



d" 



Z 850 per shaft 
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Hence, in tone, a reasonable limiting approximation could be given as 

7T = if per shaft 

As regards X we do have some latitude of variation in that it is solely 
a function of geometry, i.e. angular rise of forefoot o( and angle between 
the centerline plane and the normal to the vessels side at the bow . 

From figure AX' we see that for any given value of / 3 the value of X 
increases with decreasing tx . It would seem, therefore, desirable to decrease 
c\ . Unfortunately, however, we see from figure £& that a decrease in <X 
beyond the value of A - 26° we have already chosen will materially reduce 2 . 
We show in section G-l that the effect of the value of X upon the vessels 
capability of breaking ice by ramming is greatly in excess of the effect 
caused by varying X and we stated that if any choice were made it would be 
best to make X as large as possible even if at the expense of a reduced 
value of X • It was at that time that we selected x s 26° and d s 45° 
as our optimum values. Consequently, if we were now to dictate a decreased 

/ 

value of so as to increase our vessels ability to break ice in continuous 
motion we would seriously effect its capability for breaking by ramming. 

Since the ramming consideration will dictate the thickest ice that the vessel 
will be able to penetrate it is felt that it should have overriding considera- 
tion and our values of c\ and X> remain at 26° and 45° respectively. 

In like manner, when we discussed the vessels ability to retract from the 
ice once i*b was beset we stated that from the point of view of retraction it 
might be desirable to inorease the value of X to about 30°. If this is done 
however, our value of X decreases from 1.26 to 1.09 or a 13.5$ decrease. 
















I 

i 













- I*? - 



Tbis in turn would mean a 13.5% decrease in the vertical bow reaction force 
developed in (124) above - a considerable decrease. Since vessels of the 
icebreaker type are normally equipped with heeling and tanks to aid in 
breaking the adhesion bond between the ice and the vessels hull if the vessel 
becomes beset and since these tanks are so very effective in breaking the 
bond allowing the reverse propeller thrust to retract the vessel, it is felt 
that our present consideration of the vessels continuous icebreaking ability 
should take precedence and that should not be increased. 

Consequently, it would seem that our previous choice of X * 26° and 
'S z 45° is still the most desirable as regards the ships general operational 
ability in the ice. Based on these values we have 



\ r 1.26 



and our expression (124) becomes 
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(125) 



If we make use of our approximation of propeller loading ratio we may express 
this as 



A = £>. y tons (126) 

Pin 7~3 

Figures o'S to &>& show plates of the approximate number of inches of ice 
through which a vessel may proceed continuously based on expression (125) 
above and as a function of ice salinity and surface temperature. 

There is one other item which we should consider here. That is, when 
selecting the actual value of o( and j to use in a proposed design we must 
in all instances insure that X > o . To do otherwise would make the ship 
incapable of breaking ice except by ramming wherein it made use of its kinetic 
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energy. Consequently, not only do we desire ,) as large as possible, but 

as an actual limit we must insure that 

, , &o-V c* -i>a 5 - - / a 

A - - - ^ 

-j~ b ^ <X / c^y~& 3 

or that 



(X s'* ^ (127) 

In section D-3 we discussed friction and stated that the values of the 
dynamic coefficient of friction for sea ice on steel plate normally ranged 

D /<~t ajS / 

from about 0.1 to 0.2. These values are rather low and for most designers 
past and present, we can show that 

However, this is a consideration and the designer of any proposed vessel should 
insure that his selection of •> and i are not so radical as not to satisfy 
(127). Figure & / shows a plot of co*f<x <w j with the dotted lines indicat- 

i 

ing the probable upper limits of the dynamic and static coefficient of friction. 
As long as the selected values of A and j fall above these two dotted lines 



our above condition is satisfied 
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Preliminary Design Considerations 
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Part 1 

Preliminary Estimate of Vessel Proportions 

Let us consider ‘our expression for the vertical force developed by the 
bow of our vessel, i.e. 



where 
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If we solve this expression for vessel development, we have 
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(84) 



Mow when we investigate the non dimensional quantity A we find that it may 
be considerably simplified for the purposes of a preliminary displacement 
estimate. If we examine past American built icebreaker type vessels we find 
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that the value of Jl. i , is relatively constant, varying from about 1.01 to 
1.034. In like manner the value of for these vessels is also relatively 

constant, varying from about 0.266 to 0.272. We may, therefore, select 
representative values of ^ and » tending in our selection to favor the 

more recent icebreaker types. Doing this, we may say 

4 . = /• e> J / 

&.V- - <£>• zi, } 

Therefore, the ratio 



k. 



4T 



3.83 



which value varies by only about 0.5% from the average value for the vessels 
considered. Based on this ratio, our expression for A becomes 



A - 



1 ,, s.&? 7 
" "cj / 



(85) 



Due to the peculiar requirements for the underwater hull shape of icebreaker 
type vessels, it would seem that the non dimensional ratio would 

remain relatively oonstant with vessel size. This ratio was investigated 
for seventeen past and present icebreaker vessels of the polar type of all 
sizes and it was found to vary from about 0.625 to 0.73. We may select a 
mean value of 




- 0.6735 



as representative of modern icebreaker practice. Therefore, expression (85) 
becomes 

A= o.(>73fJ/+ J ( 86 ) 

This selected value of c V r is further justified if we examine the five 
most recent American built icebreakers. In this instance we have a variation 
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of front 0.626 to 0.683 where the higher value is more representative of the 
larger types such as Glacier and Mackinaw. 

In order to arrive at a value of Co for use in our first estimate of 
vessel displacement we investigate icebreaking vessels built since 1940 
and determine a trend in C. M as a function of speed length ratio. The results 
of this investigation are shown as figure oZ , The mean value, if used for 
our preliminary estimate, will only involve a maximum percentage difference of 
4Jl from extreme values. This percentage difference is felt to be within 
acceptable limits for our purposes. 

A further investigation was made of over thirty five past and present 
icebreakers so as to determine the variation of design speed length ratio. 

The values varied from about 0.849 to 1.20. This seems a wide variation, 
however, closer inspection reveals that for these thirty five vessels, depend- 
ing upon the year in which they were built, the following variations exist: 



Year Built 

1920-1929 

1930-1939 

1940-1949 

1950-1959 



Range of Design 
0.829 - 1.12 
0.849 - 1.22 
0.946 - 1.044 
0.865 - 1.052 



The trend is towards a higher average vessel design speed. Now it may be 
inferred that future designs will go to even higher design speeds due to 
increased powering and hence higher speed length ratios. This is somewhat 
doubtful, however, as explained in section z~~Z- since vessel draft limita- 
tions will severely restrict the maximum size of propeller which may be used 
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to absorb these increased powers. This in essence becomes an economic pre- 
blem of how much power can a given propeller absorb efficiently. This 
coupled with the faot that a definite limitation as to vessel length does 
exist so as to insure adequate maneuverability in the ice field, it seems 
hardly reasonable that the values of design speed length ratio of future 
vessels will ever go materially above u/v ir/ , it is true that in earlier 
vessels values up to 1.22 were used, but these vessels were extremely short 
and light and were able to use much larger propeller diameters per given 
shaft horsepower than is true today, as will be true in the future. Conse- 
quently, it would seem that the values listed for the 1950-1959 period are 
the most realistic and representative of possible future trends. Hence for 
our preliminary purposes let us assume a representative value of 

- 0.978 

which represents an average of the ten most recent icebreaker designs for 
which data was available, i.e. 



Name 


* la r 


Tear built 


hap it an Class 


0.986 


1954 


Labrador 


1.013 


1954 


General San Martin 


1.030 


1954 


Glacier 


0.941 


1955 


Montcalm 


0.899 


1957 


O &£ 


1.052 


1957 


Moscow 


0.956 


Building 
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Name 


Vn: 


4 

Tear built 


Canadian hull 620 


0.928 


Building 


Finnish design 


0.998 


Proposed 


American design 


0.978 


Design Study 



If we assume this representative average value of spaed length ratio 
at design speed, we have from figure 



c 



0.7426 



Therefore, our expression (86) becomes 



C 
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= 4.0 



(87) 



and expression (84) becomes 

(88) 

Let us now investigate the expression A and t? which appear in (88) above 
and which are given as 
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C^ro J c <y4 c* 



CcrO/3 c -«^ 

ccr? 3 

Each expression is a function of the dynamic coefficient of friction between 
steel and ice, , the angular rise of the forefoot, U , and the angle 

between the centerline plane and the normal to the shell at the bow // . 

Let us consider the coefficient of friction first. 

As outlined in section D-S, little reliable information exists concerning 
an absolute value of the dynamic coefficient of friction. It is known that 
it varies with ice temperature and salinity, however, due to insufficient 
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data, accurate predictions of this variation is impossible. We can, however, 
use what data is available in order to arrive at a realistic representative 
value which we shall consider constant with temperature and salinity. This, 
of course, admits to some inherent error, however, such error is small enough 
not to materially effect our results. Section D-3 describes the reasoning 
involved in the selection of this representative value of -p z> which was given 
for salt water ice as 

L = °* 15 

It must be recalled that the above assumption is based on dry salt water ice. 

If the ioe is wet this coefficient of dynamic friction approaches zero as a 
limit and there is a marked increase in the value of V and € as shown by 
the dotted lines in figures / and &> s' . This fact is one of the principal 
arguments given by the proponents of tne bow propeller who maintain that such 
a propeller can he used to keep the ice directly forward of the ship wet, hence 
increase the vessels breaking ability. For polar operation, however, the thick- 
ness of ice encountered and its hardness are such as to cause propeller damage 
which far outweighs the icebreaking advantage. Besides, when the normal polar 
breaker with no bow propeller rams the ice and then backs off at full power, 
sufficient backwash is generated so as to cause that ice surrounding the bow to 
beoome awash, hence partially achieving results similar to the bow propeller. 
Despite the benefit acruing from a reduced coefficient of friction, we shall 
consider only dry ice as the most severe condition and that which should be 
considered in the preliminary design. 
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The angle (b is a difficult one to consider in that it may have a 

wide range of values and still produce a reasonable and effective design. 

„ fi 

Since the value of ft effects solely the friction force developed and since 
for minimum friction force /] should be small, it would be advantageous to 
use as smell a ft in our design as possible. From an investigation of 
twenty five previous icebreakers we see that /j has varied from 39° to about 
68° over the years. Now the smaller the angle 3 , the more spoon shaped is 
the bow section and the greater the area over which the developed bow reaction 
force F must be distributed. In addition, experience shows that an excessive 
spoon shaped bow causes a rapid reduction of vessel forward energy when breaking 
ice due to the tendency of ice slabs to adherp to the bow and be pushed forward 
with the vessel. Hence a smaller flare angle at the bow or sharper stern would 
be desirable. This in essence means a larger value of ftZ . This, however, 
is in conflict with the requirements of icebreaking force in that for maximum 
icebreaker effectiveness it would be desirable to have y and t? as large as 
possible and from figures 6 V and we see that both y and <? increase 

with decreasing J . Consequently, our selected value of ,J must be a com- 
promise, and based on past operational experience of American breakers it would 
seem that a value of 

w 

3 » 45° 

is an acceptable choice. 

Based on our above value of J , we can proceed to select a value for 
the 'Use of forefoot so as to give as large a value of and z? as possible 
so as to optimize icebreaker effectiveness. A review of figures <4 and 
shows that these quantities arc in opposition, i.e. as c\ increases z? increases 



- Z - 



to a maximum then begins to decrease while Y decreases continuously. As 

regards / , therefore, the smaller the value of Y the better for our ice- 
% 

breaking problem. However, there is a practical limit as to the minimum 
acceptable value of Y . If <Y is very small and the vessel climbs up 
on to the ice but the ice does not break, the component of bow reaction 
force F acting normal to the ice surface will be a major portion of F hence 
causing a large friction force and an increased possibility of the vessel 
becoming beset. If c\ is larger, the component of the vertical bcrw raafiter 
force F normal to the ice surface is smaller giving a smaller friction force, 
and the component parallel to the stern is larger increasing the tendency of 
the vessel to retract. Hence, as regerds X our selection of X must be a 
compromise. However, if we consider, our expression for vortical bow reaction 
force F, i.e. 



f - 0.9/ x ' o A / o-PaSx 4 t- ~--Y . — 

A 3» . 

we note that the expression * * ^ is very large as compared to either 
o-9/ x /T or 8 * '71 , therefore, in order to insure maximum F it 

would seem logical to select c\ so as to make j 2 a maximum even if such is 
not ideal as regards A . From figure we see that for -S s 45°, the cx 






for maximum A is 



ex' = 25° 



From figure t,*' we see that for r 45° and for maximum X and the hipest 
possible X we should use an x which is a little bit greater, say 

^ - 26<> 



Based on our selected values of y and J we have 



I 



- Z??~ 



** = 0.557 
X = 1.26 

Hence our expression (88) becomes 
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(89) 
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If we use this expression in conjunction with those for vessel impact speed, 
i . e . 

° acceleration in open (63) 

water 

acceleration in ice (68) 

clogged channel 

we have two expressions which we may use for the determination of a preliminary 
estimate of vessel displacement. 

Previously we stated that in order to commence a design study there were 
certain items of information which should be known. Among these are limiting 
draft or minimum vessel beam, thickness of ice desired to be capable of over- 
coming, mean air temperature for the proposed area of operation during the 
earliest month of planned ice penetration and possibly the desired design free 
route speed. Based on these known or estimated conditions we may proceed in 
our preliminary ship estimates. 

Investigation of past and present icebreaker type vessels with twin stern 
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screws indicate that the maximum useable propeller diameter may be represented 
as a function of vessel draft or load waterline extreme beam. Figures ^ ^ and 
show a plot of these variations. Hence based on either known limiting 
draft or minimum desired load waterline beam we may select a reasonable value 
of propeller diameter. In addition, based on this known limiting draft or 
minimum load water line beam, we can, by means of figures ^"through S& determine 
a first rough approximation of vessel length between perpendiculars and either 
vessel draft or extreme beam whichever was unknown above. It must be reiterated 
that these estimates of vessel length and beam are only rough approximations 
since with limiting draft fixed, the length and beam must reflect desired dis- 
placement so as to result in a reasonable block coefficient. 

In section F-2 we discussed a propeller loading factor T °/^ which we said 
could be used to estimate the powering of a vessel. If we investigate past 
and present icebreakers we can determine a characteristic variation of * 
as a function of vessel length. A plot of the variation is shown as figure s'o . 

Based on our estimated vessel length we may select a reasonable propeller 
loading ratio and hence, estimate the bollard pull which we can reasonably 
expect under present and foreseeable future standards of icebreaker powering. 
Knowing the expected bollard pull per shaft and the approximate propeller dia- 
meter, we may, from figure 9 or /z determine the shaft horsepower per shaft 
necessary to produce this propeller thrust and from figure or '3 the shaft 
RPM necessary for its attainment. If on the other hand, the proposed installed 
shaft horsepower is known in the preliminary design stage, the above approach 
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may be reversed so as to employ figures ? and /z - to more accurately deter- 
mine the ballard pull which will be available. Note, however, that in either 
case, the estimated value of ballard pull will represent that obtainable from 
a Troost series propeller and must be modified as shown in section E-2 to 
aocount for the larger hub diameters and thicker blade sections normally used 

on icebreakers, i.e. 

7~ £. 7 

4 C r<^/ 7 L 7£oo± r' 

where 





r 0.97 


solid propeller 




s 0.86 


built up propeller 




- 0.685 


3 blade d propeller 




s 0.980 


4 bladed propeller 



In order to determine the vertical bow reaction force which we must be 
capable of developing, we make use of the estimated mean air temperature in 
the probable area of vessel operation during the earliest month of expected 
penetration (table 1 and 2) and the desired thickness of ice which must be 
overcome. Based on this mean air temperature and assuming ice with no snow 
cover we may estimate the ice surface temperature with our expression (12). 
Knowing the ice surface temperature and ice thickness we can, from figures 
36 ,37 , or 38 , dependent upon type of ice involved, dbtennine an estimate 
of the vertical bow reaction force F which we must develop. 

If we assume that the vessel will accelerate through one ship length in 
ramming the ice, wjiich is reasonable from practical experience, we can say 
that 
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In addition, based on our representative design speed length ratio estimate 
for polar icebreakers we can estimate our vessel design free route speed as 

\y — 0.978 

We now have reasonable estimates of our expected bollard pull ( rz ) , 
vessel length (L), design draft (H) , maximum DWL beam (B) , design free route 
speed (V), vertical bow reaotion (F) and vessel accelerating distance. 'Con- 
sequently, the only unknowns in expressions (90) and (63) or (68) are vessel 
displacement (W) and vessel impact speed ( at) . If we first evaluate expression 
(90) for various assumed values of vessel impact speed, we can obtain a plot of 
displacement versus impact speed. If we then evaluate expression (68), as the 

most severe case, for the same assumed values of vessel impact speed we can 

' / 

determine a variation of X with impact speed. Knowing the variation of a 

we can determine the variation of vessel displacement with impaot speed from 
our expression (67), i.e. 

2/-eS 

(*'7 

where F(X) is determined from figure /3 for our Values of A . The resulting 
variation of vessel displacement as a function of impact speed can also be 
plotted as before and where the two plots intersect we have a simultaneous 
solution for vessel displacement. This shall be our preliminary displacement 
estimate. 

We should now correct our vessel length and beam approximations so as to 
reflect a reasonable block coefficient. Figure 67 presents a plot of prismatic 
coefficients as a function of speed length ratio and displacement length ratio. 
Based upon our first estimate of vessel length and displacement and our assumed 
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average value of speed length ratio we may determine a reasonable V for our 
vessel. Figure 6 8 presents apLot of block coefficient as a function of speed 
length ratio and prismatic coefficient. Both figure and <=8 are represen- 
tative of current icebreaker practice and suitable as a guide in the selection 
of a representative block coefficient. Since block coefficient 

c e - — 

and sinoe our vessel draft is fixed, we have 

z ^ - £1^ 

c B * 

But for the large polar icebreaker it is desirous that 

L / 8 v yo" 



hence 

and 



L - /^. // 



L\J / 

/ /c 8 /s 



B s 0.2225 L 



We now have corrected estimates of vessel length and maximum load waterline 
beam compatable with a reasonable block coefficient. Based on these second 
estimates of vessel length and beam we can re-evaluate our estimates of 
vessel powering and displacement, which in turn will allow a second estimate 
of prismatic and block coefficient. Convergence toward reasonable and cora- 
patable values of length, beam, powering, displacement, prismatic and block 
coefficient will be rapid and require only about two trial calculations as 
described above. The resulting values shall form our final preliminary esti- 
mates. In addition, since midship section coefficient 



- 



<- 6 _ 

C P 



we can also, at this time, select a reasonable value of midship section 



coefficient either by direct applioation of our selected values of and c a 












Il 



- 



or by use of figure 6 ? which shows a variation of midship section coefficient 
as a function of speed length ratio. 
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Part 2 

Transverse and Longitudinal stability Evaluation 

Icebreakers are normally vessels with extremely large beams. Consequently, 
they normally have large metecentric heights and resulting high intact stability. 
Because of this fact it is normal for the designer not to be too concerned with 

the problem of stability early in the design stages. Unfortunately, however, 

, A 

adequate stability may not be assumed solely by virtue of a high initial met/a- 
centric height. 

Icebreaking operations are such as to make a large initial metecentric 
height desirable so as to insure adequate stability under all operating condi- 
tions. Probably the most serious instance in which transverse stability may 
be Jeopardized is vihen the vessel has completed a charge at the ice resulting 



in her bow being raised a considerable distance out of the nrffter on to the 
ice. Such a situation is made even more serious if there is considerable iced 
top A* m . i 

Consequently, it would seem desirable to be able to investigate this pro- 
blem of stability early in the design stages while vessel characteristics and 



parameters are still in a state of flux. Let us, therefore, Just consider 
the initial transverse and longitudinal stability of the vessel after which we 
can investigate the effect upon this stability of the icebreaking process. 

The height of the center of buoyancy above the baseline can be easily, 
and quite accurately, estimated for the icebreaking type vessel by the well 
known expression 



f /£J 7 ~ 




( 128 ) 
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where all coefficients are, as throughout the paper, based on length between 
perpendiculars. Earlier in the paper we presented methods and graphical aids 
for the determination of first approximations of vessel displacement, length, 
maximum load waterline beam, draft, prismatic coefficient and block coefficient. 
In order to determine a reasonable value of waterplane coefficient to use in 
our above expression wo make use of figure 7 <> which shows a variation of 

4 

waterplane coefficient as a function of prismatic coefficient. Again, the 
values plotted are very typical of current and past icebreaker design practice. 
The transverse metecentric radius or height of metecenter above center of 



f buoyancy may be given as 



8» t - Z- - 
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C e t-d/T 1 
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// 



(129) 



where Q is the transverse inertia coefficient. An investigation of icebreaker 
type vessels indicate that the transverse inertia coefficient can be expressed 
as a function of load waterplane coefficient. The results are shown in figure 



7/ • 

We now have all the information necessary to estimate the location of 

ft 

transverse metecenter from (128) and (129). Therefore 

t - /C >3 /• 

In a similar manner we can examine longitudinal stability. Since the 
vertical distance from the baseline to the center of buoyancy in the longi- 
tudinal direction is equal to that in the transverse direction we have 

■ . / d /• 

where KB is given by expression (128). The longitudinal meti?centric radius 
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or height of met^center above center of buoyancy is given as 








( 130 ) 



where is the longitudinal inertia coefficient. An Investigation 

similar to that conducted for the transverse condition reveals that the longi- 
tudinal inertia coefficient may also be expressed as a function of load water- 
plane coefficient. The results for the icebreaking type vessel are 6hown in 
figure . Hence in a manner similar to that used in the case of transverse 
stability we may with reasonable accuracy determine the location of longi- 
tudinal metscenter. For this longitudinal case there is one qualification 
we can make. In ioebreaker type vessels the vertical height of center of 
gravity varies from about 0.95 to 1.2 times the full load draft, therefore, 
we may assume that the center of gravity is at the lead waterline and say that 



from the principal dimension and form characterisitcs of the vessel in the 
full load or design draft condition, it would be desirable now to be able to 
predict the effect on stability of typical operating conditions. When 
operating in an ioefield the ice conditions are often such as to require the 
vessel to overcome the ice by ramming. When an icebreaking vessel rams solid 
ice of considerable thickness the kinetic energy of the vessel in conjunction 
with the developed propeller thrust act to force the vessel up on to the ice 
due to the inclined stern. This act of climbing on to the ice results in a 



sy?_ - t /3 /■ - // 

Having estimated the position of longitudinal and transverse metecenter 



reduction of vessel draft and a change -in vessel twau 

t£)m 
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If wo consider the transverse stability first we can see that this change 
of draft may materially reduce the vessels stability due to a reduction in height 
of mettfcenter* In addition* since the vessels stem ie sitting on the ice shelf 
exerting a vertical force F, stability is still further reduced due to the 
virtual rise in KG caused by the partial vessel support by the ice sheet. This 
later consideration is similar to that encountered in drydocking a vessel with 
trimo The effect of change in vessel trim can normally be neglected if the trim 
is not excessive (more then about 1 % of the length). However* even in those 
instances where the trim does exceed the above limits due to the vessel riding 
up onto the ice, the vessel will be trimmed by the stern and due to vessel shape 
the effect will be to increase KM. Consequently, we can normally ignore the 
effect of triip and consider as the more severe case only the effect of change 
in vessel draft and grounding of the stem. As regards the effect of this vessel 
motion on longitudinal stability we maj state that the distance from the keel or 
baseline to the longitudinal metecenter is normally so large that any effect due 
to a change in vessel draft or vessel trim will be quite small and insignificant. 
Hence, we can safely ignore the longitudinal consideration so long as we have 
insured that the intact longitudinal metiscentric rediua is adequate. 

In considering the effect of change in vessel draft on the transverse 
stability it is useful to consider the concept of vertical prismatic coefficient 
for the representation of vertical distribution of displacement. This is giVen 



Since both block coefficient and waterplane coefficient vary by the same factor 
with changing drafts for the normal type vessel, we may justifiably assume that 



as 



0 V “ Cb / Cv 



and has been indirectly employed in our expression 




( 128 ) 
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vertical prismatic coefficient remains constant with draft. In order to deter- 
mine the ordinary range of values of vertical prismatic coefficient for ice- 
breaker type vessels an investigation was made of past and present vessels. 

The results of this investigation are shown in figure 73 wherein vertical 
prismatic coefficients are plotted as a function of full load prismatic 
coefficients and vessel length between perpendiculars. Consequently, we see 
that by using expression (128) and the vertical prismatic concept which is con- 
stant with draft we may determine the vertical center of buoyancy at any draft. 

In section F-2 we showed that the change in vessel draft due to riding up 
on to the ice could be approximated as 



where 7 ^ is the waterplane area at full load condition. Here the assumption 
was made that change in vessel draft would be small (less than 10$) hence the 
waterplane properties would not be materially effected. This is not really 
true since waterplane area will vary. The use of (131), however, will not 
lead to too serious an error in but if greater accuracy is required we 
can improve our approximation. 

There have been a considerable number of attempts to investigate the 
variation of form characteristics with varying draft but the most consistent 
and simple to apply seem to be those published by Dyer in his The Dr aft - 
Displacement Curves and their Derivitives wherein he formulates 
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(131) 
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(132) 



and 
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(133) 
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where the subscript one refers to any intermediate draft as distinguished 
from the full load or design draft. Consequently, in order to improve our 
estimate of we can use the expression as given in (131) as our first 

estimate to determine 

//, r /✓ - A J ( 134) 

Based on this first estimate draft X, , we can determine the final water- 
plane area, /7^ , for this draft by use of (133), Then we can correct (131) 
to get a more aocurate estimate of change in draft as 

A 7 ° ^ 

, //- *j, i 135 > 

This process can be repeated until the desired accuracy is achieved, however, 
for most instances in the preliminary design study it is sufficient to use 
either (131) or (135) to estimate the change in vessel draft. 

We must next consider the effect of this change in vessel draft upon 

C /j w 

0'#* = 



<-« A/ 



(129) 

This is not a6 easily done since we cannot predict the variation in vessel 
beam with draft. We can, however, circumvent this difficulty. Let us con- 
sider the expression for , i.e. 

/c - /c o' y- /J r-y/ f 



- (' a/J 



Ci & 



If we divide through by vessel draft this becomes 



a/ 
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3ince the flare at aide for a normal icebreaker type veaeel is carried well 
above the normal full load or design draft, vessel beam will decrease with 
decreased draft. Consequently, the beam to draft ratio will remain approximately 
constant. We may, therefore, justifiably assume that the beam to draft ratio 
is constant with draft for the range of drafts in which we are interested. 

Hence in (1*6) above, since C v is also constant with draft, transverse inertia 
coefficient and block coefficient ere the only variables we must consider as 
functions of changing draft. From expression (1JJ) we may dfctermine C w for 
our new draft and from figure 71 the transverse inertia coefficient corr- 

eaponding to this new waterplane coefficient. Expression (155 ) will ®lao 
allow us to estimate the block coefficient at our new draft. Consequently, 
we have 



W* 



I 

Cy + I 



t 




( > 3 ?) 



where B/H end C v are constant and evaluated at the full load draft and 

items with the subscript t.o refer to evaluation at the new draft f-l* : H ~ . 

Hence, we can estimate the height of transverse metfcenter at our new draft, 

A 

jtM*. » and reasonably predict the effect of draft redeuction upon metacentaic 
fadius. 

When the vessel aits with its bow on the ice shelf, part of the vessel 
weight is borne by the ice sheet and pert of it is wpter borne. 
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Thet portion of vessel weight (equal to the vertical bow reaction force P) 
supported by the ice sheet effects stability in a manner similar to the 
removal •£ a weight of magnitude F from the point of contact between the 
vessel and the ice sheet. Because of this, there will bee e rise in the 
vertical center of gravity. This virtual height of center of gravity may 
be expressed as 



K Cr„ 



Y(/ • C( r - F 
W - F 



where it is assumed that the average height of the point of contact of 
the stem end ice sheet ia H2 and KB refers to the intact height of center 
of gravity at full load of deaign draft. Therefore, the virtual metecentric 
height of the vessel aa it sits on the ice sheet ia 

Q V) v - - K (r v- 

By the above described methods we have been able to predict the effect 
of the icebreaker ramming process upon stability during the early stages 
of deaign. We have considered thie ramming and climbing up onto the ice 

t 

as the most serious operational evolution which could materially effect 
vessel intact stability. HoweTer, we must realize that our approach was 
idealized 
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ill* igMilil I )i i riJiBiV. In actual vessel operation 
in arctic or antarctic waters there will often be the perplexing and serious 
problem of topside iceing. The. presence of such ice on the upper reaches of 
the ship can be extremely serious and especially when breaking ice by ramming. 
The possible reduction in transverse metecentric radius when the vessel rides 
up on to the ice in conjunction with the often serious rise in vessel center 
of gravity due to top side ice can act to produce a dangerous reduction in 
vessel transverse stability. It is impossible to predict the actual quantita- 
tive effect of this topside ice upon stability since the amount formed and 

% 

its location above the main deck is so definitely a function of unpredictable 
weather and sea conditions as well a6 the amount of spray which is being taken 
by the vessel. Since this iced top ca n have such serious effect 

upon vessel intact stability it must be considered by the designer of an ice- 
breaking type vessel, and he must insure that his vessel has sufficient intact 
stability so as to be able to withstand the substantial use in vessel center 
of gravity caused by topside ice. Invariably, the prudent ship operator will 
attempt to keep this ice formation to a minimum, however, it may not be assumed 
since weather and sea conditions may often make this topside ice removal too 
hazardous or impossible to accomplish. Consequently, it behooves the ship 

designer to assume the probable presence of this topside ice and provide 
% 

adequate intact stability to withstand it. 

One other item which will effect ship stability and which should be con- 
sidered early in the design stage is the question of vessel depth or freeboard. 
The selection of a proper freeboard for the vessel will greatly contribute to 
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itB performance both, in an ice field and in the open sea. 

Since icebreakers are normally such short fat ships with near round 
bottoms, they are invariably bad performers in a heavy sea. Their reputa- 
tion for heavy rolling is well known and in any type of a head sea they 
pitch quite violently. Many types of ships which pitch badly are provided 
with a high freeboard at the bow plus a steeply rising sheer and strongly 
flared sections so as to help keep solid green water off the deck in heavy 
weather. Icebreakers could well use such a configuration of bow, especially 
so when in the higher or polar latitudes where air temperatures are low and 
spray coming over the bow causes serious ice formation on the topside struc- 
ture. Several disadvantages to such a bow configuration have, however, been 
raised in the past. The steeply rising sheer makes it difficult for men to 
keep their footing when the main deck forward is iced over and, consequently, 
makes ice removal or handling of ground tackle difficult and dangerous in 
polar regions. It is possible, by use of a raised forecastle deck to provide 
necessary freeboard forward and yet eliminates all, or about all, sheer. The 
use of this large freeboard forward, whether in the form of a raised fore- 
castle deck or steep deck sheer, will, to some extent, hamper the conning 
officers view of the ice just forward of the bow. Some designers raise serious 
objection to such a restriction of vision, however, it is felt that any such 
restriction of vision will not be severe and that the bridge deck is normally 
hi$i enou^i above the water to provide the conning officer with as much 
freedom of forward vision as he need have or even desires to have. The very 
slight reduction in visible area forward caused by a raised forecastle deck 
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or high freeboard is certainly negligible when compared to the seriousness 
of heavy topside iceing forward. 

As mentioned previously, icebreakers roll fantastically. In heavy seas 
a 30° roll is commonplace, a 40° roll not unusual and 50° or more not unheard 
of. In such situations, the danger to intact transverse stability due to deck 
edge immersion must be realized. It would be beneficial, therefore, if as 
generous a freeboard as is possible be provided over the midships half length. 
In this regard it might be beneficial or desirable to combine the need for 
generous freeboard forward as well as amidships and provide for a completely 
enclosed damage control or main deck to extend aft to about the after quarter 
point. This would provide adequately generous freeboard and, in addition, 
since the covered main or damage control deck is that for which floodable 
length is computed, it would be feasible to provide for continuous fore and 
aft inside passage ways so that personnel need not traverse weather decks 
during heavy weather. In the selection of vessel depth amidships, figure ? Y 
can serve as a useful guide. It shows the representative variation of depth 
amidships as a function of length between perpendiculars for varying ratios 
of length between perpendiculars to draft. The variation shown is typical of 
present day icebreaker practice. 

The choice of freeboard at the stem is again influenced by stability 
considerations. When a vessel breaks ice by ramming and her bow climbs up 
on to the ice there is a tendency for the stern to depress. If the deck 
becomes awash due to this change in vessel then the transverse mete- 

centric height will be materially reduced. This will aggravate the loss of 
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metecentric height already incurred due to change in vessel draft when riding 
up on to the ice and intaot transverse stability could be endangered. It is 
desirable, therefore, to provide sufficient freeboard aft to keep the deck from 
becoming awash when the vessel experiences a change in defined by the 

angular change of Section F-2, i.e, 

Cy 

In order to estimate the value of ^ , the distance from the intersection of 

stem and load waterline and the longitudinal center of flotation, we can say 



that 



3 - 



where yj • distance from amidships to longitudinal center of floatation with 
aft positive and forward negative. Figure 7 ^ shows the variation of sr? 
for ioebreaking type vessels as a function of waterplane coefficient. For 
values of < 0.715 we may assume that /h ~ o with negligible introduced 

error . 



In summation, therefore, the effects of freeboard are increased vessel 

I 

weight, rise in vessel center of gravity and increased lateral plane, area. 
Weight is no problem in an icebreaker, in fact, it is in a way desirable so 
aB to increase ice breaking effectiveness. As regards vertical center of 
gravity, it is true it will rise but the transverse metecentric height of an 
icebreaker is normally very large and some slight reduction due to an inoreased 
f<7 Q- oan be tolerated and probably absorbed without noticeable ill effect. 
Since winds in the polar regions are very strong, the increased lateral plane 



area will cause increased wind resistance to motion. How serious an increase 
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is difficult to predict but it is doubtful if the magnitude will be such 
as to cause serious effects upon vessel acceleration ability. This increased 
lateral plane area could also cause a slight effect upon vessel directional 
stability or course keeping ability due to wind action. However, since ice- 
breakers are usually such deep drafted vessels, the noticeable effects should 
be negligible. 
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Part 3 



Illustrative Example of Methods 
A. Preliminary Design Estimates 

Consider an icebreaking vessel to be built for polar service whose 
primary area of operation is to be the Antarctic with possible alternate 
summer deployment in the Arctio. The vessel is to be a first line breaker 
capable of early independent penetration of the pack ice under normal ice 
conditions during the month of October in the Kanian Bay, Little America area. 

Based on a study of the probable areas of vessel operation in the Ant- 
arctic, let us assume that it is necessary to restrict vessel design full 
load draft to 30’ -0*. In addition, let us assume that this area study has 
revealed that the probable maximum ice thickness to be encountered during 
late winter or early spring is 15’ -0” of one season ice with some possibility 
of older ice. 

The items known at the beginning of , our study are, 

Maximum draft - 30’ -0” 

Maximum ice thickness - 15’ -0” 

Month of First Pene- 
tration s October 

From table 1 we see that the mean air temperature for the Kanian Bay, 
Little America area for October is 

^ * 



- 30° F. s - 30.9° C 



- 3 ° f - 



The ice surface temperature, assuming no snow cover may be approximated as 

^ = 0.86 ( 0^4 0.1?) - 0.4 °C. 112) 

• -33.7° C say - 34° C. 

Since there is a possibility of old ice as well as younger one season ice 

in the area we shall consider perennial sea ice as the most severe case. 

Based on our limiting draft of 30* -0" we Tiave from figure vj" 

propeller diameter _ q 
D raft 



Therefore 

« 

Propeller Diameter D a 18* -5" 

In addition, from figure V 7 

IBP s L - 398* -0” (First Estimate) 

Based on this estimate of vessel length, the extreme beam at the design water 
line is from f igure " V& 

B* s 79* -3** (First Estimate) 

From figure so W e note that characteristic propeller loading factors 
■vf / D w are not given for vessel lengths in excess of 300* -O'*. In order 
to determine a value for use with our problem we must try to interpret the 
indicated trends. As can be seen, the ratio ^ DV tends to rise rapidly 
for vessels beyond 275* - O'* in length. In section F-2, however, we indicated 
that from an economic viewpoint the value of this ratio approaches a limit 
of about 850. In our problem, therefore, we may assume the limiting value, 

850 lb/ft 2 per shaft 
289000 lbs. 

578000 lbs. a 258 tens 






/cT S 



per shaft 
TOTAL TROOST 
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Assuming the conventional practice of a built up three bladed propeller 
with thick strength seotions, we have 

^ 7~r(l OOS7 

- 196 tons 

Knowing the ideal propeller thrust (Troost) per blade and propeller diameter, 
we can, from figure f determine an approximation of shaft horsepower 
required as 

SHP/shaft ■ 15800 

SHP Total = 27600 

and from figure the propeller RPM as 

RPM = 150 

If the thickness of ice we wish to overcome is 180" and surface temperature 
is -34° C. we can estimate from figure 36> for perennial sea ice that vertical 
bow reaction force F required as 

F = 4900000 lbs. s 2185 tons 

Let us assume that in normal breaking by ramming the accelerating distance 
is equal to one skips length, therefore 
_5 ' - L s 398 » - 0" 

In addition, based on our vessel length, a first approximation of vessel 
free route speed is 

y = 0.978 / 398 = 19.5 knots 

We can now apply the expression for icebreaking force and vessel impact speed 
so as to obtain a reasonable first estimate of vessel displacement. Dealing 
first with icebreaker force we have the modified expression (90), i.e. 

/S~' 2 - — 

j/' / r ' L /y - z* z c r 7^ z 7 f/ 



UJ 



- 7 * 7 - 



Substituting our known values this can be simplified to 

200000 

10 ■ — 

<r 

Assuming various values cf we can obtain a plot of vessel displacement 
versus vessel impact speed as shown in figure ?& • Now consider vessel impact 
speed due to acceleration both in open water and in an ice clogged channel. 
Considering first acceleration in open water we have 

_ t a . im / 

XT S 3.1 X 72 

therefore, substituting known values 

.*• ' s 0.0121 *r 

Assuming various values of vessel impact speed at we can obtain a variation 
of Y as a function of ■ / v . From figure 7/ we obtain values of F(x') 
corresponding to each value of X^. These values of F(X^) used in our expression 
for open water 

2.86 V 2 _ 307 

Fix** - ) x> 

will give a variation of vessel displacement as a function of impact speed 
as shown on figure 7<&> . Since this curve and the one above give a simul- 
taneous solution to our problem, we have from the intersection of the two 
curves a first estimate of displacement of 
W s 10200 tons 

This, however, is an optimistic estimate since the vessel will seldom 
accelerate in open water. Instead, the working channel will soon become 
full of ice pieces and the vessel will, under normal circumstances be 
accelerating in an ice filled channel. For acceleration In an ice channel, 
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Sample Calculation 



uj - 



ip y -*-*y r . xr * 



Uj - 



—— i /( z / sr ) ( jo ) ~(a e ;#±' X i0 J 



2 < p > & o c>o 

y\T 



open water 

, , / O 6 1 t 

at r 3 .fx r m 

yyr 



ice channel 



^ - A /^> *71 



> . ** 2 / 



x p 



(j. /)(/?<.) 

x o.o j z/ 'iT~ 

i / " 



» 6. i f X P 






i / 



. ( z . tfiX ' f- s )j _ 

jo ? 

3 c * ■) 




F(X )( ' ><*J 



t, 3 o o 



* v ; 



ov ; 
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2 1 


100000 


0.0242 


0 .0003 


1024000 


0.113 


0.004 


1573000 


4 


50000 


0.0484 


0.0013 


236500 


0.226 


0.025 


252000 


6 


33350 


0.0726 


0.0025 


123000 


0.339 


0.06 


105000 


8 


25000 


0.0969 


0.0045 


68300 


0.452 


0.12 


52500 


10 


20000 


0 . 121 "- 


0.0075 


41000 


0.566 ' 

* 


0.20 


31500 


12 


16700 


0.145 


0.0105 


29300 


0.679 


0.32 


19700 


14 


14300 


0.1693 


0.0142 


21600 


0.792 


0.515 


12230 


16 


12500 


0.1936 


0.0193 


15900 


0.906 


0.95 


6630 


18 


* 

11100 


0.218 


0.0245 


12530 








20 


10000 


0.242 


0.0310 


9930 


1 












- — — 1 











nr 


U) 


Jf ' 


/?(*•) 


L*J 


■X ' 


S ( a ') 


LXJ 






* 



i 

i 
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therefore, 



and 



IT « 1.12 X ' 71 0,524 
21.25 V 2 



= 



F(r') 



— * . *> */ S 



Substituting our known values these expressions become 
,Y ' s 0.0566 ^ 
s 6300 

FU) 

By a method identical to that used for acceleration in open water and with 
figure */ 3 , we can obtain a new curve of vessel displacement versus impact 
speed as shown in figure 76. Again, the Intersection of this curve and our 
first curve gives a simultaneous solution to our problem. Hence a more 
reasonable displacement estimate is 
W - 15200 tons 

In order to now consider the compatability of our estimated dimensions 
with a reasonable block coefficient, we proceed next to an estimate of vessel 
form parameters. Figure ^7 shows a plot of prismatic coefficient as a 
function of speed length ratio and displacement length ratio. For our case 
Vnr = 0-978 

X-/= 241 

Since ours will be a diesel electric vessel, we see that figure & 7 does 
not extend to as low a displacement length ratio as ours. We can, however, 
use the curve for the lowest value listed, i.e. («-ou) « 300. Thus 

Cp s 0.696 

Based on this value of Cp and our speed length ratio we can determine a 
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representative value of block coefficient from figure £8 , i.e 



C B = 0.545 



and 




Now 



P 35W 

°b = Jm 



Assuming that for our vessel we desire 



L/ b : 4.5 



which is reasonable for a large polar icebreaker, we have 



L z 12.57 ■/ c Jh - 384 ft 



hence our first estimate of length was a bit high and for our second estimate 
we shall assume 

LBP = 384 ft 

Therefore our second estimate of vessel maximum beam becomes 
B r 0.2225 L = 85.4 ft 

Now we shall investigate the effect of these modifications in vessel length 
and beam on our other selected vessel characteristics. 

Since our vessel draft is still fixed at 30* - 0" our propeller diameter 
cannot be modified. Our altered vessel length would normally modify the ratio 
' r,, /£> w , however, since our vessel is still considerably in excess of 

300' - 0" we are justified in continuing to assume that powering will be at 
or near our discussed economic maximum 




= 850 



I 



, 







I 
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Therefore, our previously determined values of actual propeller thrust, shaft 
horsepower and shaft HPM remain unchanged. Vessel vertical bow reaction force, 
F, necessary to overcome ice 15* - 0 H thick is not a function of vessel length 
or beam, hence it too remains unchanged. 

Again we assume 

$ rz L = 384 ft 



Therefore our estimated vessel free route speed becomes 

y z 0.978 y r 384 = 19.2 knots 

If we now reapply our expression for vessel displacement 



W s 



z 



/ 



V - 2 • z y /7 



and that for vessel impact speed in an ioe channel 

o . r *- V 



-tr 



- 1.12 x K 



w = 



./ ■ e r / 



- I * v /. 



• o y ) 



we obtain a new graphical solution for vessel displacement. Doing this, however, 
we note only a minor increase in our displacement estimate due solely to the 
modification of vessel free route speed with length since all other items remain 
constant, i.e. 



W r 15300 tons 

Since the effect of this minor change in displacement will produce a negligible 
effect upon our estimated prismatic and block coefficients we now have a reason- 
able and compatible preliminary estimate of vessel proportions. 

v 

To further determine the adequacy of these proportions, we should now 
investigate vessel stability. Based on our value of Cp we can from figure 7 o 
determine a reasonable value of load waterplane coefficient, i.e. 



3' 3 



C w = 0.7955 



From expression (128) we can determine height of center of buoyanco as 






Cm + <- 



(H) s 17.8 ft 



The height of transverse metercenter above oenter of buoyancy is given as 
(129) 






a 



C rt /r 



where transverse inertia coefficient C t ■ determined from figure 7/ as 

j . s 0.0525 



therefore 



fdrrO t - 23.4 ft 

and 



l ,r! ^ - 16 ■ 0 1 -V] t — 41.2 ft 

Since we have shown that for icebreakers the vertical height of center of 
gravity may be taken at the design water line, we have 

A(r - 30’ -0" 



Therefore metecentric height for our vessel is 

7 ,y ? c ~ f- y >7 tr - ,4 ( r B 11.2 ft 

Considering now the longitudinal condition 



b>v L = C - c 

c 9 v 

where C, is obtained from figure 77 as 

C- s 0.0478 



therefore 



( 130 ) 



b>r). 



Z 431 ft 



] 



I 

I 
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Let us now consider the effects of one of the most severe operating 
conditions on this stability, i.e. the vessel sitting with its bow up on the 
ice shelf. Since resulting vessel is small we can neglect its effect 

and consider only change in vessel draft. A first approximation of this 
change in draft was given as (131) 

4-0 

where 

/L r A /J C u, = 26050 ft 2 

hence 



This, however, 
mation 

where 

No?j 

hence 

and 

Therefore 



4 A * 2.94 ft 

assumed constant watorplane area, hence as a closer approxi- 



a4, = 



4 ^ / 4 ^ 






-- ! C, 
>/ 



(133) 






= 27.06 ft 



<03 
(- o 



= 0.685 



4. 



24900 ft 2 



A A , = 3.01 ft 









V - 



26.99 ft 






' 
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Since vortical prismatic coefficient, C v , is constant with draft, the 
new height of center of buoyancy is 

(Ml) ' /6 ° u 

we have shown that * 



KM 



tv 



C J± ( 6 \ 

C B . ^ 



Cv 1 > 



where C v and B/H are constant with draft and 

/-<V 



c *. 



therefore 



and frtom figure 71 



Hence 



, M \— y 
CV, ' l / , 



C 8i z O. SiO 
C». r o • O 



c[ t = O ■ © v <f?/ 

_ j 3 <?. 5 ~ 

= 36.2 # 



and 



( /*>>> 






** 

The vertical location of virtual center of gravity is given as 

^ W. AT, a A*//! 

A- Cr^ ~ — * 

u/ 

hence 

GM ti = jL 

V 

This certainly is an adequate QM^ for the inclined condition and 9eema 



JO, <?■ 5' // 

- S . ?3T H 



aufficient to care for any rise in center of gravity due to possible 
topside ice. Our preliminary vessel proportions, therefore, seem adequate 



end reasonable. 
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At this point there are several other items we a an investigate based 
on oar preliminary proportions. These are: size of vertical forefoot and 

location, vessel retraction ability and vessel icebreaking potential without 
ramming. 



Considering first the forefoot we have that the vertical forefoot will 
not strike the ice sheet if 

A ^ c 

where c is given as 



C = / 






U) 1 



/• <0.0 *~ 



*' 7T 



/.* </r~ 



// 



A 



o o oo s' vy 




Substituting we have 

C s 0.1583 

Bence, in order for the vertical forefoot not to strike 
Jl + 0.1683 




or sinoe vessel draft is 30*“0* F the height of vertical forefoot 



JiW ^ 4.76 ft 



- 3J7 ■ 



0 

Based on our value of K z 26 , therefore, the vertioal forefoot should 
be 4* -9* high and located about 51' -8" aft of the forward perpendicular « 
Considering now the ability of the vessel to retract if once beBet 
we have for acceleration in an ice channel 

-z-IC-L • uj''*'"'' 

t; > 0.0000234 - ■ - 



in order that the vessel be able to retraot under normal conditions. For 

• s» 

our vessel 



O. fru 



= G>l) 



o. fSl, 



r 153 tons 



hence since 



J ooa-cy Zi / 



to jT 
// 



/ , ) (/->'iooJ 

rr [o oooo 23 V/ r* 7?, J 



Z^w 



0.0000234 



^ V 
// 



our vessel should be capable of extracting itself under normal breaking con- 
ditions. For the more severe case where the vessel has been moored into the 
ioe for a period of time and the stem has adhered to the ice we must satisfy 
the condition 

60 '" / r ' 

0.000229 — — 

/v 

. 153 tons <1 756 tons 

In this instance, therefore, the vessel will not be able to extract itself 
and heeling tanks or some other auxiliary method must be used to free the 
vessel. 



Turning now to the vessels icebreaking ability without ramming, l.e. 
continuous breeking, we see from section F-7 that the forces involved may be 

•4 

expressed as 



/- = 2.52 



3 j c? ~ 



Based, on our value of bollard pull, the Vortical bow reaction force con- 
tinuously available without ramming is 

— 2 • 52 x 196 * 494 tons 

. 1107000 lbs 

If we assume ice conditions as in the statement of our original problem we 
find from figure Sg that the maximum thickness of ice which our vessel can 
break continuously is 

£ " - 86" = 7.17 ft 

At this point we have arrived at pretty reasonable estimates of vessel pro- 
portions and parameters. V/e may now proceed to a determination of vessel 
lines. Reiterating those items which we have thus far determined we have 



LBF = 


384' -0" 


W 


m 


15300 tons 


H r 


30* -0" 




a 


0.696 


B - 


85' -5" 


L & 


— 


0.545 


D » 


18* -5" 




= 


0.784 


V 9 


19.2 knot6 




& 


0.7955 


SHP « 


27600 




- 


30° 


EPLS e 


150 


i 


s 


45° 


Blades = 


3 




9* 


0.1583 




I 



/ 



B. Evaluation of an Existing Vessel 

The method outlined in this paper may be used to evaluate the effective- 
ness of an existing vessel as well as provide a basis for preliminary design 
estimates. In this regard it gives the designer a means whereby he can compare 
several existing vessels to determine their relative icebreaking effectiveness. 

Contrasted to the preliminary design consideration where only relatively 
few items concerning the proposed vessel are known, we now have the condition 
where we know all vessel dimensions and parameters 0 As an example of our 
approach in this instance let us consider the potential icebreaking abilities 
of the Glacier as representative of a typical present day polar icebreaker. 

The characteristics of the vessel axe as follows: 



LBP 


* 


290’ -O'* 




s 


0.537 


B 


r 


74* -O'* 




3 


0.800 


H 


m 


25* -9” 


^ P 


m 


0.698 


V 


• 


16 knots 




= 


149.6' 


SHF 


z 


21000 


1 7*? t 




9.95’ 


D 


s 


17* -6" 




= 


35.6* 


W 


= 


7750 tons 


/c d 


a 


15.4' 


KG 


S 


25.65' 


X 


- 


30° 




X 


20.2’ 


J 


- 


cr> 

o 




& 


283 




s 


0 *466 



Since this is a twin screw vessel, the shaft horsepower per shaft is 



one half the total given above. The vessel mas designed for three bladed 
propellers, hence from figure ? based on SHP/ shaft and propeller dia- 
meter 



SHP/shaft s 10500 



D = 17»-6" 



T TR00ST o 235000 lbs/ shaft 
T 

TR00ST = 470000 lbs total 



= 210 tons 



Correcting for blade thickness and increased hub diameter due to use of built 
up propellers 



T ACTUAL = Ki k 2 troost 

• (0.86) (0.885) (210) - 160 tons 



£ _ * \ ' 

w ^ /o 



1.031 



L 



/ 




) 



{){ o rJiXc . y) 

( o S 1 ( ^ yoj 



0.257 



H - 





' ] w .ti) 

1 r/f.TvJ 



4.051 



For acceleration in open water, based on an accelerating distance of one 



ship length 






I, 
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* 4.0&1 



For acceleration in open water, based on an- a cc e l e r a t ing^ distance of one 
skip length 



W = 



2.86 V‘ 



FIX 1 ) 77 0,24 



FIX 1 ) . . o.028. 

(7750)(160) 0,24 



0.621 



. . X A s 0.229 

V = 3.1 x 

s (3.1 ) (0.229 ) ( 160) 
s 16.6 ft/sec 

l 

If acceleration is in an ice channel 



0.621 



W = 



21.25 V s 



F(X X ) 77 0,047 

.*. FlX 1 ) » (21.25) ( 16) 2 - 0.552 

( 7750) (160) 0,047 

.*. X 1 a 0.806 

0.524 



V a 1.12 X 1 0 
T s (1.12)(0.806)(160) 

= 12.9 ft/sec 



0.524 



o 0 

r 30 J = 46 



. J a 1.09 



v - 0.548 

C ”* 



I 



- 3 <: z - 



Assuming for the moment that the vertical forefoot does not strike the ice 
shelf, the vertical bow reaction, F, which may be developed is 
Open water acceleration: 



F 



0.91 x r; 



o- r.y/ V 



J 



*4 /V 



0 ■<//)(, 3>y\_/(>o) 









(s?y°) ( / » ■ o. y-vt-J 

r */■ « ?>-> 



- 1808 . 5 tons 

a 4040000 lbs 



Ice channel acceleration: 
F - 1443.5 tons 
s 3230000 lbs 



There is, however, a possibility that the vertical forefoot will strike the 
ice sheet and hence restrict the vertical bow reaction force which may be 
developed. In order for the vertioal forefoot not to strike 



- ■_ 

therefore assuming acceleration in an ice channel as most characteristic 



a • > i 

- / - - / 



*■ / ' z ~ x > • '7~ 

— v- ru l ' 

/ -u v 



y — 

! ,, 



/ 



o Ci r* 



jt ' . T-/-***' 






L 



7 ?J O 






rvoj C / *oJ / ° V ' r 
4 r. 



° O o o o \ 



3 ■ irac.)'" r 



C S 0.083 
-*L a 0.46 



. . . 4. ,> 



Hence the vertical forefoot will strike and the vertical bow reaction 
force, F, which may be developed is limited. Therefore, the vertical bow 
reaction force we can develop is 

s " 

^ ^ . o // 

, — A-n ° ~> O > y*z > r-~ 

'1 / 1 

- v/ ; 5 ^ 

, 7 , 

I 4. }~si Jl*/- I I °° U ' !TO’0J i/toj i^-v, ^ o 



s*s 



A 






This represents a considerable loss of potential icebreaking ability. 
However* if we assume, as we did in Part A, that the vessel is to be used 
in the Antarctic, Little Amerioa area during the month of October we have 

r 3 ~ 34^ C • 

o 

Hence the amount of ice which can be broken by this vessel is, from figure 

r~ - 101 * 

s 8.43 feet 

Of course, this assumes i°/io coverage or a continuous unbroken icefield. 
Let us now investigate the effect upon stability of the vessel sitting on 



the ice sheet. 



- ' 



As a first approximation of change in vessel draft 



AS = 



_ 35F 



A 



where 



A u s LBCw » 17200 ft 



« » ” 1*37 ft 

This assumes constant waterplane area. For a closer approximation 



v/here 



hence 



and 



Therefore 



Now 



&J. * 



o - 



A 

x/* 



A, / A 



/A I A 

i ^ ' 



/V - /-/ - A J - 

- _ A 



c. 



o . b ? / 



A W1 . 16750 ft* 
= X - 39 ft 






M 



V - Aj = 24.36 ft 



G, , " 



where C v and b/H are constant with draft. Also 



A 

A 



A 



/- c. 

'A ) A, 
A ' 
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Therefore 



o< 3 



Cy 



- 



end from figure 71 



Hence 



* 

~*r 

KM t 



o . 77 if 



o . o S 03 



- /. 3 96 



3 4 . o j-i; 



The vertical location of virtual canter of gravity 

K <r-v- - - ZS. S ji 



W - /=• 



Hence 



- < 9-2 U 



CyMt t = ILMt, - ^ 

This is certsinly adequate even allowing for topside icing. It is interesting 
to note that this negligible loss of stability is due solely to the serious 
restriction in bow reaction force F and vessel climbing potential due to the 
vertical forefoot striking the ice sheet. 
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Section H 

Table of Characteristics of Past and Present 
Icebreaking Vessels 



K 



- 3Z?~ 



* Mote - all values given are for the design condition as far as possible. 

Under item $27 SHP refers to diesel electric installation, IHP to 
steam installation 



1. 


Nam© 


ABEGWEIT 


ALEXANDER 

HENRY 


ALGONQUIN 


— 

ALTE 


APU 


CACTUS 


D'lHERVIL 


2. 


Tear Built 


1947 


1960 


1934 


1926 


1899 


1942 . 


1953 


3. 


nationality 


CANADA 


CANADA 


U.S. 


SWEDEN 


FINLAND 


U.S. 


CANADA 


1 ' ' 
















4. 


IBP 


348-0 




150-0 


194-7 


141-9 


170-0 


285-0 


5. 


LOA 


372-6 


210-0 


165-0 


204-1 


144-0 


180-0 


310-0 


6. 


LWL 


355-0 












302-0 


7. 


B (Extreme) 


62-0 


43-6 


36-0 


55-9 


35-6 


37-0 


66-6 


8. 


B (LWL) 


60-0 




36-0 


53-2 


34-0 


35-0 


64-6 


9. 


Depth 


24-9 


21-0 


21-0 


28-9 

t 




17-5 


40-0 


10. 


Draft - H 


19-0 


16-0 


12-0 


V 

19-8 


18-0 


12-0 


27-6 


'll. 


Displacement - W 


7600 


2400 


1000 


2464 


800 


940 


8700 


12. 


Design Speed - V 


16.5 


13.5 


12.8 


15.5 


12 


13 


15 


13. 


Block Coeff - C b 


0.513 


0.650 


0.538 


0.418 




0.461 


0.580 


14. 


Prismatic Coeff - Cp 


0.571 




0.632 


0.571 




0.565 


0.690 


16. 


sfcr/t>v 

Coeff - C x 


0.90 




0.852 


0.733 




0.817 


0.840 


16. 


Waterplane Coeff -C* 


0.791 






0.662 




0.728 




17. 


ZB 












7.7 




18. 


za 












14.71 




19. 


“t 












17.57 




20. 


ZM t 












215 




21. 


0M t 












2.86 




22. 


“t 












9.87 




23. 


LCF from 












1<>85 

aft 





- / & s : 



24. 


LCO from <& 








- 




3.7 

aft 


25. 


type propulsion 


Diesel 

Electric 


Direct 

Diesel 


Steam 


Steam , 


Steam . 


Diesel 

Electric 


26 . 


Number of shafts- aft 


2 


2 


1 


1 


1* . 


1 




Fwd 


2 


0 


0 


1 


0 


0 


2 ?. 

f 


SHP or IHP (total) 


15400 


3550 


1500 


4000 


1350 


1000/ 

1200 


28 . 


Wetted Surfaoe 


* 










6563 


29. 


Propelled dia - D 


13-0 


8-6 






' 


8-6 


30.. 


Number of blades 






‘ 






5 


31. 


Propeller pitch 


13-6 


6-0 








7-0 


32. 


EPM at design speed 


128/ 

155 


300 


140 


100/ 

125 


84/ 

90 


208 


33. 


Mean width Ratio 




4 








0.333 


34. 


Blade thickness 
fraction 












0.0458 


35. 


Thrust deduction 
factor 












0.240 


3$. 


f 

Wake fraction 






4 




4 


0.205 


37. 


Propeller top 
clearance 












0.6 


38. 


Rise of forefoot <X 


23.5° 


26° 


33° 


24° 


20°- 


30° 


39. 


ft (approx) 














40. 


Flare Amidships 


go 


0° 


4° 


16° 




20® 


41. 


Rudder Area 


142.5 






• 


• 




62. 


Lateral Plane Area 


6750 












43. 


Propeller developed 
area ratio 


0.528 













Steam 

2 

0 

10800 



14-0 

4 

13-6 

145 



2.35 

30* 

57° 

60 

172.6 

8400 



I 



0.50 



- 



EDWARD 

CORNWALLIS 

1949 

CANADA 



ELBJORN 

1953 

DENMARK 



ERMAK 

1899 

RUSSIA 



ERNEST 
LA POINTE 

1940 

CANADA 



FEODOR 

LITKE 

1909 

RUSSIA 



GAIDAMAK 

1898 



GENERAL 
SAN MARTIN 

1954 

ARGENTINA 



GOTA 

LEJON 

1933 

SWEDEN 



GLACIER 
1955 
„ U.S. 






240-0 

259-0 

43-6 

20-6 

18-0 

3700 

13.5 

0.68 



157-6 

167-6 

39-5 

39-4 

16—2 

1400 



310-7 

320-0 

71-6 

70-1 

42-6 

24-0 

7875 

14 

0.524 

0.610 

0.859 

0.721 



172-0 

184-0 

36-0 

17-0 

15-0 

1675 

13 

0.632 



250-0 

265-0 

47-11 

47-6 

26-6 

16-9 

2570 

18.5 

0.452 



163-0 

198-7 

42-0 

41-6 

18-0 

12-2 

1030 

12 

0.438 



242-10 

277-10 

252-6 

62-4 

61-0 

32-3 

21-4 

4830 

16 



173-11 

183-9 

49-5 

46-11 

28-3 

20-0 

1900 

0.408 



290-0 

309-7 

74-0 

72-6 

38-0 

25-9 

7750 

16 

0.537 

0.698 

0.770 

0.800 

15.40 

25.65 

35.6 
283 
9.95 
20.2 

4.6 aft 



1.4 aft 



3 3o - 



Steam 

2 

0 

3500 

10-6 

* 

11-9 

130/ 

140 



25° 



1 

1 

2700 



160 



Steam 

3 

1 

9000 



105 



25° 



67. 5 C 



20 ° 



Steam 

2 

0 

2000 

10-0 

10-9 

140 



20 



0.446 



Steam 

2 

0 

7000 



125 



52 c 



Steam 

1 

0 

1300 



83/ 

94 



30 



Diesel 

Electric 

2 

0 

6500 



135 



Steam 

1 

1 

3800 



120 



23 



20' 1 



Diesel 

Electrio 

2 

0 

21000 

25920 

17-6 

3 

12-1 

175 

0.322 

0.049 



0.75 

30° 

46° 

o 

20 

173 

6670 

0.382 



f 
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HOLUEH 

DANSK 


HUDSON 


f 

HULL 614 j 


HULL 620 


INTO 


ISBJOEN 


ISBRTTAEEN 


ISBRYTAREN 

II 


JAAKAHHU 


1942 


1934 


j 






1923 


1895 


1914 


1926 


DENMARK 


U.S. 


CANADA j 


CANADA 


FINLAND 


DENMARK 


SWEDEN 


SWEDEN 


FINLAND 


210-4 


104-0 


f 


280-0 


254-1 


156-2 


134-0 


190-7 


246-0 


226-0 


110-6 


223-6 


! 

315-0 


273-10 


170-7 


139-9 


200-0 


263-0 » 








307-3 












56-6 


24-0 


48-0 


70-0 


63-6 


40-2 


35-6 


55-10 


63-1 


53-5 


24-0 




69-0 


61-2 


39-4 


34-0 


52-6 


60-2 


24-6 


12-7 


21-0 


41-0 




21-0 




28-8 


31-10 


17-10 


9-6 


16-0 


28-0 


20-4 


18-6 


18-10 


20-9 


21-10 


3000 


308 


2700 


8974 


4415 


1330 


720 


2350 


4850 


15 


12 


13.5 


15.5 


S 16.5 


12.5 


12 


14 


13.5 


0.533 


0.456 


0.650 


0.560 


0.482 


0.408 


0.410 


0.397 


0.525 


0.634 


0.618 












0.527 


0.666 


0.842 


0.738 












0.755 


0.789 


0.730 


7.95 














0.741 




12.40 


















14.12 


















106.95 


















1.72 


















6.17 


















2.53 

aft 
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Diesel 

Electric 


Diesel 

Electric 


Diesel 

Electric 


Diesel 

Electric 




t 

* 

j 

. Steam 


i 

Steam 


) 

1 

Steam 


2 


1 


2 


3 


2 


1 


1 


l 

1 


1 


0 


0 


0 


! 2 


0 

. 


0 


1 


6610 


800 


4250 


15000 


| 


j 2500 


1200 


2600 




7-6 


11-6 


13-6 


13-10 










3 
















6-3 


12-9 


13-0 




f 






123 


250 


145 


136/ 

170 

♦ 




j ' 90 


102 


130 

i 








2.5 

| 










87° 


30° 


26° 


| 30° 




20° 


31° 


20° 


46.5° 




j 

( 


■50° 










16° 


220 


0° j 


6° 






11° 


• 24° 






t 


204.1 








* 








'8820 














j 

0.519 


i 

j 

1 














% 











Stesm 

2 

1 

7840 



117 



26° 

20 ° 
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KAPITAN 

CLASS 

1954 

RUSSIA 


KARNU 

GLASS 

1959 

FINLAND 


KICKAPOO 

1922 

U.S. 


KIROV 

CLASS 

1940 

RUSSIA 


KRIS JAN IS 
VALDEMARS 

1925 

RUSSIA 


LABRADOR 

1954 

CANADA 


LEONID 

KRASIN 

1917 

RUSSIA 


LENIN 

1917 

RUSSIA 


LENIN 

1959 

RUSSIA 


265.0 


224-0 


151-3 


332-0 


185-0 


250-0 


297-0 


273-0 




272-10 


243-3 


157-0 


357-8 


196-6 


269-0 


323-3 


281-0 


440-0 


63-6 


55-7 


35-0 


73-2 


55-10 


63-6 


70-10 


63-10 


90-6 


63-0 


54-9 


35-0 


69-6 


54-0 


56-8 


70-6 






31-2 


28-10 


16-8 




28-6 


37-9 


41-0 


31-11 


52-11 


23-0 


19-0 


12-8 


23-10 


22-0 


29-1 


26-0 


19-0 


30-3 


5360 


3370 


750 


8330 


2800 


6490 


8730 


5074 


16000 


16 




11.5 




15.2 


16 


15 


16 


18 


0.469 


0.506 


0.392 


0.532 


0.446 


0.469 


0.562 






0.597 










0.622 








0.819 










0.754 








0.710 










0.730 
























6.24 


















% 



- 33 '/- 



Diesel 

Electric 

2 


J Diesel 
} Electric i 

i \ 

1 2 | 


i | 

■ Steam * 
1 


2 

10500 


2 

i ! 

j 7500 


0 

1 

850 


120 






22° 


1 

25° 


29° 


65° 






20° 




0° 



Diesel 

Ileotric 


Steam 


Diesel 

Eleotric 


Steam 


Steam 


Atomic 

Turbo-Electrl 


3 


1 


2 


3 


2 


3 


0 


1 


0 


0 


1 


0 


12000 


5200 


10000 


10000 


7500 


39200 




16-0 


15-0 






' 




i 

4 


3 












12—6 






♦ 




122 


105/ 

145 






205 




33° 


30° 


25° 


18° 




- • 




50° 










17° 


20° 


20° 


25° 








137 












6430 


• 
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MACKINAW 

1944 

U.S. 


MBOJLA 

SELIANINOV- 

ITCH 

1916 

RUSSIA 


MONTCALM 

1957 

CANADA 


MOSCOW 

1960 

RUSSIA 


I 

MURTAJA 

1890 

FINLAND 


nadeshni. 

1897 

RUSSIA 


NAUGATUCK 
r CLASS 

1939 
| U.S. 


N.B. 

MCLEAN 

1929 

CANADA 


NOBTHLAl 

1927 

U.S. 


280-0 


275-0 


201-10 


368-9 


137-5 


180-0 


105-0 


260-0 


200-0 


290-0 


292-0 


220-0 


400-7 


156-0 


192-0 


J 

110-0 


277-0 


216-7 






209-10 














74-4 


57-6 


48-0 


80-4 


36-1 


42-11 


26-5 


60-4 


39-0 


70-0 




48-0 




35-4 


42-0 


25-0 


59-4 


39-0 




28-1 


32-0 


21-0 




24-11 


25-5 


15-0 


31-0 


• 

24-9 




19-2 


19-3 


16-0 


34-5 


15-7 


17-9 


11-4 


19-6 


13-8 




5252 


4950 


2756 


15000 


825 


. 1710 


361 


5034 


1785 




16 


15.5 


13 


18.3 


12.5 


10/14 


12.5 


15 


11.7 




0.497 , 


0.574 


0.690 




0.382 


0.443 


0.417 


0.582 


0.611 


| 


0.612 


► 










0,572 




0.635 


! 

1 


0.812 












0.728 




0.963 


1 


0.728 












0.738 


! 






11.13 












7.13 


1 






16.52 


| 

i 


I 

i 








11.45 








34.68 


i 


! 








13.45 








338 












105 








18.16 












2.0 








23.55 












6.32 








1.0 




















aft 




















1.1 












1.23 








aft 


1 








1 


fwd 
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Diesel 

Electric 


Steam 


Steam 


Diesel 

Electric 


Steam 


Steam 


Diesel 

Eleotrio 


Steam 


Diesel 

Electric 


2 


2 


2 


3 


1 


1 


1 


2 


1 


1 


0 


0 


0 


0 


0 


0 


0 


0 


10000 


8000 


4000 


22000 


1200 


2475 


1000 


6500 


1 1025 



18975 

17-0 

3 



170 

0.327 

0.0559 

0.27 



11-6 



8-4 



106 



14-6 



11-7 




11-3 1 






5-0 


18-0 


136/ 


103 


145 


83 


95/ 


265 


105 



132 



0.207 










0.75 




1.78 


| 




30° 


22° 


26° 


0 

33 


0 

30 


30° 


43° 




61° 


o 

to 

tO 


39° 


20° 


10° 


0° 


Y-> 

O 

O 


20° 






88.2 










3360 


1 








0.509 


1 




1 



22 ° 

15° 



0.461 



27° 



- J 3 ?- 



t 


ODEN 


PETRI 

VELIXIJ 


SAMPO 1 


SADREL 


SISU 


SOIA 

MARU 


STADT 

reval 


STALIN 

CLASS 


STEPHAN 

MACAEOW 


V 

• 


1957 


1912 


1898 


1929 


1938 


1932 


1895 


1939 


1917 


• 


SWEDEN 


RUSSIA 
1 1 1 


FINLAND 
- 


CANADA 


FINLAND 


JAPAN 




RUSSIA 


RUSSIA 


' 

1 


261-6 


1 

170-7 


191-6 


200-0 


194-9 


310-0 


147-4 


335-0 


236-0 




273-0 


182-1 


202-0 


212-0 


210-6 


327-10 


148-8 


350-0 


248-0 


1 


63-9 


50-10 


43-0 


42-0 


47-6 


46-6 


38-10 


76-0 


57-0 


, 


62.5 


48-7 


42-0 




46-6 




37-1 


74-6 




1 


31-3 


27-5 


21-9 


21-0 


18-6 


30-0 


18-8 


30-0 


30-4 


It 


21-11 


19-1 


18-3 


14-0 


16-1 


15-2 


14-0 


26-3 


22-0 




5020 


1610 


1850 


1860 


2000 


3770 


865 


9300 


4570 


V 


16/17 


14.4 


10/13 

l 


15 


16 


17 


11.5 


16.5 


15 


1 


0.483 


0.378 


0.441 




0.482 




0.396 


0.520 






0.592 


0.533 












0.610 






0.817 


0.708 




1 








0.850 






0.706 


0.723 






1 






0,720 








1 

I j I I 

I 



' JSS- 



1 

1 


| 

Steam 


Steam 


Steam 


Diesel 

Eleotric 


Steam 


Steam 


Steam 


Steam 


2 


1 


1 


2 


2 




1 


3 


2 


2 


1 


1 


0 


1 




0 


0 


1 


10340 


3200 


3000 


3600 


4005 

r 


5500 


1600 


10050 


5550 








11-6 

♦ 


- 


12-6 














11-9 




17-8 








I 

| 


130/ 

143 




155 


140 


112 




125 


90/105 


I 

23° 


30° 


22° 


23° 


24° 




* 

30° 


30° 


18° 


65° 


1 




61° 








69° 


20° 


20° 


I 

20° 

i 




10° 
















I 


0.408 




1 

i 


1 

1 1 
1 




1 
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• 


STORE 

BJORN 


STORIS 


SUNDEW 


SUUR-i 

TOOL 


TARMO 


THULE 


PROPOSED 


TROUVOR 


VACTIONLAND 


1931 


1942 


1944 ! 


1914 


1907 


1953 




1896 


1952 


DEMARK 


u.s. 


U.S. 




FINLAND 


SWEDEN 


RUSSIA 




U.S. 




- 180-5 


220-0 


170-0 


236-4 


210-5 


187-0 

♦ 


363-5 


161-5 


348-0 


196-10 


230-0 


180-0 


247-4 


219-6 


204-0 


390-0 


166-0 


362-0 


• 

49-2 


43-0 


37-0 


57-1 


47-0 


52-9 


80-5 


40-3 


73-6 


46-11 


41-0 




56-1 


46-5 


49-11 


77-1 


39-0 


68-0 


• 23-11 


19-3 


17-3 


22-10 




27-0 


46-0 


22-3 


25-4 


18-2 


14-0 


12-8 


18-8 


18-2 


17-6 


31-2 


18-9 


16-6 


2500 


1916 


1027 


3562 


2300 


2100 


12800 


1450 


6740 


. 13.5 


14 


14.8 


12/14 


13 




19 


9.J/13 


13/16 


0.533 


0.478 




0.505 


0.454 


0.462 


0.507 


0.430 


0.603 


1, 0.628 


0.572 








0.570 


0.600 




0.706 




0.848 


0.835 








0.810 


0.850 




0.854 


i 


0.734 


0.761 










0.650 




0.808. 


* 




15.64 




















20.08 




















322 




















4.44 


3.16 












! 


;* 1 


\ 

. 


3.58 

aft 












i 


i 

i 






2.33 

aft 


i 


I 






1 

i 

1 










- 3 */ ° - 



1 

Steam 


Diesel 

Electric 


Diesel | 
Electric 


Steam 


1 

Steam 


Diesel 

Electric 


Diesel 

Electric 


| 

Steam 


Diesel 

Electric 


2 


1 


1 


2 


1 


2 


3 


1 




1 


0 


0 


1 


1 


1 


0 


0 




5400 


1800 


1675 


4500 


3500 


5040 


22000 


2000 


9440 




10-6 


9-0 
















5 


5 
















7-4 


6-4 














125 


160/ 

195 


205 


100 


89/ 

102 


145 


330 


110 


162 




0.236 


0.236 


















0.045 















24° 


30° 


18° 


20° 


23° 


26° 


25° 










64° 


66° 




o 

r-i 


' 20° 






ro 

o 

o 


18° 





- SV'' 



• 


VO IMA 


VOIMA II 


WALTER 1 
FOSTER 


WILLIAM 

PARSON 


WIND 
CLASS 1 


YMER 


, PROPOSED 


PROPOSED 


ALATNA 


♦ 


1924 


1953 


1954 


1955 


1944 


1932 




1 


1957 


. 


FINLAND 


FINLAND 


CANADA 


CANADA 


U.S. 


SWEDEN 


U.S. 


FINLAND 


U.S. 



' \ 



200-10 


254-4 




325-0 


250-0 


246-1 


340-0 


224-0 


290-0 


210-7 


274-0 


229-2 


351-0 


269-0 


257-11 


361-0 


243-6 


302-0 

% 


46-7 


63-7 


42-6 


68-0 


63-9 


63-4 


75-6 


57-4 


60-11 


45-11 


61-5 




67-0 


62-0 


61-0 


74-0 


54-10 


60-0 


18-9 


31-2 


19-6 


34-9 


37-9 


33-0 




28-10 


26-6 


16-9 


21-0 


16-7 


18-5 


25-9 


21-0 


29-0 


19-0 


19-0 


2070 


4415 


2850 


7108 


5300 


3465 


10500 


3280 


5700 


14 




13.5 


15 


16.5 14/16 


18 


16/17 


12 


0.469 

1 


0.4B5 


0.65 


0.615 


0.465 


0.476 


0.50 


0.489 


0.606 




0.594 






0.618 


0.593 




1 

0.601 


i 

0.654 




0.817 






0.752 


0.801 




0.814 


0.926 




0.701 






0.724 


0.725 


0.80 

I 


0.725 


0.775 


r I 








15.4 




1 




10.8 


; 

1 








22.63 








16.7 










31.04 








26.2 










224 




360 




314 










8.41 








7.5 










15.64 








15.4 










1.0 




5.0 




4.0 










aft 




aft 




eft 










1.22 








0.5 








i 


aft 








fwd 
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Steam 


Diesel 

Electric 


1 

Steam 


Diesel | 
Electric 


Diesel 

Electric 


Diesel Nuclear 
Electric Turbo-Elec 


Diesel 

Electrio 


Diesel 

Eleotric 


1 


2 


2 


2 


2 


2 


2 




2 


1 


2 


0 


1 


0 


1 


0 




0 


3700 


10500 


2500 


1 

13333 


10000 


9000 


30000 


8800 


2530 
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